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AIR/AIR  EJECTOR 

By 
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Chairman:   Vernon  P.  Roan 

Major  Department:   Mechanical  Engineering 

A  detailed  experimental  investigation  of  a  two-dimensional,  Mach  1.8  air-primary, 
Mach  0.3  air-secondary  ejector  at  high  Reynolds  number  has  been  performed,  from  which 
a  nonintrusive  method  for  whole-field  visualization  using  turbulent  wall-pressure  has  been 
developed. 

The  experiments  were  conducted  using  mean  and  time-accurate  wall  pressure 
measurements,  impact-pressure  measurements  using  a  traversing  probe,  and  Schlieren  and 
shadowgraph  visualization  techniques.  The  time-accurate  pressure  measurements  were 
recorded  using  a  sealed  Kulite  miniature  pressure  transducer  with  a  0.7  mm  diameter 
sensing  diaphragm.   For  all  except  the  optical  methods,  measurements  were  taken  from 


XX 


the  initial  flow  interface  to  about  13  hydraulic  tube-diameters  downstream  in  the  constant- 
area  mixing  section. 

From  the  mean  measurements,  values  of  stagnation  pressure,  density,  velocity, 
static  pressure,  Mach  number,  and  dynamic  pressure  were  developed  and  are  presented. 

Using  the  time-accurate  pressure  measurements,  a  color  contour  plot  of  the  rms 
pressure  was  developed  that  definitively  shows  the  regions  of  the  flow  in  agreement  with 
the  other  measurements.  Additionally,  probability  density  functions,  skewness,  and 
kurtosis  were  calculated.  Peak  values  of  skewness  (S)  and  kurtosis  (K)  on  the  centerline 
at  about  2.5  hydraulic  diameters  are  S=1.85  and  K=11.5.  The  inlet  rms  pressure  values, 
normalized  by  freestream  dynamic  pressure  for  the  primary  (-0.001),  were  found  to  be 
in  fair  agreement  with  previous  experimental  values;  however,  those  in  the  secondary 
were  much  higher  (-0.2),  apparently  due  to  the  acoustic  radiation  from  the  primary. 

Fourier  analysis  of  the  time-accurate  pressure  measurements  show  that  the 
autospectra  contain  k',  k'^^\  and  k'"^^  pressure  spectrum  functions  as  predicted  by 
prevailing  theory  for  the  overlap  layer,  turbulence-turbulence  interaction,  and  turbulence- 
mean-shear  interaction,  respectively.  It  is  believed  that  this  is  the  first  experiment  in 
which  the  k'"^^  spectral  slope  has  been  observed,  perhaps  as  a  result  of  the  high  Reynolds 
number. 

Two-point  longitudinal  and  transverse  measurements  are  presented  that  show  the 
development  of  multiple  large  eddies  in  the  flow. 
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CHAPTER  1 
INTRODUCTION 


1  ■  1   Background 

To  meet  the  challenge  of  increasing  the  performance  of  modem  gas  turbines, 
designers  require  knowledge  of  the  smallest  details  and  behavior  of  the  flow.  These 
effects  cannot  be  modelled  using  computational  methods  without  some  knowledge  of  what 
they  are  and  where  they  originate.  Where  they  originate  is  dependent  upon  several  factors 
including  turbulence  behavior,  turbulence/shock  interaction,  the  physical  shape/design  of 
the  enclosure,  and  the  mean  dynamics  of  the  flowfield.  The  range  of  details  from  mean 
to  Kolmogorov  flow  scales,  as  well  as  wall  roughness  and  slip  at  the  wall,  cannot  be 
closely  modelled  even  with  some  of  our  best  computational  techniques,  such  as  Direct 
Numerical  Simulation  (DNS).  Turbulence  alone  can  be  modelled  using  DNS,  but  DNS 
is  limited  by  machine  speed  and  memory  so  that  current  DNS  flows  have,  at  best,  a 
Reynolds  number  of  the  order  of  1,000  and  a  spatial  resolution  that  can  resolve  details 
up  to  about  1000  wall  units,  far  short  of  what  is  required  for  complete  modelling  of  a 
complex  flow. 

Experimental  noise  analysis  techniques  using  piezo-electric  pressure  transducers 
(e.g.,  Willmarth  1958b,  Kistler  and  Chen  1963)  or  microphones  (Schewe  1983)  have  been 
used  for  about  40  years  to  measure  turbulent  wall-pressure,  sparked  by  the  need  to 
investigate  and  control  the  effect  of  pressure  on  aircraft  fiiselages.  Throughout  the  1950s 
and  1960s,  most  fiindamental  experimental  research  focussed  on  turbulent  wall-pressure 
measurements  under  a  thick  flat-plate  boundary-layer  (e.g.,  Willmarth  1956).  This  led  to 
the  familiar  result  of  the  turbulent  energy  cascade  from  large  eddies  (integral  scale)  to 
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small  eddies  (dissipation  scale)  via  the  inertial  subrange,  which  was  in  accord  with 
Kolmogorov's  first  and  second  hypotheses.  Following  the  end  of  the  Apollo  and 
supersonic  transport  (SST)  programs  of  the  1960s,  experimental  research  in  this  area  has 
focussed  on  fundamental  understanding  of  the  small  scale  disturbances  in  the  flow  near 
the  wall. 

An  area  of  current  interest  that  can  be  used  to  further  understand  some  of  the 
fundamental  details  associated  with  internal  mixing  flows  is  to  study  the  ejector  shear- 
layer,  which  comprises  a  supersonic  primary  jet  and  a  subsonic  secondary  flow. 
Presently,  the  model  of  a  developing  ejector-flowfield  is  typically  divided  into  three 
regions  preceded  by  uniform  primary  and  secondary  flows,  and  followed  by  a  region  of 
fully  mixed  flow  (see  Fig.  1.1).  Modelling  the  ejector  using  a  computational  process, 
which  will  not  be  addressed  in  this  work,  involves  a  marching  scheme  that  starts  at  the 
point  of  initial  contact  (start  of  region  1). 

Region  1  begins  with  this  initial  contact  and  ends  where  the  potential  core  of  the 
primary  ends.  Region  2  also  begins  with  the  initial  contact  and  ends  where  the  secondary 
shear  layer  reaches  the  wall.  Region  3  starts  at  the  end  of  region  1  or  2  (whichever  is 
longer)  and  continues  until  the  streams  are  fully  mixed. 

Recent  experimental  work  at  the  University  of  Florida  (Petersen  et  al.  1992,  Roan 
et  al.  1992,  1993,  Dufflocq  et  al.  1992,  1993,  Benjamin  et  al.  1993)  has  focused  on 
traversing  an  impact  probe  across  both  axisymmetric  and  two-dimensional  planar  flows 
at  several  axial  positions,  which  has  enabled  some  understanding  of  the  flow  inasmuch 
as  pumping  performance  was  measured  and  compared  to  theory,  growth  rates  were 
estimated,  and  the  lengths  of  regions  1,  2  and  3  were  measured  or  inferred  when  possible. 
This  experimental  work  has  led  to  the  development  of  a  semi-empirical  model  for 
predicting  the  lengths  of  regions  1  and  2  (Benjamin  and  Roan  1993).  In  the  axisymmetric 
case,  the  experiments  identified  large  coherent  structures  in  the  primary  core  that  probably 


3 

result  from  vortex/shock-wave  interaction,  whereas  such  structures  are  not  apparent  in  the 
two-dimensional  planar  case.  Figure  1 .2  shows  a  comparison  of  an  axisymmetric  and  a 
two-dimensional  planar  case  at  similar  flow  conditions.  The  result  of  the  vortex/shock 
interactions  is  seen  in  region  1  of  the  axisymmetric  case. 

Most  experimental  research  of  high  speed  internal  flows  has  relied  mainly  on 
techniques  such  as  laser  doppler  velocimetry  or  insertion  of  some  type  of  pressure  probe 
into  the  flow.  The  current  research  was  undertaken  using  a  different  approach  to  provide 
additional  detailed  measurements  necessary  to  better  understand  and  enhance  the  ability 
to  measure  and  predict  the  behavior  of  ejector  flowfields.  The  research  was  also  designed 
to  provide  more  insight  into  the  structure  of  this  type  of  internal  flow  and,  hence,  enable 
designers  to  better  optimize  their  designs. 

This  research  investigates  compressible  shear-flow  regions  that  are  sandwiched 
between  two  essentially  inviscid  streams  of  supersonic  and  subsonic  fluid,  all  of  which 
are  surrounded  by  boundary  layers.  Far  downstream  where  the  shear  layers  have  reached 
the  walls,  the  flowfield  comprises  only  shearing  fluid  in  the  freestream,  which  is  again 
surrounded  by  boundary  layers.  There  is  no  reason  to  expect  that,  in  terms  of  the 
turbulent-pressure  field,  the  inviscid  flow  regions  behave  any  differently  from  previously 
reported  subsonic  and  supersonic  turbulent  boundary-layer  flows,  although  jet-noise 
contamination  may  be  present.  The  behavior  of  the  shearing  regions  is  undetermined  and 
results  from  this  research  will  characterize  them.  On  the  small  scales  in  discrete  ranges 
of  the  universal  equilibrium  range,  local  isotropy  should  generally  still  apply  because 
compressibility  simply  acts  as  a  source  of  dissipation  in  addition  to  the  viscous  stresses 
(Hinze  1975,  p.  312).  Also,  if  there  is  significant  acoustic  radiation  emitted  or  absorbed, 
a  change  in  the  energy  transfer  rate  will  occur  that  will  cause  a  change  in  the  spectral 
slope.  Hence,  a  knowledge  of  the  physics  or  behavior  of  the  following  is  required  to  lay 
the  foundations  for  ejector  flow: 
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(1)  the  fundamental  premises  of  turbulent  flow, 

(2)  the  turbulent  boundary-layer  in  subsonic  and  supersonic  flows  under  the 
influence  of  zero-,  adverse-,  and  favorable-pressure  gradients, 

(3)  mean  freestream-pressure  and  velocity  fields, 

(4)  ejector  shear-layer  growth  rates,  and 

(5)  acoustic  phenomena  of  jets. 

The  effect  of  pressure  transducer  size  on  resolution  of  all  the  pressure  scales  in 
the  flow  is  a  concern  that  has  been  studied  to  a  fair  degree  and  is  also  outlined. 

1.2   Ejectors 

Ejector  systems  have  been  investigated  for  many  years  due  to  the  interest  in  their 
fluid-pumping  ability.  These  devices  are  of  particular  interest  to  the  aircraft  industry, 
where  they  have  been  studied  and  employed  for  their  enhancement  of  thrust  augmentation, 
lift  augmentation  and  engine-noise  reduction  (Braden  et  al.  1982).  Other  possible 
applications  of  ejectors  include  jet  pump  compression,  extraction  of  secondary  fluids, 
ventilation  and  air  conditioning,  innovative  thermal  cycles,  and  high  energy  chemical 
lasers  (Power  1994). 

The  pumping  ability  of  ejectors  is  primarily  the  result  of  momentum  transfer  from 
the  high-velocity  primary  fluid  to  the  slower  secondary  fluid.  Ultimately,  the  overall 
performance  of  these  devices  is  influenced  by  a  number  of  physical  processes  that 
individually  and/or  collectively  govern  the  development  of  the  associated  flowfields. 
Among  the  processes  of  interest  are  (a)  details  of  mixing  of  primary  and  secondary 
streams,  (b)  viscous  dissipation  due  to  mixing,  compression  waves,  and  solid  boundaries, 
and  (c)  the  effects  of  geometry.  Maximizing  ejector  efficiency  translates  into  enhancing 
the  mechanisms  that  result  in  efficient  mass  and  momentum  transfer  and  inhibiting  those 
that  have  a  detrimental  effect. 
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An  ideal  ejector  is  one  yielding  ideal  primary-to-mixed  out  momentum  transfer  as 
determined  from  a  control  volume  analysis.  In  a  real  ejector,  the  closest  to  ideal  that  can 
be  obtained  is  a  configuration  where  the  shroud  of  the  ejector  is  long  enough  to  allow  for 
near  completion  of  the  mixing  process  and  yet  short  enough  to  minimize  wall  friction 
losses.  (Additionally,  there  are  area  ratio  trade-offs.)  These  are  conflicting  ideals  that 
have  been  well  substantiated  by  ejector  studies.  Therefore,  ejector  design  compels  trade- 
offs where  the  goal  is  to  find  the  optimum  balance  between  these  opposing  factors  for 
given  requirements  and  input  parameters. 

Traditionally,  theoretical  ejectors  studies  have  been  aimed  at  comparing  predicted 
overall  performance — usually  in  terms  of  thrust  augmentation,  overall  pressure  ratio,  or 
entrainment  ratio — ^with  results  obtained  experimentally  (Braden  et  al.  1982,  Petersen  et 
al.  1992).  In  many  cases,  performance  is  predicted  using  one-dimensional  conservation 
equations  to  find  flow  properties  at  the  "fully  mixed"  exit  plane  of  the  ejector.  Such 
models  show  the  general  trends  that  are  expected  under  ideal  flow  conditions  (uniform 
flow  and  no  losses).  However,  they  generally  predict  performance  superior  to  those  of 
real  ejectors,  and,  hence,  their  value  as  design  tools  is  limited. 

Similarly,  semi-empirical  models  have  been  used  to  predict  real  ejector  flowfields 
in  which  undesired  dissipation  effects  are  unavoidable.  Here,  factors  such  as  wall 
frictional  losses,  nonuniform  inlet  and  exit  flow,  imperfect  mixing  of  the  two  streams,  and 
primary  nozzle  efficiencies,  for  example,  are  taken  into  account  using  best  estimates.  The 
approach  in  this  case  has  usually  been  to  adjust  parameters  within  the  model  ("calibrate") 
so  that  predicted  flowfields  match  the  experimental  results.  These  models  are  only  valid 
for  the  range  of  calibrated  ejector  geometries  and  flow  parameters,  and,  therefore,  their 
potential  as  design  tools  is  also  limited. 

Progress  in  ejector  technology  largely  depends  on  being  able  to  predict  ejector 
flows  accurately.  This  ability  can  only  be  achieved  through  a  better  understanding  of  the 
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Figure  1.1        Definition  of  ejector  regions  and  lengths. 
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Figure  1.2 
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(a)  axisymmetric  configuration  M^-U-O^, 

(b)  two-dimensional  planar  configuration,  M^  0.41. 
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LITERATURE  SURVEY 

.  ,.„„. .  r„ — ...  *7- "::::::  -„.. ,-. 

.0...  ...  .„  ..0,  aave.e,  and  f—  P^ --';;;^^^^^^^^^^^  ,,  ^,„, 

e.ene.  — .s.e  ....  .  ..  .ea  ...  .a.  ^n  .a       -e  .  ^^^ 

iw,  r,.!neciallv  supersonic)  boundary  layer  now,  ii 
studies  n  compressible  (especially  suH  excellent 

.  ,  ,o  the  inherent  difficulty  of  the  measurements.     Excellent 
limited  in  scope  due  to  the 

comprehensive  reviews  of  turbulent  wall-pressures  in  incompressible 
.ompiled  by  W.llmarth  (.75)  and  Be—  (1^0,  ^^_^^^^^^^ 

-  -Ksround  on  turbulence,  a"-—  ^J;       ^„.,.  ,e  motivation 
hypotheses  (Kolmogorov  1941a,  b,  hereinafter  referred  to  as  K41),  P 
,    the  expenmenta,  methods  and  results  presented  m  the  present  research. 


10 


2j_rurbiUenceThe^^ 
p.ese„.e.,  foUo^ed  by  a  d.css.on  of ..  app.,cat,o„  to  specm.™  ^nct,o.. 

-^^r^^r^es.  .e  ......= -.,0.  one.  .^^^^^^^^^ 

..  Ju..p.   — ce  .s  *  .™o.„eo.  ...  *e  spa.Ue— 0, 


mean 


turbulence  quantities  are  zero,  e.g. 


-mA^.0  (2,) 


The  nrst  hypothesis  can  be  summanzed  as  follows^    f <"•  ' 
ine  iiisi    JK"  ,         .       .u.  ua-lmknce  is  slalislicaily 

.uriulence  ,>,ere  e.is,s  a  range  omh  .a^enumbers  .here  ,he  iMer^e 

the  kmematic  viscosity,  V.   The  state  oj  eqm  ,  ^f  , eternal 

,„.We«c.  c.«  on,y  U  a  resuU  o,,,.  e^c.  of.^  paran.e.ers  s  ar^  v. 


11 

To  put  this  hypothesis  into  mathematics  (Sreenivasan  1991),  the  hypothesis  can 
be  restated  such  that  the  multivariate  probability  distributions  of  the  velocity  difference 
in  a  direction  /,  Au.=u.(x)-u.(x+f),  are  functions  only  of  r=|f|,  s,  and  v.  Specifically, 
the  nth-order  structure  functions  have  the  form 

^i"-(errarM)  (2.2) 

Here,  T]={y^ley'^  is  the  Kohnogorov  length  scale  and  the  functions^  are  universal;  ti  is 
on  the  order  of  magnitude  of  the  energy  dissipation  scales. 

The  second  K41  hypothesis  states  that,  at  very  high  Reynolds  numbers,  there  is 
a  range  of  scales  T]<r<L  ,  called  the  inertial  subrange,  in  which  the  structure  Junctions 
become  independent  of  viscosity;  that  is,  /^(oo)=X^,  where  K„  are  universal  constants 
independent  of  the  flow 


The  K41  hypotheses  arise  fundamentally  from  dimensional  arguments,  and  their 
validity  rests  on  experimental  support  (Sreenivasan  1991). 

The  idea  of  inter mittency  was  first  pointed  out  by  Landau  in  1944  (Landau  and 
Lifshitz  1987)  and  arises  due  to  spatial  variations  off  over  physical  volumes  larger  than 
the  inertial  subrange,  i.e.,  the  energy  cascade  process  in  the  inertial  subrange  undergoes 
fluctuations  caused  by  larger  than  inertial-subrange  scales.  Therefore,  intermittency 
causes  nonuniversality  of  the  constants,  K„,  because  the  large  scales  are  flow-specific. 
Recent  experimental  measurements  of  high-order  velocity  structure  functions  (Anselmet 
et  al.  1984)  indicate  that  the  scaling  function  exponent,  n,  shows  strong  deviation  (i.e., 
non-Gaussian  statistics)  from  the  K41  theory  for  structure  functions  higher  than  second- 
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order.  There  is  currently  much  discussion  and  controversy  about  the  fundamental  cause 
for  intermittency,  and  the  reader  is  referred  to  Sreenivasan  (1991)  and  Anselmet  et  al. 
(1984)  for  more  detailed  discussion. 

2.1.2   Inertial  Spectrum  Functions 

Using  Kolmogorov's  second  hypothesis,  it  can  be  shown  from  dimensional  analysis 
(e.g.,  George  et  al.  1984)  that  the  velocity  spectrum  function  in  the  inertial  subrange  is 
given  by 


1    -1  (2.4) 

E{k)=ae'k  ' 


where  E(^)  is  defined  by 


E{k)=UfF,,{k)d<^) 

(2.5) 


where  the  integral  is  over  spherical  shells  a  of  radius  k.   The  cross-spectral  densities  of 
the  velocity  moments  are 


n#")^;r^///-/^'Xr)''=^- 


P'r  (2.6) 


with  the  velocity  moments  defined  as 

5,.r)=«.(f)u.(f.r-)  (2.7) 
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Tennekes  and  Lumley  (1987)  recommend  that  experimental  data  indicate  a=1.5, 
approximately.  Figure  2.1  shows  the  form  of  the  velocity  spectrum  function  over  all 
wavenumbers. 

George  et  al.  (1984)  derived  models  for  turbulent  pressure  fluctuations  by  directly 
Fourier  transforming  the  integral  solution  to  the  incompressible  Poisson  equation  for  a 
constant-mean-shear-flow  (the  integral  solution  is  presented  in  the  next  section).  They 
found  that  the  spectrum  function 

EJhy-ffFj1^)doif)  (2.8) 

is  composed  of  three  parts  as  follows: 

E^iJc)=irJk)*7Tjk)^7r,{k)  (2.9) 

where  tt^  and  v^  represent  the  second-  and  third-moment  turbulence-mean-shear 
interaction  spectrum  functions  respectively,  and  rr,  represents  the  turbulence-turbulence 
interaction  spectrum  function.  The  inertial  subrange  for  each  of  these  spectrum  functions 
must  be  determined  by  only  the  parameters  k,  e,  and  the  mean-shear  velocity-gradient  K 
By  dimensional  analysis,  they  arrive  at 


±'rrJk)=a,K'eh''^ 
P  (2.10) 
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\:r^,{k)=a,Ksk-' 
P  (2.11) 


1  -    -- 

J-'Tr{k)=a  B'k  ' 

2     P  P 

P  (2.12) 

where  c^,  aj,  and  o.^  are  constants.  These  constants  are  found  by  assuming  that  the 
turbulence  is  isotropic,  although  the  turbulence-mean-shear  interaction  is  not.  The  values 
of  constants  calculated  by  George  et  al.  (1984)  for  the  one-dimensional  spectrum  (which 
can  be  measured)  are 

)32=0.307a  (213) 

i33=0  (2.14) 

/3,=0.283a2  (215) 

From  their  calculations  George  et  al.  (1984)  suggest  that  the  shear  contribution  is 
impossible  to  isolate  from  spectral  data,  even  in  flows  with  moderate-to-high  shear  rates 
relative  to  the  turbulence  shear-rate  u/l,  and  conclude  that  this  is  the  reason  that  such  a 
range  has  not  been  previously  observed.  They  also  note  that  there  is  no  direct  dissipation 
of  pressure  fluctuations. 

Bradshaw  (1967b)  argues  by  dimensional  reasoning,  or  by  direct  substitution  into 
the  Poisson  equation,  that  the  one-dimensional  spectral  density  is  proportional  to  k'  in  the 
overlap  region  of  the  turbulent  boundary  layer.  Panton  and  Linebarger  (1974)  calculated 
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wall  pressure  spectra  using  the  Law  of  the  Wall  and  Coles'  wake  function,  which  also 
revealed  an  overlap  region  of  A:"'.  Following  Panton  and  Linebarger  (1974),  P  denotes  the 
spectrum  in  outer  variables  and  T  the  spectrum  in  iimer  variables,  n  is  Coles'  pressure 
gradient  parameter.   Then 

P(ki6;  n)  =  <j)(k,)/Xw'  6,        !P(k,v/u^  =  <t)(kOu,  Ix^'  v,        Re  =  u,6/v  (2. 16) 

There  is  an  intermediate  wavenumber  where  these  expressions  are  equal:  P=!P/Re. 
Additionally,  if  there  is  a  range  of  wavenumbers  where  the  expressions  are  equal,  then 
the  equality  may  be  differentiated  for  this  range  such  that 

(k,6)2  P'(k,6;  O)  =  (kiv/u^'  !P'(k,v/u^  =  const.  (2. 1 7) 

For  this  expression  to  be  true,  P'  must  be  independent  of  n  in  the  overlap  region,  and  the 
spectra  have  the  form 

P(k,6)^  =  A(k,6)-' ,     !P'(kiv/uJ  =  A(k,v/u^->  (2.18) 

This  gives  a  straight  line  slope  of -1  on  a  log-log  scale.  Panton  and  Linebarger  (1974) 
conclude,  along  with  Bradshaw  (1967b),  that  the  overlap  region  is  beyond  the  reach  (too 
high  frequency)  of  most  measurements. 

Recent  measurements  by  Farabee  and  Casarella  (1991)  show  that  in  the  very  low 
wavenumber  region,  the  spectrum  follows  a  k'  behavior,  as  suggested  by  Bradshaw 
(1967a). 
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2.2   Turbulent-Pressure  Investigations  in  Boundary-Layer  Flows 
The  flows  in  ejectors  have  regions  of  zero,  adverse  and  favorable  pressure- 
gradients  that  will  affect  the  wall-pressure  measurements.     Hence,  presentation  of 
pertinent  literature  for  such  pressure-gradient  flows,  as  well  as  the  effects  of  rough  and 
smooth  walls,  is  covered  in  the  following  three  subsections. 

2.2.1    Incompressible.  Zero  Pressure-Gradient  Flows 

The  incompressible  Navier-Stokes  equation  may  be  given  by 


du.    duu.      I  dp         d^ 

i+ LJ.  =  - il  +  v Uj 


and  taking  the  divergence  of  Eq.  (2.19)  gives  the  relationship  between  the  pressure  and 
velocity  fields  (Lilley  and  Hodgson,  1960): 


^^i  (2.20) 

where 

q= '- 

dx^dXj  (2.21) 

The  right-hand  side  of  Eq.  (2.21)  can  be  manipulated  (Thomas  and  Bull,  1983)  to  give 


dU.  duj'     d^  I  /  /  -r-,\ 


^      dXj  dx.     dx.dx}  '  '      '  "  (2.22) 
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These  equations  show  that  the  pressure  field  is  related  to  the  velocity  field  through  the 
turbulence-mean-shear  and  turbulence-turbulence  shear  stress  interactions. 

The  solution  for  the  fluctuating  pressure  on  the  wall  p^  using  the  integral  solution 
for  the  Poisson  equation  in  the  half-plane,  and  neglecting  surface  integrals  (Lilley  and 
Hodgson,  1960)  can  be  written 

2V>o      T-^^l  (2.23) 

This  indicates  that  sources  over  a  region  of  flow  (theoretically,  a  semi-infinite  region)  will 
contribute  to  the  wall  pressure  fluctuations  and  that  the  contributions  from  various  source 
regions  will  fall  off  rapidly  with  their  distance  from  the  point  under  consideration 
(Thomas  and  Bull,  1983). 

Kraichnan  (1956a,  b)  was  the  first  to  report  theoretical  estimates  for  the  mean- 
square  wall-pressure  and  spectra.  He  overcame  the  difficulty  of  having  to  treat  the 
turbulent  boundary  layer  as  anisotropic  by  assuming  that  the  turbulent  flow  is 
homogeneous  in  planes  parallel  to  the  wall.  The  source  terms  for  the  turbulence-mean- 
shear  interactions  were  estimated  based  on  a  rough  approximation  (-50%  accuracy)  to  the 
mean  shear  and  Kraichnan's  (1956a)  results  for  the  Fourier-transformed  velocity  and 
pressure  fields  for  homogeneous  anisotropic  turbulence.  The  prediction  of  the  ratio  of 
root-mean-square  (rms)  wall-pressure  fluctuations  to  mean  wall  shear  stress  was  of  order 
six,  much  more  accurate  than  would  be  expected.  He  also  found  that  for  the 
homogeneous  and  anisotropic  cases,  the  rms  turbulent  pressure  for  the  turbulence-mean- 
shear  interactions  was  8.25  times  larger  than  that  of  the  turbulence-turbulence  interactions. 

Incompressible  measurements  by  Willmarth  and  Wooldridge  (1962)  under  a  thick 
boundary  layer  showed  that  the  rms  wall  pressure  was  2.19  times  the  wall  shear  stress. 
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Static  pressure  fluctuations  over  a  flat  boundary  in  the  atmospheric  boundary  layer  by 
Elliott  (1972)  were  found  to  produce  turbulent  pressure  values  of  about  2.6  times  the 
mean  stress  at  heights  up  to  6  m.  Hinze  (1975)  recommends  that  a  reasonable  average 
from  all  incompressible  experimental  data  available  is  /7^«3t^.  Willmarth  (1975) 
concluded  that  the  ratio  of  turbulent  wall-pressure  to  freestream  dynamic  pressure  has  a 
value  of  0.0035  for  incompressible  flow. 

Corcos  (1964)  made  some  calculations  based  on  pressure  and  velocity  experiments 
and  found  that  the  inner  part  of  the  Law  of  the  Wall  region  (y^slOO)  seemed  to  be 
substantially  free  of  pressure  sources  (except  perhaps  at  very  high  frequencies),  and  within 
that  region  (a)  the  pressure  can  be  given  in  terms  of  its  boundary  value,  and  (b)  the  local 
velocity  field  is  dependent  upon  but  unable  to  affect  appreciably  the  turbulent  pressures. 

Using  the  method  of  space-time  correlation,  Willmarth  (1958a)  discovered  that  the 
random  pressure  fluctuations  under  a  flat-plate  boundary  layer  were  convected  with  a 
speed  of  approximately  0.8f/„.  More  detailed  experiments  (Willmarth  and  Wooldridge 
1962,  Bull  1963)  showed  that  the  convection  velocity  varies  with  streamwise  spatial 
separation  of  the  measuring  stations  and  that  for  small  spatial  separation  the  convection 
velocity  is  low,  0.56C/„,  but  increases  to  0.83 f/^  for  very  large  spatial  separation.  The 
increase  in  convection  velocity  with  streamwise  separation  of  measuring  points  was 
attributed  to  the  more  rapid  decay  of  the  smaller  pressure  producing  eddies.  Wills  (1970) 
measured  the  wave-number/phase  velocity  spectrum  of  the  wall-pressure,  which  showed 
the  above  relation  for  different  wavenumbers  (eddy  sizes).  Willmarth  and  Wooldridge 
(1962)  also  concluded  that  a  pressure-producing  eddy  of  large  or  small  wavelength  X 
decays  and  vanishes  after  travelling  a  distance  of  approximately  6X. 

The  effect  of  wall  roughness  was  investigated  by  Blake  (1970),  who  found  that 
there  was  no  effect  on  the  ratio  of  rms  wall  pressure  to  wall  shear  stress.   However,  the 
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magnitude  of  the  rough-wall  rms  pressure  increased  by  at  least  a  factor  of  two  over  the 
smooth  wall  value. 

From  an  analysis  of  several  experimental  investigations,  Thomas  and  Bull  (1983) 
proposed  that  the  dominant  contribution  to  the  rms  pressure  fluctuations  comes  from  the 
source  term  associated  with  the  turbulence-mean-shear  in  agreement  with  Kraichnan  and 
hence,  for  two-dimensional  flow,  Eq.  (2.23)  becomes 


27r^J„  dy  dx  [f-x^l  (2.24) 


Schewe  (1983)  found  that  the  skewness  of  the  turbulent  pressure  has  a  value  of 
-0.2  and  a  kurtosis  of  4.9,  compared  to  values  of  0  and  3,  respectively,  for  a  Gaussian 
distribution.  Analysis  of  the  time  records  showed  that  the  difference  from  the  Gaussian 
was  caused  by  wavetrains  or  pulses  that  convect  at  0.53U„,  and  can  therefore  be  attributed 
to  pressure  structures  in  the  buffer  layer  at  12;Sy^:s21  based  on  the  measured  velocity 
profile.  The  measurements  also  confirmed  that  there  were  more  negative  pulses  than 
positive  ones  as  indicated  by  the  negative  skewness.  The  characteristic  wavelength  of  the 
wavetrains  based  on  their  fi-equency  of  occurrence  and  their  convective  velocity  was 
X*='kuJ\=l45.  Hence,  a  transducer  would  have  to  be  smaller  than  ^'^=145  to  begin  to 
resolve  the  wavetrains.  Additionally,  Schewe  calculated  from  the  probability  density 
function  of  events  pth>3|pn„J  that  they  contribute  40%  to  the  rms  pressure,  although  the 
events  only  occur  for  1%  of  the  time.  Karangelen  et  al.  (1993)  found  that  the  same 
threshold  was  exceeded  for  5%  of  the  time  and  contributed  49%  to  the  rms  pressure. 
Hence  it  can  be  seen  that  the  pressure  pulses  can  vary  significantly  between  different 
flows.  The  measured  time  between  burst  events  for  channel  flow  (Tiederman  1990)  was 
found  to  be  90  scaled  on  inner  variables  and  was  independent  of  Reynolds  number.  For 
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boundary  layer  flow  the  time  between  events  is  closer  to  300  (Blackwelder  and 
Haritonidis  1983).  Using  these  values  as  lower  and  upper  bounds,  Karangelen  et  al. 
(1993)  conclude  that  the  average  time  between  bursts  corresponds  to  the  pressure  events 
range  of  2sp,^/|p,^|s3,  which  suggests  a  direct  correspondence  between  bursts  and  large 
amplitude  wall  pressure  events. 

Dinkelacker  and  Langeheineken  (1983)  studied  the  relation  between  intermittent 
high  positive  wall-pressure  fluctuations  and  velocity  fluctuations  in  a  50  mm  diameter 
pipe  with  Re]3=9600.  The  results  were  conditioned  over  several  thousand  events  and 
clearly  show  that  the  high  wall-pressure  events  are  related  to  a  steep  increase  in  <u'> 
measured  between  y*=5  and  y'^=52.  Preceding  this  event  is  a  comparatively  long  period 
in  which  <u'>  is  below  average  and  followed  by  a  comparatively  long  period  in  which 
<u'>  is  above  average.  From  the  duration  of  the  periods  the  streamwise  length  of  the 
pattern  was  calculated  to  be  twice  the  pipe  diameter  (Ax "^^l 000).  Similarly,  the  <v'> 
measurements  showed  that  the  flow  moves  away  from  the  wall  and  then  towards  the  wall 
either  side  of  the  turbulent-pressure  peak,  but  occurs  over  1/6  the  period  of  the  <u'>  event 
and  to  a  height  of  at  least  y*=100.  The  disturbance  in  the  azimuthal  direction  was  found 
to  have  a  span  of  order  At]"'=50.  It  was  concluded  that  these  events  were  "bursts,"  which 
were  defined  by  other  investigators  to  be  instantaneous  high  levels  of  turbulence  at  about 
y^=15. 

Kobashi  and  Ichijo  (1986)  measured  wall  pressure  fluctuations  in  relation  to  the 
coherent  motions  of  the  turbulent  boundary  layer.  The  low  frequencies  were  found  to  be 
related  to  the  large  scale  motions  of  the  outer  layer  and  prevail  through  and  outside  the 
boundary  layer,  while  the  high  frequencies  were  found  to  be  related  to  the  bursts  in  the 
wall  region.  They  concluded  that  the  large  scale  motions  are  composed  of  periodic  and 
nonperiodic  components  that  are  initiated  by  the  instability  of  the  mean  flow  and  rotate 
in  the  direction  o/the  mean  flow  shear  (cf  Fig.  2.2).    The  bursts  are  vortical  in  nature 
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and  rotate  against  the  mean  flow  shear.  The  origin  of  the  bursts  is  from  the  interaction 
of  the  periodic  large  scale  motions  with  the  wall.  The  burst  is  bound  by  two  pairs  of 
counter-rotating  vortices  that  are  inclined  forward,  and  the  ejection  and  sweep  that 
characterize  these  burst  phenomena  can  be  explained  as  the  flows  induced  by  the  paired 
vortices.  Kobashi  and  Ichijo  (1986)  also  hypothesized  that  the  paired  vortices  are  formed 
from  streamwise  vortex  pairs  that  appear  in  the  large-scale  periodic  motions. 

The  coupling  between  high-amplitude,  positive  wall-pressure  peaks  and  flow 
structure  in  the  near  wall  region  was  studied  by  Johansson  et  al.  (1987).  They 
characterized  wall-pressure  fluctuations  in  a  two-dimensional  boundary  layer  by  making 
simultaneous  measurements  of  high-amplitude,  positive  turbulent  wall-pressure  peaks  with 
u'  and  v'  .  The  results  indicated  that  the  amplitude  of  the  pressure  peak  is  linearly  related 
to  the  amplitude  of  the  turbulent  velocity  peak,  which  in  turn  indicates  that  the  generation 
of  high-amplitude  pressure  peaks  is  predominantly  governed  by  the  turbulence-mean-shear 
interaction.  Also,  these  pressure  peaks  were  found  to  be  caused,  or  the  cause  of,  shear- 
layer  structures  ("bursts")  in  the  buffer  region  (5<>'^<30).  Large  negative  turbulent  wall- 
pressure  peaks  were  found  to  be  primarily  associated  with  sweep-type  motions — flow 
parallel  to  the  wall. 

2.2.2   Compressible,  Zero  Pressure-Crradient  Flow 

Measurements  on  the  wall  under  a  turbulent  Mach  0.6  boundary  layer  by  Serafmi 
(1963)  show  that  the  convection  velocity  varies  for  different  streamwise  transducer 
spacing  in  exactly  the  same  way  as  for  the  incompressible  boundary  layer.  Lilley  (1963) 
showed  that  the  wall  pressure  fluctuations  are  the  result  of  fluctuations  in  both  the 
vorticity  and  sound  modes.  Kistler  and  Chen  (1963)  found  that  the  rms  pressure  on  a 
solid  surface  for  flows  in  the  Mach  number  range  1.33-5.00  was  proportional  to  the  local 
skin  friction.  The  value  of  Pnns/x„  increases  fairly  linearly  from  the  incompressible  value 
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of -3  to  ~5  for  M2;2.  They  also  found  that  the  convection  speed  falls  from  0.8  at  M=1.33 
to  0.6  at  M=5.  The  effect  of  Reynolds  number  was  found  to  be  P,^/<2~-^^«-  •  Th® 
integral  scale  of  the  wall-pressure  fluctuations  changed  from  0.166  at  M=1.33  to  0.00066 
at  M=4.54.  Additionally,  they  found  that  the  peak  value  of  the  correlation  coefficient 
drops  to  one  half  for  a  spatial  separation  of  the  measuring  points  of  about  two  tenths  of 
the  boundary  layer  thickness.  Laufer  (1964)  measured  the  radiated  pressure  intensity  from 
a  supersonic  boundary  layer  and  found  it  to  be  two  orders  of  magnitude  less  than  that 
measured  on  the  wall. 

2.2.3   Incompressible.  Adverse  and  Favorable  Pressure-Gradient  Flows 

Adverse  pressure-gradient  flows  differ  from  zero  and  favorable  pressure-gradient 
flows  in  that  they  are  not  self-preserving,  and,  hence,  the  value  of  every  parameter  is  a 
function  of  local  position.  Burton  (1973)  compared  strong  adverse  and  favorable 
pressure-gradient  flows  with  smooth  and  rough  walls  to  a  zero  pressure-gradient  flow. 

For  a  favorable  gradient,  it  was  found  that  rms  wall  pressure  intensity  varied  in 
proportion  to  mean  wall  shear-stress  with  the  same  value  as  the  zero  pressure-gradient 
case.  The  favorable  gradient  decreased  longitudinal  spatial  decay  rates  and  increased 
convection  velocities,  whereas  roughness  had  the  opposite  effect. 

The  adverse  gradient  slowed  convection  velocities  and  increased  spatial  decay 
rates.  Pressure  statistics  were  found  to  depend  on  local  mean  flow  parameters  and 
upstream  flow  conditions,  but  not  on  wall  roughness.  The  velocity  profiles  were 
logarithmic  over  a  small  region  close  to  the  wall,  but  were  generally  wake-like,  so  that 
outer  variable  scaling  was  appropriate. 

Insight  into  the  suitable  scaling  of  boundary  layers  can  be  seen  from  treatment  of 
the  integrated  momentum  equation  for  two-dimensional  turbulent  flow  (Burton,  1973): 


,^ 
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.8'^-  pilule)  =  -U„  '^dy-?f4^^^—\dy  'dy^s'f^^^^y 


dx       dx^       '      •'^dx  0  0^ 


dx       dy' ) 


dx'  (2.25) 


For  rough  walls  the  terms  on  the  right-hand  side  of  Eq.  (2.25)  must  be  retained, 
whereas  for  smooth  walls  the  terms  on  the  right-hand  side  of  Eq.  (2.25)  are  at  least  three 
orders  of  magnitude  less  than  those  on  the  left-hand  side  and  can  therefore  be  written  as 

^  (2.26) 

For  a  favorable  pressure  gradient  flow,  the  two  terms  on  the  left-hand  side  of  Eq. 
(2.26)  are  of  equal  magnitude,  and  their  ratio  is  a  valid  measure  of  the  magnitude  of  the 
pressure  gradient.  However,  for  adverse  gradient  flow,  the  wall  shear  stress  is  small 
compared  to  the  pressure  gradient  term  and  the  right-hand  side  of  Eq.  (2.26).  Hence, 
scaling  with  inner  wall  variables  is  not  appropriate  for  the  adverse  pressure  gradient  case. 

Burton  also  concludes  that  broadband  spatial  coherence  of  the  wall-pressure 
fluctuations  is  only  moderately  affected  by  imposition  of  a  favorable  pressure  gradient, 
when  scaled  on  outer  variables.  Coherence  is  markedly  improved  in  the  adverse  gradient 
flows,  which  are  dominated  by  disturbances  that  are  large  compared  to  the  displacement 
thickness  and  therefore  decay  more  slowly.  In  only  adverse  gradient  flows,  broadband 
convection  velocities  were  found  to  be  higher  between  oblique  separations  than  between 
longitudinally  separated  points.  This  was  also  the  case  for  narrowband  convection 
velocities.  Phase  velocities  were  only  moderately  affected  by  favorable  gradients,  but  the 
effect  on  phase  velocities  for  adverse  gradient  flow  was  strong  and  complex,  but  always 
lowered  the  phase  velocities  considerably,  the  lowest  values  being  about  0.25U„. 
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2.3   Finite  Transducer  Size  Effects 

The  problem  of  inadequate  spatial  resolution  of  a  pressure  transducer  was 
recognized  by  Willmarth  (1956).  The  finite  size  of  a  transducer-sensing  element  limits 
its  space  resolution  of  a  convecting  pressure  field  associated  with  a  turbulent  flow. 
Consequently,  a  lack  of  resolution  in  space  causes  an  apparent  inability  to  resolve  in  time. 
Corrections  were  first  attempted  by  Corcos  et  al.  (1959),  followed  by  an  improved  method 
by  Corcos  (1963).  This  method  was  based  on  new  experimental  measurements  of  the 
longitudinal  and  lateral  space  correlations  and  predicted  that  the  attenuation  at  high 
frequencies  was  much  greater  than  predicted  from  the  previous  method. 

Willmarth  and  Roos  (1965)  used  several  different-sized  pressure  transducers  and 
extrapolated  the  results  to  obtain  the  vanishingly  small  value  of  rms  pressure.  The  result 
was  p^T„=2.66,  which  would  account  for  Willmarth  and  Wooldridge's  (1962)  results  to 
be  increased  by  approximately  13%.  They  also  concluded  that  Corcos's  (1963)  correction 
can  be  used  at  low  frequencies  but  not  at  high  frequencies. 

Schewe  (1983)  experimentally  showed  that  the  dimensionless  transducer  diameter 
d'^=duyv=19  is  sufficient  to  resolve  the  essential  structures  of  the  turbulent  pressure 
fluctuations.  If  this  criterion  is  assumed  to  hold  for  adverse  gradient  flow,  then  the 
condition  d*=19  is  easier  to  obtain  in  adverse  gradient  flow  due  to  the  reduction  in  t„. 

2.4   Shear-Layer  Characteristics 
The  equations  of  motion  for  the  two-dimensional,  incompressible,  turbulent 
mixing  layer  have  been  solved  by  Gortler  (in  White,  1974).  The  boundary  layer  equations 
are  solved  with  the  following  antisymmetric  boundary  conditions. 


u{-oo)=u,       Ui^co)=u,       uiO)=hu,.u,) 


(2.27) 
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The  solution  is 


u-\{u,^u„ 


1  + 


X 


(2.28) 


where  a»13.5. 

No  analytical  solution  has  yet  been  found  for  the  compressible  mixing  layer. 
However,  fiinctional  forms  of  the  mixing  layer  growth-rate  have  been  suggested  by  using 
order  of  magnitude  estimates.   Using  this  method.  Brown  and  Roshko  (1974)  suggested 


Ml  U  J 


(2.29) 


Bogdanoff  (1984)  suggested  that  the  convective  velocity  of  the  largest  structures 
in  a  mixing  layer  could  be  estimated  by  equating  the  stagnation  pressures  from  both  sides 
of  the  mixing  layer  with  respect  to  the  velocity  of  the  large  structures. 
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For  the  case  when  Yi^2.  the  above  equation  can  be  solved  explicitly  for  U^ 


^c- 


(2.31) 


and  the  convective  Mach  numbers,  M^j  and  M^j'  defined  as 
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u,-u  u-u, 

^i  '        ^2  (2.32) 

are  equal  and  can  be  written  as 


^i^^2  (2.33) 


It  can  be  seen  that  the  relative  Mach  number  (defined  in  Eq.  B.21)  Mj=2Mj  when  the 
average  speed  of  sound  is  used  to  calculate  M,. 

Papamoschou  (1986)  proposed  that  the  compressible  shear-layer  growth  rate  is 
functionally  related  to  the  incompressible  growth-rate  at  the  same  density  and  velocity 
ratios  such  that 

g/  _    Jh.(r,  s,  M^J 

So     -^•('■'  ^'  ^c,=0)  (2.34) 

This  isolates  the  effect  of  compressibility  so  that  the  above  ratio  is  a  universal  function 
for  free  shear-layers.  The  ratio  is  given  by  Benjamin  (1990)  using  a  curve  fit  to  the 
experimental  data  as 


^0  (2.35) 


Dimotakis  (1986)  proposed  that  there  are  three  main  phases  to  shear-layer 
entrainment,  which  are  referred  to  as  (1)  induction,  (2)  diastrophy,  and  (3)  infusion. 
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Induction  is  when  fluid  in  the  vicinity  of  the  vorticity-bearing  fluid  is  set  in  motion  by 
the  Biot-Savart-induced  velocity  field,  which  is  a  kinematic,  not  diffusive,  process. 
Irrotational  fluid  sufficiently  close  to  the  vortical  fluid  will  in  fact  participate  in  the  large- 
scale  structure  motions  long  before  it  has  acquired  vorticity  of  its  own.  These  motions 
appear  at  the  low  wavenumber  part  of  the  turbulent  spectrum,  and  although  irrotational, 
can  be  considered  part  of  the  turbulent  flow. 

The  second  stage  is  diastrophy,  where  the  irrotational  fluid  is  strained  until  its 
spatial  scale  is  small  enough  to  put  it  within  reach  of  the  (viscous)  diffusive  processes. 
Viscosity  then  takes  over  and  causes  cascading  to  the  Kolmogorov  scale. 

The  third  stage  can  be  associated  with  other  possible  diffusive  processes,  such  as 
molecular  mixing  or  heat  conduction,  and  may  or  may  not  precede  diastrophy,  depending 
on  the  relative  magnitude  of  the  corresponding  molecular  diffusivity  to  that  of  the 
kinematic  viscosity.  In  the  case  of  gas-phase  entrainment,  it  would  be  difficult  to 
distinguish  between  infusion  and  diastrophy  because  the  corresponding  diffusion 
coefficients  are  of  the  same  order. 

Dimotakis  (1986)  also  derived  the  incompressible  shear-layer  growth  rate  and 
accounted  for  unequal  entrainment  from  either  side  of  the  shear  layer 
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(2.36) 


where  e  is  a  constant. 

Nixon  et  al.  (1991)  calculated  the  transverse  mass  flux  around  a  vortex  core  using 
transonic  flow  theory  and  proposed  that  the  normalized  growth  rate  for  relative  Mach 
numbers  less  than  1 .4  is 
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8\r,s,M^)    _(l-^M,f 


8\r,s,M^=0)      |l-iM/)  (2-37) 

where  M^sAU/c.  Nixon  et  al.  also  found  that  a  considerable  amount  of  the  incoming  flow 
energy,  which  would  normally  be  converted  into  rotational  energy  of  the  eddies  within 
the  mixing  layer,  is  expended  compressing  the  gas. 

Papamoschou  (1989)  attempted  to  disrupt  the  structures  that  he  observed  in 
Schlieren  photographs  by  using  various  splitter  plate  trailing-edge  devices  (vortex 
generators,  trip  wires,  and  saw-tooth  extensions).  He  found  that  these  devices  had  little 
effect  on  turbulent  structures  or  shear-layer  growth  rate. 

Goebel  and  Button  (1990)  made  detailed  velocity  measurements  of  fully 
developed,  compressible  shear-layers  (0.40^  M^  ^1.97)  using  laser  doppler  velocimetry 
(LDV).  They  observed  that  the  length  for  full  development  generally  increased  in  the 
order  of  mean  streamwise  velocity,  streamwise  turbulence  intensity,  transverse  turbulence 
intensity,  and  kinematic  Reynolds  stress.  Fully  developed,  streamwise  mean-velocity 
profiles  were  well  approximated  by  an  error  function  profile,  even  for  the  more 
compressible  cases.  Transverse  turbulence  and  normalized  kinematic  Reynolds  stresses 
both  decreased  like  the  normalized  growth  rates  with  increasing  M^.  Mixing  lengths  and 
other  turbulence  quantities  with  a  transverse  fluctuation  component  also  decreased  with 
increasing  M^.  This  suggests  that  the  primary  effect  of  compressibility  on  the  turbulence 
field  in  a  shear  layer  is  the  suppression  of  transverse  velocity  fluctuations.  Streamwise 
turbulence  intensities  remained  fairly  constant  with  Mp  which  means  that  the  anisotropy 
of  the  turbulence  increased  significantly  with  M,  They  also  found  that  the  flowfields 
were  reasonably  two-dimensional  based  on  velocity  profiles  measured  at  two  parallel 
planes  in  the  spanwise  direction. 
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2.5   Jet  Acoustics 
There  has  been  a  great  deal  of  research  in  the  area  of  acoustics  due  to  the  major 
increase  in  commercial  air  travel  over  the  last  few  decades.    The  three  areas  relevant  to 
the  current  research  are 

(a)  acoustics  generated  by  turbulent  flow, 

(b)  acoustics  generated  at  nozzle  lips,  and 

(c)  acoustics  generated  from  imperfectly  expanded  nozzle  flow. 

A  subset  of  all  these  categories  is  the  interaction  of  any  one  source  with  a  solid  wall. 

The  clearest  way  to  understand  acoustics  generated  by  turbulent  flow  is  to  consider 
the  basic  equation  first  derived  by  Lighthill  (1952).  This  equation  is  the  basis  for 
virtually  all  work  on  noise  and  can  be  derived  from  the  mass  and  momentum  conservation 
equations: 


dt      dx.  (2.38) 


and 


dpu.    dpuu.    dp-- 
dt         dx.       dx. 


J  J 


(2.39) 


where  Q„  is  the  rate  of  mass  introduction  per  unit  volume  and  F;  is  the  external  force  per 
unit  volume.  The  nonhomogeneous  wave  equation  for  propagation  of  sound  in  a  uniform 
medium  at  rest  due  to  sources  of  matter,  external  forces  or  applied  fluctuating  stresses 


^P     ^id'pJQm      dF.d'T,. 

dt^     ^  dxf      ^t      d^i    ^^6  (2.40) 

where 


T.j=puu.^p,j-cip8.. 
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.  U,e  Ush*..,  equation,  c,  .s  *e  sound  speed  in  *e  n,edi„m,  and  5„  ,s  *e  Kronecker 
delta.  T^.  s^ess  tensor  for  a  Stokestan  gas  is  gtven  n,  .en.s  of  the  velocty  field  by 


P,rPh*l^ 


J«,_a»,^2Kl>  (2.42) 
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Effects  such  as  convection  of  sound  by  the  turbulent  flow,  or  .he  variations  in  the  speed 
of  sound  wtthtn  it.  are  taken  into  account  by  .ncon«rat,on  as  e,uivalent  stresses  (LighthtU 

1'53).  ^      ,.„       , 

.      ,     J      c^„   n  4(11  sives  the  effect  of  a  different 
Each  term  on  the  right-hand  side  of  Eq.  (2.40)  gives  me 

acoustic  source  mechanism.    The  fi.t.  aQ^a,  gives  the  effect  of  mass  mtroduction. 

Examples  include  pulse  jets,  tip  jet  rotors,  and  «te  random  mass  flucmations  that  can 

occur  across  the  exit  plane  of  a  jet  exhaust  (Pao  and  Lowson  1969).   The  second  ten.. 

aF/ax,  gives  the  effect  of  externa,  fluctuating  forces  that  can  act  on  the  fluid.  Examples 

include  compressors,  propellers,  helicopter  rotors,  and  the  fluctuating  forces  that  exist  on 

a  nozzle  lip  or  on  any  body  in  a  mrbulent  flow-stream.    The  third  teim,  aT,/ax,ax, 

incorporates  several  different  effects.   T„  is  called  the  "acoustic  stress  tensor."   At  low 

Mach  number,  the  most  significant  fluctuation  of  the  acoustic  stress  tensor  m  a  free  jet 

will  be  caused  by  the  turbulent  velocity  fluctuations,  which  affect  the  pu,u,  product.  The 

effect  of  viscous  stresses  and  thermodynamic  fluctuations  are  contained  in  the  remaining 

terms  of  Ty. 

curie  (1955)  investigated  the  influence  of  solid  boundaries  using  the  above  theory 

and  showed  that  the  fundamental  frequency  of  dipoles  near  the  wafl  is  one  half  of  that 
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generated  by  quadrupoles.  This  is  because  the  quadrupole  strength  per  unit  volume,  pUjUj, 
being  essentially  proportional  to  (velocity)^  will  have  double  the  frequency  of  the 
fluctuating  velocity.  On  the  other  hand,  the  fluctuating  force  exerted  on  the  fluid  at  the 
solid  boundaries  will  have  the  same  frequency  as  the  velocity  fluctuations.  The  frequency 
relation  can  be  shown  if  one  considers  a  flow  in  which  the  velocity  in  the  Xj-direction  is 
given  to  a  first  approximation  as 

u.  (f)«/4  .(f)cos(nf + q  (2-43) 

(where  ^j  is  also  a  function  of  x),  then  the  quadrupole  strength  per  unit  volume  is 


T.j^  puJUj^  pA.AjCOs(nt + C,)cos(m + Q 
=  o  P^/4f  os(2«f + C- + Cy) +cos(C,.  -  Cy)}  (2.44) 


Also,  at  sufficiently  low  Mach  numbers  the  contribution  to  the  sound  field  from  the 
dipoles  should  be  greater  than  from  the  quadrupoles.  Exactly  how  small  the  Mach 
number  must  be  before  this  occurs  will  depend  upon  the  flow  in  question. 

As  the  velocity  in  a  free  jet  is  increased,  the  maximum  angle  of  farfield  noise 
emission  starts  to  move  from  normal  to  the  jet  boundary  to  a  downstream-pointing  angle 
(Lighthill  1954).  For  jet  velocities  considerably  greater  than  the  surrounding  speed  of 
sound,  the  directional  maximum  becomes  the  Mach  angle  of  the  local  convective  Mach 
number.  Lighthill  also  noted  from  far-field  measurements,  that  low  frequencies  are 
emitted  in  a  more  downstream  direction  than  high  frequencies,  and  the  high  frequency 
sound  is  associated  with  the  mixing  region  just  downstream  of  the  no2:zle  where  the 
quadrupole  strength  is  due  to  lateral  transport  of  mean  momentum  across  the  shear  layer. 
The  low  frequencies  originate  from  the  ftilly  turbulent  region  downstream  (cf.  Fig.  2.3). 
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I„  superso„>c  jets  ^^^^^^^  ^^  ^„„,  „„ 

^^-^"rtrr::!, ... ... .  a„  a...... .o...  .e 


(2.45) 


.      ,       Tv.^  value  of  d  is  the  smaller 
^p  k  the  critical  value.  The  value  oi  u 

of  the  nozzle  dimensions.    Glass  U^  ^H  rate  and  decay  of  the  mixing 

_,,e*ac.ca„..U.w..evana..o„s.„.e.o.^-^^^^^^^^^^ 

,.......a...a.au.„.oo„e^^^^^^^^^ 

f„...bu,e„c=pass,ng.hroushashocUoff,n.tes.re„g*  ^^^^ 

.        ,-       The  values  range  from  0.0036%  at  M-l.Ol  to  62/. 
density  rat.a    The  values  r    B  „  Q82-|  studied  the  charactensUcs 

,       ,  M  7.  is  1  39%    Tam  and  Tanna  (1982)  stum 
number  lire  value  a.M.2  ^^  ^^^  ^_^^_^^  ^^  ^,„,, 

°'  — ^-"^'^"  TTtJ  K  :^^f  --   s  the  .et  Mach  number  and 
associated  noise  to  be  proportional  to  ^M,   M,  )  , 

M,  is  the  desi^  Mach  number.                                               ^^  ^^^^^^_^,^  .^ 

Tarn  (19.)  found  that  the  noise  of       ea^lv          V  ^^  ^_^^^^^  ^^_^^ 

en,anates..omt.orather.oca,i.dre.onsof.      ewh                   ^  ^^^^^.^^^^ 

.^downstream  oftheno.lee.it.  Hesusses^^atar^^^^^^^  ^  ^^^^^^^^^ 
■u,    fnr  transfemng  the  kinetic  energy  of  the  jet  inio 
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A  .„  interact  strongly  with  the  arttbient  Huid  there, 
the  mtxmg  layer  at  two  locations  and  to  tnteract 

hence  g.v.ng  rise  to  intense  notse  radtatton.  ^^^  ^^^^^ 

radiation  from  the  nozzle,  ^.u;  ,,uprical  waves  from 

^  Wpcjtlev  ri973)  also  observed  spherical  wdvc^ 

^■■•^"^^"^^^""^    ::rC   — dto..rhnlencew.th,n,d.a.eter 
the  nozzle  of  a  supersonic  hekum  jet  Th  J  „  „.  g;  kHz,  and  directed 

.„a  radtated  wea.  shoc^  fron,  this  locat.on  w,th  a  frequency  of       ^ 
primarily  along  a  cone  60"  from  the  jetaxts. 
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Figure  2.1 
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Figure  2.3        Flow  regimes  in  the  ejector. 


CHAPTERS  PROCEDURES 

1  1    GeneralDescriEtion 

3„.c„.  seconda^  w,nd  W  «*        -P  ^^^^^^^^  „,  ,„,,..  by 

pn^a.  an.  seco„.a..  .  -  «Pe— J-      ^^  ^^^  ^^^  ^^,^  ^^,,  .. 
.He  C3as  Dynamics  Ubo..oo,  ^■""7;^;;  ^^,,,,  („,  ,es,gn)  Mach  nu»be.,  the  *ee 

,.0„aa..pp>y.  ----^""^";  t      j.en..ca,  spU«e.  ,a.es,  an.  n,.e. 
_.  .e  .0.^.  -  eon.c.  -^'^^ ^     ,,,  „a.  e.  of  *e 

,p.a.s  .s  connecea  fo  a  —  '^^  ;_,,^„  .„  .  ,  ,  ....  cc-seCon 
.0.  a  SOS  n,.  .y  -.  n,.  ..an.>-o^^  _^  ^^^__  ^^^^  _^^^  ^  ,^^^^  „„,..  .. 
«teh  houses  a  butterfly  control  valve.  ^^^^^^^^^^  p^„,p,p,, 

— -°*-"°^^^^'^  ^":;;:eCip"-'---"-^*"^^ 

«.euse.tov,sua,.ethesUuctu.onh^^^  ^^^_^^^^^^^^^     ^^^^^ 

,,.,  .a  the  extent  of  tbe  — t.^^,  _  ^en  a.  ..enty-flve  st.ea.«.se 
,onies  ob.a,nea  us.n,  a  t.ve.,ns  --  ^^^^^_^^^^  ^^^^^^^  ^^^  ^^_,,  ,..e 
p„s,..o„s  «.tl,  nominal  separafons  of  20  m  ^^^^  ^^^^__  ^^  ^^^  ^^„,  ^, 

.     ,      .„,  a  30-tube  manometer  bank.    Photogr  P  ,,  „ere  obtained 
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AH  of  the  measurements  are  accompanied  by  a 

Schheren  and  shadowgraph  ^^  ^^^^^^^  ^^^^^^^ 

„easu.eme„.  «ere  aU  controlled  us.ng  a  PC-AT  co^P 

.■      „f  Aereqmred  functions  during  the  short  durat,on  of  the  runs.  Th     pe 
execution  of  fl,=  requ.re  ^_^  ^^^^_^^  ^^^_^  ^  ^^  ^  ^^„,, 

of  the  blow-down  fachty  was  mterm.ttent,  w,th 

.„,„t  and  visualization  techniques,  as  wen  a^ 
The  ejector  hardware,  pressure  measurement  and 

boundarV-layer  calculations  are  descnhed  in  detatl  below. 

J  o   r>p^ipn  Methodology, 
.He  two-dintensronal  ejector  r,gusedin.h,sstud,wasapre-ex,sti„sfac,h.y«t 

.     ,.,  ,931  in  1 992  for  dissimilar  gas,  ejector 
by  the  auftor  and  his  coworkers  (Benjamm  et  al.  1993) 

nnents     The  r.g  allows  for  vanation  m  m^ng  section  he.ght  from  38.1  mm  to 
expenments.    The  ng  „„f,2  7mm    In  this  stirdy,  tite  maxtmum 

5J.4  mm  wtth  a  f«ed-he,^t  pnmary  nozzle  of  .2.7  mm^ 
,e.ghtof51..mm  was  Chosen  sothatthemaxtmumnumberofwallmeasurem 

be  made  in  that  direction  (see  Hg.  3.1).  ^^ 

The  mmtmum  allowable  dtmens.on  of  the  nozzle  ex.t  he.^t  wa^  s. 

.  H  hv  the  tiaversing  pitot-pressure  probe,  and  (b)  the  number 
„„:  (a)  tite  resolution  provided  by  the  g  P      P  ^  _^  ^^ 

of  .averses  necessary  to  map  the  vartous  regions  of        sh         y 

■     ,  H  that  the  resolution  of  the  pitot-pressure  probe  was  no  less  titan  0.2  mm, 

ments  should  fall  across  the  layer  (m  fact,  many  more  readmgs  were 
^"^-""■""  ,  '  .abovemeantthatthelayermustgrowtoathickness 
obtained  across  the  shear  layer).  The  above  m 

u  f      *,  thickness  can  be  measured  with  the  probe, 
of  about  2  mm  before  the  tiiickness  _,„y.,p.ced  streamwise 

Secondly,  it  was  decided  that  a  mmimum  of  tiiree  equa  y   P 
Tere  n  cessary  to  identify  the  pnma^  core  region  of  the  shear  layer  (region  1). 
traverses  were  necesimy  iv  ^ 
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seconda.,  and  (.)  .he  ,aye.  ,ows  w,.  *=  sa^e  ra.  ,„  bo.h  d.ec.,o„s  ,e,a.,ve  .0  the 
„„,„,.sec.,o„  ce„.e„^e).  Th.s  ™ea„.  *a.  the  s.aHest  ™x,„g  sect.on  that  ccu.d  be  use 
„„uM  result  ,„  a  layer  that  wouM  he  about  6  nttn  .hi*  at  the  end  of  .he  core  (3  nt.  each 
..e  of  the  centerhne).  That  ,n  utm  nteant  that  the  nozzle  exU  height  had  to  be  at  least 
6  ntm.  A  factor  of  two  was  mtroduced  to  account  for  any  n,arg,n  of  uncertatnty;  hence 

a  primary  exit  height  of  12.7  mm  was  chosen. 

The  exit  height  of  the  seconda^  s.rean,s  was  selected  based  on  s.mi.ar  arguments. 
^dd,t,o„a,ly,  ,t  was  desrred  that,  at  least  m  some  test  cases,  umnbced  secondary  flow  extst 
beyond  the  potnt  where  the  pnmary  core  ends.  Usmg  these  cntena,  .t  was  decde   th. 
secondary  flows  having  ex.t  heights  a.  leas.  e,„a,  .o  *e  pnmary  jet  height  would  be 
adequate.  Hence  a  variable  ex.t  height  of  .2.70  mm  to  19.05  mm  was  chosen  for  each 

half  of  the  secondary  stream. 

in  selectmg  the  depth  of  the  m^ng  secfon  (see  F.g.  3.1)  tt  was  destred  to 
^mimtzeflowratesaswell  as  provrdeaflowfieldthatwas  approximately  twcdimensrona, 

(,  e  stdewall  boundary-layers  are  thin  compared  to  the  w.dth),  and  be  adequately  s.zed 
for  Schlieren  vrsualization.  The  latter  requrred  that  the  mrxing  section  was  deep  to 
resolve  small  densi.y-grad.ents.  Smce  these  are  opposing  factors,  an  acceptable  mtxmg- 
.ection  wrdth  that  reasonably  met  the  cnterta  was  detennined  to  be  25.4  mm.  Complete 
ae.,ls  of  the  ng  conceptton  and  des.gn  can  be  found  m  the  report  by  Roan  et  al.  (1993). 

I  ^   Pvperimental  Apparatus 
in  the  m^ing  section,  the  supersonic  primary  center  jet  is  sandwiched  between  two 
symmetrica,,  subsomc  secondary  streams  and  all  three  d.scharge  into  a  constant-area  duct 
(see  Fig  3  1).   Pnmary  and  secondary  streams  mrx  along  the  shear  layers  which  beg,„ 
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development  of  the  flowfield  is  mapped  by  traversing  a  pitot-pressure  impact  probe  across 
the  mixing  duct  at  a  number  of  axial  locations. 

The  primary  and  secondary  flows  were  supplied  by  a  high-pressure  gas  storage  and 
delivery  system.  Control  of  the  rig,  as  well  as  data-acquisition  of  the  mean  conditions, 
was  accomplished  using  the  laboratory  PC-AT  computer.  This  control  process  enabled 
short  run  times  and  identically-repeatable  instrument  reading  capability.  The  test 
apparatus  as  well  as  the  interface  hardware  and  instrumentation  are  described  in  more 
detail  in  the  following  sections. 

Fine  adjustments  of  primary  and  secondary  supply  pressures,  as  well  as  the  mixing 
section  back-pressure  had  to  be  made  for  each  test  case  to  achieve  the  desired  flow 
conditions.  To  achieve  repeatable  pressure  conditions  between  test  runs,  the  run  event 
sequence  and  timing  was  carefully  planned  and  implemented  to  ensure  stability  and 
validity  of  the  measurements  while  minimizing  run  time.  The  details  of  the  set-up,  run- 
control,  and  data-acquisition  procedures  are  also  discussed  later  in  this  chapter. 

3.3.1    Gas  Storage  and  Delivery  Svstem 

Pressurized  air  was  used  for  supply  of  the  primary  and  secondary  flows.  The 
filtered  and  dried  air  was  delivered  from  a  bank  of  storage  tanks,  and  supplied  to  the 
primary  and  secondary  plenums  through  pressure-regulating  valves  as  desired.  A 
schematic  of  the  gas  storage  and  delivery  system  is  shown  in  Fig.  3.2. 

Compressed  air  was  supplied  by  a  Worthington  two-stage  positive-displacement 
340  cfm  (150  HP)  compressor  which  delivers  a  maximum  working  pressure  of  300  psia. 
The  compressed  air  is  passed  through  several  oil  and  water  traps  and  a  desiccant  air-drier 
to  remove  any  moisture.  The  air  is  then  stored  in  a  bank  of  cylindrical  tanks  manifolded 
together,  with  a  total  volume  of  approximately  420  ft'.  In  order  to  maintain  an  almost 
constant  supply  pressure  level,  the  compressor  was  operated  throughout  test  runs.    In 
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addition,  the  pressure-regulatmg  valves  were  set  when  the  supply  pressure  reached  a 
predetermined  value;  this  value  is  determined  after  the  compressor's  second  stage  is 
disengaged  at  its  maximum  deliveiy  pressure.  Monitoring  of  the  supply  pressure  is 
performed  using  a  transducer  connected  to  the  supply  line. 

Air  was  delivered  to  the  apparatus  through  a  2  m.  schedule-80  steel  pipe  sized  to 
keep  the  pipe  Mach  number  below  0.4  for  the  flow  rates  of  interest.  A  gate  valve  placed 
between  the  tanks  and  the  test  rig  was  used  to  isolate  the  supply  and  the  test  rig. 
Automatic  on/off  control  of  the  air  flow  was  accomplished  using  a  solenoid  valve 
(Omega,  model  SV207)  placed  downstream  of  the  gate  valve.  The  solenoid  valve  could 
be  controlled  either  manually  with  a  switch,  or  electronically  using  a  digital  output 
chamiel  from  the  computer.    The  air  could  then  be  delivered  to  the  primary  and/or 
secondaiy  supply  Unes.  On  the  secondary  branch,  control  of  the  secondaiy  flowrate  was 
effected  utilizing  3/4  in.  and  1/4  in.  needle  valves  connected  in  parallel.  On  the  primary 
branch,  a  3/4  m.  needle  valve  and  a  3/4  m.  globe  valve  connected  in  parallel  controlled 
the  flow  rate  of  primary  air.   On  each  branch,  the  large  and  the  small  valves  were  used 
for  gross  and  fine  flow-rate  adjustment  as  required. 

^  •^  7   Primarv  and  Slecnndarv  Lines 

The  gas  flows  reqmred  for  the  primary  and  secondary  were  delivered  to  the 
respective  plenums  from  the  pnmary  and  secondary  branches  descnbed  m  the  previous 
section.  The  pnmary  line  was  constnicted  from  3/4  m.  schedule-80  steel  pipe  and 
attached  to  the  primary  plenum  by  a  flexible  metal  hose  with  an  inside  diameter  of 
19.05  mm.  This  attachment  method  isolated  the  ejector  from  any  mechanical  vibration 
due  to  the  flow  through  the  supply  piping. 

The  secondary  flow  was  supphed  to  the  secondaiy  plenums  through  four  steel 
reinforced  PVC  hoses.    These  hoses  were  connected  to  a  manifold  consisting  of  four 
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sections  of  1  in.  pipe  welded  at  90°  from  each  other  around  the  capped  end  of  the  2  in. 
secondary  supply  line.  A  schematic  of  the  air  supply  hardware  described  above  is  shown 
in  Fig.  3.2. 

3.3.3   Plenums 

The  primary  and  secondary  plenums  consist  of  flow-conditioning  sections  followed 
by  settling  chambers  (Figs.  3.1  and  3.3).  The  walls  of  these  sections  were  machined  from 
aluminum  stock  and  formed  continuous  flow  passages  of  rectangular  cross  sections.  The 
primary  plenum  was  in  the  center,  sandwiched  between  two  identical  secondary  plenums. 
Two  large  rectangular  plates  12.7  mm  thick  served  as  side  walls  for  the  primary  and 
secondary  plenums.  The  remaining  walls  were  machined  from  25.4  mm  thick  aluminum 
plate  which  were  (a)  the  internal  walls  separating  the  primary  and  secondary  plenum 
cavities,  (b)  the  rear  wall  for  all  plenums,  and  (c)  the  outer  walls  of  the  secondary 
plenums.  Two  stainless  steel  dowel  pins  inserted  through  each  wall  in  (a)  and  (c)  and 
into  location  holes  drilled  in  the  two  side  walls,  accurately  position  each  part  of  the 
plenum  assembly.  All  plenum  walls  are  bolted  together. 
3.3.3.1    Primary  plenum 

The  primary  plenum  consists  of  a  flow  conditioning  section  and  a  settling 
chamber.  It  has  a  rectangular  cross-section  59.7  mm  high  by  25.4  mm  wide.  These 
dimensions  approximately  constrain  the  plenum  flow  velocity  between  3  m/s  and  27  m/s 
for  a  wide  range  of  operating  conditions,  as  suggested  by  Pope  and  Goin  (1978). 

The  flexible  metal  hose  which  delivers  the  primary  gas  to  the  plenum  fastens  onto 
a  thick- walled  aluminum  nipple  with  a  25.4  mm  inside  diameter.  This  inlet  nipple  is 
welded  at  right  angles  to  one  of  the  plenum  side  walls.  The  gas  then  passes  through  two 
different-sized  perforated  plates  machined  from  aluminum,   a  section  of  aluminum 
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honeycomb,  and  two  different-sized  stainless  steel  w,re  meshes  (ef.  F.gs^  3.1  and  33). 
These  flow-ntanagement  deviees  help  ensure  un.fonn  flow  and  low  turbulence  intens.ty. 
The  perforated  plates  are  3.9  tnm  th.ck  and  spaced  12.7  mm  apart  (.6.6  mm  from 
center  to  center).    The  plates  fit  into  slots  machined  ,„  all  four  walls  of  the  plenum 
passage.  The  first  plate  is  27.4  mm  from  the  m.et  mpple  centerHne.  and  has  a  12x5  array 
of  4  42  mm  diameter  holes;  the  center-to-center  spacing  of  the  holes  is  4.98  mm  and 
5  08  mm  along  the  long  and  short  d.mensions  of  the  plenum,  cross  sect.on,  respecfvely. 
The  second  plate  has  a  .4x6  array  of  3.68  mm  d.ameter  holes;  the  spacmg  in  th,s  array 
is  4  27  mm  and  4.23  mm  along  the  long  and  short  dimens.ons,  respectively. 

The  honeycomb  sectron  is  located  34.9  mm  downstream  from  the  second  plate. 
This  section  ,s  82.6  mm  long,  and  has  a  3.18  mm  hexagonal  cell  s.ze,  and  fits  ,n  a  recess 
machmed  in  all  four  walls  of  the  plenum  (cf  Figs.  3.1  and  3.3). 

The  first  of  the  two  mesh  screens  rs  approximately  25.4  mm  downstream  of  the 
honeycomb,  and  has  a  mesh  of  30x30  wires/in.  and  a  wtre  d.ameter  of  0.165  mm 
resulting  in  a  35.2%  obstructron  and  a  wire-spacmg  to  wire-d,ameter  ratio  of  5.13.  The 
second  mesh  is  .2.7  mm  downstream  from  the  firs,  one,  and  has  a  mesh  of  40x40 
wires/in.  and  a  0.165  mm  wire  diameter,  resulting  in  a  45.2%  obstruction  and  a  wire- 
spacmg  to  wire-diameter  ratio  of  3.85.  THe  screens  are  soldered  to  brass  frames  which 
fit  in  slots  machined  in  the  plenum  walls.  Both  dte  screens  and  frames  are  shown  m  Figs. 


3.1  and  3.3. 


I  -J.J. 

The  modular  design  of  the  flow  managements  devices  enabled  easy  inspection, 

cleaning  or  replacemem  of  any  component. 

The  prima.,  gas  enters  the  settling  chamber,  which  extends  44.8  mm  downstream 
of  the  second  screen.  A  1.59  mm  diameter  T-type  thermocouple  is  located  between  the 
first  and  second  screens  to  measure  total  pnmary  temperature,  and  a  0.79  mm  d.ameter 
pressure  tap  located  on  one  side  wall  25.4  mm  downstream  from  the  second  screen  is 


43 
used  to  measure  primary-plenum  stagnation  pressure.'  The  outside  of  the  wall  is  tapped 
for  mounting  a  pressure  transducer. 
3.3.3.2   Secondary  plenums  and  flow  passages 

Each  secondary  plenum  consists  of  a  flow  conditioning  section  and  a  settling 
chamber  similar  to  those  of  the  primary  plenum.  Their  rectangular  flow  areas  are 
71.9  mm  high  and  25.4  mm  wide,  which  results  in  a  local  velocity  of  26  m/s  for  the  flow 
conditions  used  in  this  investigation  (see  Chapter  4  for  flow  conditions). 

The  secondary  gas  is  delivered  to  each  plenum  through  two  flexible  hoses  clamped 
to  inlet  nipples  that  are  welded  at  90°  on  opposite  sides  of  the  plenums.  The  nipples  are 
made  of  aluminum  and  have  inside  diameters  of  25.4  mm.  In  each  plenum,  the  gas 
passes  through  flow  management  devices  nearly  identical  to  those  used  in  the  primary  (cf. 
Figs.  3.1  and  3.3).  Two  different  perforated  plates  machined  from  aluminum,  a  section 
of  aluminum  honeycomb,  and  two  different  stainless  steel  wire  meshes  ensure  uniform 
flow  and  low  turbulence  intensity  of  the  secondary  stream. 

The  perforated  plates  are  3.94  mm  thick  and  are  12.7  mm  apart  (16.7  mm  from 
center  to  center).  The  plates  fit  in  slots  machined  in  all  four  walls  of  the  plenum  passage. 
The  first  plate  is  27.38  mm  from  the  inlet  nipples  centerline,  and  has  a  14x5  array  of 
4.45  mm  diameter  holes;  the  center-to-center  spacing  of  the  holes  is  5.18  mm  and 
5.13  mm  along  the  long  and  short  dimensions  of  the  plenum's  cross-section,  respectively. 
The  second  plate  has  a  17x6  array  of  3.66  mm  diameter  holes;  the  spacing  in  this  array 
is  4.24  mm  and  4.23  mm  along  the  long  and  short  dimensions,  respectively. 


'  Although  this  method  actually  measures  the  static  pressure  and  static  temperature,  the 
plenum  velocity  is  approximately  29  m/s  resulting  in  deviations  of  0.5%  from  the  true 
stagnation  pressure,  0.15%  from  the  true  stagnation  temperature,  and  0.2%  below  the 
true  Mach  number. 
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The  honeycomb  section  and  screens  are  the  same  as  those  described  above  for  the 
primary  plenum  except  for  the  overall  height,  and  are  located  at  the  same  streamwise 
positions. 

The  secondary  gas  enters  the  settling  chamber  which  is  97.5  mm  long.  A  pressure 
tap  located  25.4  mm  downstream  of  the  second  screen  allowed  for  measurement  of  total 
pressure  in  each  secondary  settling  chamber.  This  consisted  of  a  0.79  mm  diameter  hole 
drilled  through  the  side  wall  of  each  secondary  plenum,  and  tapped  on  the  outside  for 
mounting  pressure  transducers. 

Following  the  settling  chamber  is  a  converging  channel  of  rectangular  cross- 
section  which  accelerates  the  flow  into  the  mixing  section  (cf.  Figs.  3.1  and  3.3).  It  has 
a  constant  depth  of  25.4  mm  and  the  height  varies  monotonically  from  71.9  mm  to 
18.9  mm.  The  passage  is  formed  by  the  plenums  side  walls,  an  "inner"  contoured  wall, 
and  an  "outer"  contoured  wall.  The  "inner"  contoured  wall  corresponds  to  the  outer 
surface  of  the  nozzle  block.  The  machining  of  the  "outer"  wall  provided  an  accuracy  of 
approximately  ±0.05  mm.  Both  contoured  walls  have  a  shape  defined  by  fifth-order 
polynomials  with  zero- valued  first  and  second  derivatives  at  the  beginning  and  end  of  the 
contraction.  The  length-to-height  ratio  of  the  contraction  is  2.  Each  secondary  stream 
discharged  into  the  mixing  section  at  the  end  of  the  contraction. 

3.3.4   Primary  Nozzle 

The  primary  gas  exits  the  settling  chamber  and  discharges  into  the  mixing  section 
through  a  two-dimensional  supersonic  nozzle.  The  nozzle  was  machined  from  aluminum 
to  provide  the  required  flow  Mach  number  of  1.8  for  gases  with  a  constant-pressure 
specific  heat  ratio  of  1 .4. 

The  nozzle  consists  of  two  identical  blocks  bolted  to  the  internal  walls  of  the 
plenum  and  clamped  between  the  side  walls  (Figs.  3.4).    The  blocks  are  mounted  on 
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opposite  sides  of  the  plenum  centerline,  and  provide  the  desired  change  in  flow  area. 
Two  stainless  steel  dowel  pins  inserted  through  each  block  and  into  location  holes  drilled 
in  the  plenum  side  walls  position  the  blocks  accurately.  Each  block  has  two  contoured 
surfaces;  the  "inner"  surface  corresponds  to  the  contoured  walls  of  the  primary  nozzle, 
while  the  "outer"  surface  serves  as  the  wall  for  the  secondary  passages.  One  of  the 
plenum's  side  plates  is  split  into  two  sections  to  facilitate  easy  exchange  of  different 
nozzle  blocks. 

The  initial  portion  of  the  nozzle  blocks  is  straight  and  serves  as  an  extension  to 
the  walls  of  the  constant-area  settling  chamber.  After  this  section  the  blocks  provide  a 
contraction  region  where  the  flow  begins  to  accelerate.  The  shape  of  this  region  is 
defined  by  a  fifth-order  polynomial  with  zero-valued  first  and  second  derivatives  at  the 
beginning  and  end  of  the  contraction.  The  contraction  region  has  a  length-to-height  ratio 
of  2.5  and  ends  at  the  nozzle  throat.  This  profile  design  provides  uniform  flow  and 
ensures  that  flow  separation  does  not  occur  (Papamoschou  1986). 

Following  the  contraction  is  the  diverging  region  of  the  nozzle,  the  shape  of  which 
was  obtained  using  the  two-dimensional  method  of  characteristics  technique,  and 
neglecting  boundary  layer  corrections.  The  intersection  between  the  converging  and 
diverging  regions  provide  a  sharp  edge  at  the  nozzle  throat. 

The  nozzle  blocks  have  a  trailing-edge  thickness  of  approximately  0.5  mm.  Once 
installed,  the  nozzle  extends  177.8  mm  beyond  the  second  screen  at  the  entrance  to  the 
settling  chamber.  The  contoured  walls  of  the  nozzle  blocks  were  machined  to  within 
±0.05  mm.  The  installed  throat  height,  2/j„  is  8.79±0.03  mm  compared  to  the  design 
height  of  8.84  mm,  and  the  installed  nozzle  exit  height,  2h,  is  12.65±0.03  mm  compared 
to  the  design  height  of  12.70  mm.  The  installed  dimensions  result  in  a  calculated  value 
of  Mj=1.80  based  on  1-D  isentropic  gas-dynamics  area  relations.  The  length  of  the 
diverging  section  of  the  nozzle  is  16.2  mm. 
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3.3.5   Mixing  Section 

The  primary  and  secondary  streams  come  into  contact  at  the  exit  plane  of  the 
plenums  (cf.  Figs.  3.1  and  3.3).  At  this  point  the  streams  begin  to  mix  in  a  constant-area 
section  of  rectangular  cross-section  (Figs.  3.1  and  3.5).  The  section  is  25.4  mm  wide  and 
51.4  mm  high.  It  is  formed  by  two  side-walls,  a  top  wall  and  a  bottom  wall  which  boh 
together  providing  the  desired  flow  passage.  The  section  is  bolted  at  one  end  to  the 
plenum,  and  it  extends  500  mm  beyond  the  exit  plane  of  the  plenums.  For  pitot  traverse 
and  wall  static-pressure  measurements,  both  side-walls  were  machined  from  aluminum. 
For  turbulent  wall-pressure  measurements,  the  measurement  side-wall  was  machined  from 
brass.   The  top  and  bottom  walls  were  machined  from  brass. 

The  side  wall  used  for  mean,  static  wall-pressure  measurements  has  29  holes 
drilled  through,  all  of  which  are  1.59  mm  diameter.  Twenty-seven  of  these  holes  are 
spaced  in  the  stream  wise  direction  every  10  mm  along  the  centerline  of  the  wall,  with  the 
first  one  positioned  so  that  it  lies  at  the  exit  plane  of  the  plenums  on  the  primary 
centerline.  This  hole  is  used  to  measure  the  primary-exit  static  pressure.  The  other  two 
holes  also  lie  at  the  exit  plane  of  the  plenums  and  are  located  13.2  mm  on  either  side  of 
the  primary  centerline,  in  the  secondary  flow.  These  holes  are  used  to  measure  the 
secondary  streams'  exit  static  pressure.  A  short  section  of  stainless  steel  surgical  tubing, 
inserted  and  glued  in  a  counterbore  provided  with  each  hole,  is  used  to  connect  pressure 
transducers  via  1.59  mm  internal  diameter  flexible  plastic  tubing.  The  plate  with  wall 
static  pressure  taps  is  shown  in  Fig.  3.6. 

For  the  time-accurate  wall-pressure  measurements,  an  aluminum  frame  holding 
interchangeable  brass  plates  was  designed  with  the  same  dimensions  as  the  side  wall  used 
for  mean,  static  wall-pressure  measurements.  Two  types  of  interchangeable  brass  plates 
were  constructed:  solid  "blank"  plates,  and  plates  with  hole  patterns  for  static-pressure 
transducers.    The  plates  can  be  placed  at  various  locations  along  the  frame.    The  holes 
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in  the  plates  with  transducer  hole-patterns  are  reamed  to  2.64  mm  diameter.  Brass  plugs, 
or  sleeves  that  hold  pressure  transducers  of  2.615  mm  diameter,  fit  in  the  holes  and  are 
shown  in  Fig  3.7.  An  0-ring  seals  the  plugs  0.89  mm  from  the  flow  surface,  while  the 
sleeve  does  not  have  a  seal  due  to  constraints  of  size.  This  was  expected  not  to  be 
significant  because  there  is  a  metal  to  metal  seal  where  the  sleeve  head  is  fastened  to  the 
back  of  the  plate.  Acoustic  disturbances  caused  by  the  gap  between  the  sleeve  and  the 
hole,  such  as  a  Helmholtz  resonator,  have  a  calculated  frequency  on  the  order  of  150  kHz 
that  is  well  above  what  is  measured  here.  A  typical  installation  is  shown  cutaway  in  Fig. 
3.1,  while  Fig.  3.8  shows  the  dimensions  of  the  assembled  frame  and  plates  used  in  this 
research.  The  transducer-hole  spacing  is  either  4.09  mm  or  8.17  mm  in  the  streamwise 
direction,  and  3.44  mm  in  transverse  direction. 

The  top  and  bottom  walls  each  have  14  ports  machined  along  their  centers  to 
provide  access  for  the  stem  of  the  traversing  impact  probe.  The  streamwise  spacing 
between  ports  is  20  mm  and  the  first  port  is  25.4  mm  downstream  from  the  mixing 
section  inlet.  Custom-made  threaded  brass  plugs,  each  fitted  with  an  0-ring,  provide  a 
leak-tight  seal  for  the  ports.  The  nonalignment  of  the  plug  ends  was  minimized  by 
material  removal  using  in  turn,  220,  320,  400  and  600  grain  size  abrasive  paper  with  the 
plugs  fitted  in  their  respective  ports.  Final  buffing  of  the  surfaces  was  accomplished 
using  a  mixture  of  jeweler's  rouge  and  tallow  that  produced  better  than  a  0.2  \im  surface- 
roughness  finish.  The  resulting  nonalignment  of  the  plug  ends  can  thus  be  confidently 
stated  as  being  less  than  0.015  mm  from  the  surface  of  the  wall.  Two  additional  plugs 
were  used  in  connection  with  the  traversing  probe.  These  plugs  have  a  2.29  mm  diameter 
hole  drilled  through  their  center  in  which  the  stem  of  the  probe  slides.  Figure  3.9  shows 
the  top  wall  and  Fig.  3.10  shows  the  brass  plugs. 
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-<  ^  6   Travet^tnp  Tm pact-Pressure  Probe 

The  traversing  total-pressure  impact  probe  constructed  for  this  experimental  study 
is  shown  in  Fig.  3.5.  As  a  result  of  the  vanable  probe-interference  effect  noticed  durmg 
previous  tests  (Roan  et.  al,  1992),  a  traversing  probe  arrangement  that  results  in  a  nearly 
constant  obstmction  of  the  flow  was  constructed.  TTiis  was  accomplished  by  attaching 
an  extension  to  the  probe  stem  (i.e.,  a  dummy  stem  beyond  the  sting).  The  stem  and  the 
extension  slide  through  the  holes  dnlled  in  the  plugs  that  are  used  with  the  probe. 

To  simphfy  the  procedure  necessaiy  to  change  the  probe  position,  a  probe  was 
designed  that  allows  for  removal  of  the  stmg.  This  was  accomplished  with  a  special 
coupling  machined  from  stainless  steel.  The  sting  of  the  probe  was  soldered  to  the  male 
part  of  the  coupling.  The  male  part  was  screwed  into  the  female  part  to  which  the  stem 
and  the  dummy  stem  were  soldered.  A  Teflon  bushing  sealed  the  clearance  between  the 

parts. 

A  probe  sting  was  constructed  from  hard  tempered  stainless  steel,  hypodermic 

tubing  having  outside  and  inside  diameters  of  0.90  mm  and  0.58  mm  respectively.   The 

sting  IS  25.4  mm  long,  and  the  nose  of  the  sting  has  a  screwdriver  wedge-like  shape  with 

a  rectangular  opening  of  approximately  0.89  mm  by  0.13  mm.  To  obtain  this  geometiy, 

a  piece  of  0.127  mm  thick  shim  stock  was  inserted  approximately  4  mm  into  the  end  of 

the  tube,  and  then  clamped  together.    Care  was  taken  to  insure  that  the  tube  collapsed 

evenly  about  a  plane  passing  through  the  axis  of  the  tube.    Hie  flattened  end  was  then 

shaped  to  a  wedge  using  a  sharpening  stone,  and  then  polished  to  a  mirror  surface  with 

600  grain  size  abrasive  paper  and  348  rouge. 

The  stem  of  the  probe  and  its  extension  are  each  roughly  100  mm  long,  and  were 

made  from  hard  tempered  stainless  steel  tubing  having  outside  and  inside  diameters  of 

3.05  mm  and  2.39  mm,  respectively.   The  open  end  of  the  stem  was  comiected  directly 

to  the  pressure  transducer  with  a  compression  fitting. 
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The  probe  is  traversed  across  the  mixing  section  by  means  of  a  computer- 
controlled,  stepper-motor-driven,  leadscrew.  The  leadscrew  assembly  is  mounted  on  top 
of  the  upper  brass  wall  of  the  mixing  section.  The  assembly  slides  in  a  rail  built  into  the 
brass  wall,  and  can  be  locked  at  any  desired  streamwise  position.  The  probe/leadscrew 
assembly  is  shown  in  Fig.  3.5.  The  leadscrew  has  a  pitch  of  10.16  mm.  The  stepper 
motor  provides  200  steps  per  revolution,  while  the  motor  driver  allows  for  half-step 
motion.  This  arrangement  translates  into  a  linear  resolution  of  0.0254  mm  for  the  traverse 
of  the  probe.  The  time  required  for  the  46.99  mm  traverse  was  2.00  seconds.  The 
pressure  transducer  used  in  conjunction  with  the  probe  (see  section  3.4  for 
instrumentation)  was  mounted  on  a  bracket  bolted  to  the  traversing  nut  of  the  leadscrew 
assembly. 

3.3.7   Transition  Diffuser  and  Exhaust 

Immediately  downstream  of  the  mixing  section  is  a  760  mm  long  transition- 
diffuser  with  a  rectangular  cross-section  inlet  of  50.8  mm  by  25.4  mm  and  a  circular 
cross-section  outlet  diameter  of  102.3  mm.  The  design  has  a  4.82°  average  included 
divergence  angle,  and  an  area  ratio  of  6.4  (see  Fig.  3.11). 

The  diffuser  is  constructed  from  two  identical  sections  welded  together,  each  made 
from  4.76  mm  thick  304  stainless  steel  sheets  bent  to  the  desired  shape  using  about  20 
longitudinal  bends.  Flanges  were  then  welded  to  each  end  suitable  for  mating  upstream 
to  the  mixing  section  and  downstream  to  the  4  in.  butterfly  valve  (Keystone  992). 

The  flow  is  exhausted  to  the  atmosphere  by  way  of  4  in.  diameter  PVC  schedule 
40  pipe,  10  in.  diameter  schedule  40  steel  pipe,  and  through  a  Maxim  BR31  silencer  to 
the  atmosphere.  These  components  are  shown  schematically  in  Fig.  3.11. 
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3 .4  Instrumentation 
A  schematic  of  the  instrumentation  and  computer  interface  is  shown  in  Fig.  3.12. 
Pressure  and  temperature  data  were  measured  using  various  pressure  transducers  and  T- 
type  thermocouples.  These  devices  were  connected  to  two  16-channel  multiplexers 
(Acces  AIM-16(P))  using  shielded  twisted-pair  wire  for  the  pressure  transducers  and  steel 
braid  shielding  for  the  thermocouple  wires.  The  multiplexers  were  connected  to  a 
moderate-speed  analog  and  digital  (A/D)  computer  board  (Acces  AD  12-8),  with  12-bit 
data  resolution,  mounted  in  a  12  MHz  AT-type  personal  computer  (PC).  Using  Quick 
Basic  version  4.5  software,  the  maximum  A/D  throughput  was  approximately  9,000  single 
channel  samples/second. 

Digital  output  from  the  A/D  was  used  to  control  two  electronic  relays  (Crydom 
S440)  supplying  llOVAC  to  the  flow  solenoid  shutoff  valve  (OMEGA  SV207),  and  a 
solenoid  valve  used  for  venting  the  air  tank  to  the  atmosphere  (OMEGA  SV207). 

Three  models  of  pressure  transducers  were  used  at  the  following  locations 
(definitions  of  designations  in  parentheses): 

POl  (Primary  stagnation  pressure):  OMEGA  ENGINEERING  PX621, 
101.3-2170  kPa,  15VDC  input  20mA  maximum,  1-5VDC  output.  The 
input  power  was  supplied  from  a  Power-One  HCBB  75W  set  at  15.0VDC 
with  a  1.5 A  capacity. 

PSl,  PS2,  P1-P9  (Mixing  section,  inlet  stream,  static  pressures,  and  wall 
static-pressures  at  variable  locations  from  10  mm  to  436  mm  in  10  mm 
steps).  PSl  (primary  static  pressure)  and  PS2  (secondary  static  pressure): 
Micro  Switch  140PC,  0-207  kPa  differential,  lOVDC  input  at  20mA 
maximum,  1-8VDC  output.  The  input  power  was  supplied  by  an  Acopian 
V59D6A  with  a  capacity  of  lA,  set  at  lO.OVDC. 
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POPR  (Traverse  impact-pressure):  Validyne  DB15,  0-414  kPa  differential, 
connected  to  a  Validyne  CD  12  carrier  demodulator  with  0-lOVDC  output. 

P02  (Secondary  plenums'  average  stagnation  pressure):  Validyne  DB15,  0- 
103  kPa  differential,  connected  to  a  Validyne  CD15  carrier  demodulator 
vk'ith  0-lOVDC  output. 

DP2  (Average  differential  pressure  between  the  secondary  plenums' 
stagnation  and  static  ports):  Validyne  DB15,  0-12.4  kPa  differential, 
connected  to  a  Validyne  CD  15  carrier  demodulator  with  0-lOVDC  output. 

PSPR  (Static  pressure  on  wall  centerline  at  position  of  impact-pressure 
probe  tip):  Validyne  DB15,  0-103  kPa  differential,  connected  to  a  Validyne 
CD  15  carrier  demodulator  with  0-lOVDC  output. 

PS2bk  (Redundant  reading  of  mixing  section,  secondary  inlet,  static 
pressure):  Validyne  DB15,  0-69.0  kPa  differential,  connected  to  a  Validyne 
CD  15  earner  demodulator  with  0-lOVDC  output. 

PAIR  (Supply  air  pressure):  OMEGA  ENGINEERING  PX241,  101.3-2170 
kPa,  15VDC  input  20mA  maximum,  1-5VDC  output.  The  input  power 
was  supplied  from  a  Power-One  HCBB  75W  set  at  15.0VDC  with  a  1.5 A 
capacity. 

The  specifications  of  all  transducers  used,  including  their  calibration  curves,  accuracies, 
and  locations  are  given  in  Appendix  C,  Tables  C.l  and  C.2. 
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and  secondary  flows,  respectively. 

i^^^^^^^'^'^''^  ,„^  ^^^  England  AffUrated  Technologies  23D- 

A  5VDC,  1.0  amp  stepper  motor  (New  Eng 

r^iv^  traverse  mechanism  and  developeo  d 
,,02)  drove  the  rmpact-pressure  probe  tra  e.  ^  ^^_^^^^^,  ^^,   ,  .as 

nf  n  1 3  N-m  at  speeds  of  less  than  380  pulses  pe 
dynamic  torque  of  0.13  M  m       P  ^lewlett-Packard  6266B)  set  at 

.       n  Anvnr   0-5A  variable  power  supply  (Hewlett  fa 
powered  by  a  0-40VDC,  0  5A  ^^^  ^    ^^^  .^^^^^ 

^      ^  1  nA     During  tests  the  stepper  motor  was  run 
12.9VDC  and  1.0 A.     uunng 
.o..         -  24  2  mm/s)  to  minimize  unsteady  and  nonstationary  effects  m 

r:d:lomaLi.o..trnlet.owco„d.t.o     ..^^^^^^^^^^ 

^apezordal-ntotronvelocty  profile.  ^'^ ""      J  „.„„,,,,.  single-ax.s 

„otor  drrver  carf  Ond^tna.  Computer  Source  ^02)  an  ^^  ^^   ^^ 

.epper-motor  control  ..ard  Cdustnal  Conrputer  Sou  c  300^)  m^  ^^^^^^_ 

^     «trnller  cards  were  isolated  from  the  PC  by  usmg 
.epper-motor  dnver  and  controller  cards  ^^^  ^^^^  ^^ 

5VDC,  6.0A  power  supply  (Power-One  HCBB  75W),     P 
data  was  also  used. 

U^^mMMm.^^^^^^^^  ^  ^^^^  ^.„„ent  CO. 

-  —  statrc  pressures  ^^^^^^^^^  „,„ns  to 

t.    1   rKAr.A^\  33M1335)  which  has  30  mercury 
manometer  ban.  (Mode,  33«  ^^^^  ^^__^^  ^^^_^^  ,^^^^  „, 

to  sidewall  were  made  usmg  1/32  m.  ^^^^^^^ 

.     .  ,     1  ™      The  manometer  bank  has  fluorescent  backl.gn 
approximately  3  m.  ,  ,n5oth  second).  Before  each  test  run,  the 

reasonably  fast  photo^aphic  exposure  trmes  (-l/250th  seco 
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o  t^Qt  so  that  thermal 
.       for  at  least  30  mmutes  prior  to  running  a  test, 
lights  were  switched  on  for  at  least 

equilibrium  was  achieved. 

A     „«t^  PressureMeasureiaSSJ& 
Ui^^i^^^'"'^^^^^'''^''^'^^  sealed  gage  .in.ature 

,.„  KuUte  Se^conducto.  (Mode,  ^  ,,„^,,,,e 

„..„ce.  was  cHosen  .o  eU»,na.e  *e  po^  .  ^_^  ^^  _^  ^^  ^^  ^^  ^  ,^„^. 

^,  of  .he  transducer  has  a  max™um<i.amete,      ^  __^ 

-  .0  — e.  «e.  ^'^^::^.,...^...^ 

^="°""^™"T:::;hla.„s«es.e„.honM.a.-g.eKXV..n 

^""=^'^-  "^"dl  ^sduce.  o...de  d.a»e.  and  .he  ,ug  wa„.  ..  was 
the  0.075  mm  gap  between  ^^^ 

The  0  7  mm  pressure  sens.t.ve  area  .^^       P^^      ^^  ^^^^  .^  ^^^^ 

wheatstone  bridge  diffused  mto  a  s.hcon  '""''  ''J"  f  q  ,  52  „„  dianaeter 

,,,„,.„.e.O....«nda.e.Ucree„wh.hh..O^^ 

«ed  through  .t  on  a  0.89  mn,  d.atneter  ctrc  e.  ^^^  ^^^^^^^^  .^  ^^^^^^_ 

..,.0......=  — et^haveara^^^n  ^^_^^^^ 

,ressure,and.nanufacture.sspee~so  —  ^^^^^^  ^^^^  ^^^ 

of  .0.5%  fi.U  scale  best  nts.ra.ghth„e(FSBPSL)  ^^^  „^,„, 

,  .«  .Sn00»P,  and  a  Hat  fre,ue„.  response  ^at  ^^^  ^^^^  ^  ^^^^^^^ 

frequency  of  the  transduce.  ,s  greater  than  2  ,,  ^,  „„.p„. 

f  im  K  to  355  K,  input  impedances  ot  251. 
temperamre  range  of  300  K  to  ^^^^^^.^^  ^^^^  ^^_^^_,.^i„,3  „f 

>.„pedances  of  ,322  O  and  ,,53  .,  .^aU5.  ^^^^  ^^^_^^^^^^  ^^  ,^^ 

>---aand,.2,5.V..resp—  ^^^^^ 
manufacturer  to  standards  traceable  to  the  Nat>o 
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Matched,  custom-built  power  supplies  (see  Fig.  3.13  for  circuit)  were  used  to 
excite  the  transducers  at  18.16  VDC  nominal,  which  is  within  the  20  VDC  maximum 
excitation  level  specified  by  the  manufacturer.  Using  the  linearity  of  the  transducers,  this 
level  of  excitation  raised  the  transducer  sensitivities  to  1.450  nV/Pa  and  1.467  nV/Pa. 

Output  from  each  pressure  transducer  was  split  into  two,  one  being  indirectly  fed 
to  an  Access  multiplexer  (described  above)  for  static  readings,  and  the  other  used  for 
dynamic  readings.  A  custom-built  buffer  was  designed  for  the  split  as  noise  from  the 
multiplexer  was  found  to  seriously  contaminate  the  dynamic  side.  The  buffer  is  shown 
schematically  in  Fig.  3.14,  in  which  the  operational  amplifier  (Harris  HA-5177)  has  an 
ultra-low  offset  voltage  of  20  nV  and  a  low  noise-level  of  9.0  nV/VHz. 

The  dynamic  split  from  the  pressure  transducer  was  fed  to  a  pre-amplifier 
(Stanford  Research  Systems  Model  SR  520)  and  was  AC  coupled  using  the  pre-amplifier's 
0.03  Hz  capacitor,  resulting  in  a  high-pass  filter  with  -10  dB/decade  attenuation.  The 
pre-amplifier  gain  was  calibrated  by  the  manufacturer  to  standards  traceable  to  the 
National  Bureau  of  Standards  (±1%,  DC  to  10  kHz,  and  ±3%,  10  kHz  to  1  MHz),  and  the 
quoted  input-noise  associated  with  the  preamplifier  is  4  nV/v^Hz  for  gains  of  100  or 
greater,  10  nV/v/Hz  for  a  gain  of  50,  and  1 1  nV/v^Hz  for  a  gain  of  20.  The  gain  of  the 
pre-amplifier  was  varied  between  twenty  and  five  hundred  depending  on  the  flow 
conditions.  In  the  AC-coupled  mode,  the  pre-amplifier  has  a  low-pass  filter  capability  of 
-20  or  -40  dB/decade  attenuation,  and  high-  and  band-pass  capability  of  -20  dB/decade 
attenuation.  The  filter  cutofi^  frequencies  range  from  0.03  Hz  to  1  MHz  in  steps  of  3  and 

10. 

The  output  from  the  pre-amplifier  was  AC-coupled  (0.8  Hz  half-power  cutoff 
frequency)  to  a  two-channel,  12-bit  resolution,  Tektronix  Fourier  Analyzer  (Model  2642) 
for  time  and  frequency  domain  analysis.  The  analyzer  has  the  capability  of  measuring 
up  to  1600  spectral  lines  at  a  bandwidth  up  to  200  kHz.  The  bandwidth  was  usually  set 
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The  „o,se  floor  for  .he  analyzer  is  a.  most  4  ,  ^  ^^  ^^^  ^    ^^  ^^^^^^^  ^^^^^^  ^^ 
,„pu,  and  has  a  ^11  scale  ,nput  range  fron>  .1    ^        -  ^^^^^^__^  ^^  ,„ 

a  rate  of  2.56  .imes  the  bandwidth  to  saUsfy  the  Ny, 
__for.heant,— ~^ 

The  analyzer  was  conneCe       ,  ^^^^_  ^  ^^^^  ^^^^^  ,^^  ^„,,,,, 

„5.2k  bits/seoond  RS-232  co„nect,o„.  Us>ng  ^^^_^^ 

,„nle  and  hold,  remote,  high-speed  oau. 
was  used  as  a  simultaneous  sample  and  h  ,,„,,,  and  analyzed  at  the  useA 

u  then  be  subsequently  stored  ana  ^    ; 
The  data  in  memory  could  ^^^^^^  ___  ^^^  ^^^^^^^,  ^,  ^,„,es 

eonveme„ce.Th,s  method  was  used  for  all  da 

continuous  data  records. 

^'^^'^"'''^  employed  to  detect  the  mean  structure  of  the  flow  and  to 

j,„w  visualization  was  employed  accomplished  with 

.,„ate  the  location  of  laminar-to-turbulent  ^^^^.^    ^  ^^,.„,,  ,,. 

continuous  light  source  Schlieren  and  ^^^^^  ...eudently  on  heavy 

•    V     1  1  '^    All  the  components  were  muu 
system  is  shown  m  Fig.  3.15.  AH  m         V  components  could  be 

-  stands  at. he  same  heigh,  as  the  m«..^^- 

The  conrinuous  ligh.  sou  ce  ^^^_^^^  ^_^  ^  ^^_^^_,„  ,,  „„ 

.eunted  in  a  fan-cooled  lamp  ^^  ^^^^  ^^  ^  „„._1  Ump  Power  Supply  and 
diameter,  condenser.  The  light  source  is  powered  y  ^^  ^^  ^^^ 

.,,orCOnel  Model  C-20).   —^  ^^^  °"'°:;i,  optical  Co.).  ..e 

...„„„  F8  25  parabolic  mirror  (J.unem     f 

focus  of  the  firs.  400  mm  diameler.  ^.^^^_  ^,  ,, 

-'•-"rr:::::r;»i  ..--■— - 

focussed  onto  another  knite  eag 
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I.  L  90"  cou.e«,oc^.se  .0.  *e  honzon.,  s«ea„w,se  .reCo.  Mso,       ^.fe 
ed^e  was  removed  to  obtain  a  shadowgraph  photograph. 

■  ■  .  o„  .lass  w,„d„ws  are  schemafcally  shown  mounted  ,n  the  m«,ng 

The  mixing  section  glass  winaowb  m 
toninfi.  3  .6   They  are  made  of  prec.s.on.  opttcalquahty.  BK7-P  crown  glass,  and 
section  in  t-ig.-iio-   my  n  iiR^OA     The 

.ere  mannfaetured  hy  Schott  Glass  Technologtes  to  conform  to  Mn^0-nS3  A. 

,n^  ™m  lonl!  by  50.8  mm  wide  by  12.7  mm  thick.    The 
overall  dtmensions  measure  305  mm  long  by  ^^_^^ 

optical  specfcafons  are  surface  Harness  withm  6  waves  over  the  total  lengt    p 
:Z  1    arc-seconds,  and  scratCdig  of  S0«0.  Mountmg  n,  the  alumrnum  f^es  was 

ipltshed  usmg  a  pourable,  tw.pan  —  rubber  .evcon,  ^U^ 
,SSOO)which  cures  atroomtemperaturccunngshnnkagetsve.  small  atOOOlSm/m, 

and  when  cured,  >.  has  a  tens.le  strength  of  14.5  MPa. 

^  5   f>Qta  r.nllection 

written  in  Quick  Basic. 

N      U.A  ATFRONT  is  used  to  set  up  a  test  run 
A.        The  first  module  (set  of  programs),  called  ATFRON  i. 

for  the  first  time  and  functions  as  follows: 

^,      The  co„d,t,o„s  for  the  test  are  obtained  from  a  predefmed  test  matnx  or 

typed  in  manually. 
A2       The  probe's  axial  position  is  recorded. 
as'.      The  wall  port  corresponding  to  the  impact  probe's  axial  location  is 

recorded. 


57 
A3.       A  set  of  warnings  and  checks  regarding  the  test  rig  were  displayed  to  the 

screen  for  the  operator  to  verify. 
A4.      The  above  data  was  stored  in  a  file  called  ATFRONT.MAB  for  use  by 

other  modules. 

B.  The  second  module  is  called  ATZERO  and  measures  and  records  the  transducer 
zero  offsets  for  subsequent  use: 

Bl.       The  test  rig  is  opened  to  atmosphere  and  75  readings  from  each  transducer 

are  recorded. 
B2.       The  means  and  standard  deviations  are  calculated  and  displayed  for  the 

operator  to  verify. 
B3.       The  means  are  written  to  a  file  called  ATZERO.CAL.   If  the  file  akeady 

exists,  it  is  renamed  ATZERO.OLD  before  being  written. 

C.  The  third  module  is  called  ATMAIN  and  controls  the  actual  test  run.  The 
program  has  the  option  of  controlling  a  test  run  for  an  impact-pressure  traverse  or 
a  wall-pressure.  The  major  details  for  an  impact-pressure  traverse  are  as  follows: 
C 1 .       The  pitot  pressure  probe  is  set  up  for  a  trapezoidal  velocity-profile  traverse 

as  follows: 

1.  Initiate    traverse    using    near-maximum    motor    torque    at    an 
instantaneous  (t==0.00  seconds)  speed  of  400  pulses/second  (pps). 

2.  Accelerate  from  the  top  of  the  mixing  section  at  t=0.00  seconds  to 
953  pps  at  t=0.10  seconds. 

3.  Slew  at  constant  speed  of  953  pps  until  t=1.90  seconds. 

4.  Decelerate    from    953  pps    at    t=1.90  seconds    to    400  pps    at 
t=2.00  seconds. 
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C2. 


5  Instantaneous  stop  at  t=2.00  seconds. 

A    dna  the  keyboard.  Monitoring  of  the  supply 
Solenoid  valves  are  opened  using  the  keyooar 

oLen„.,.n,e...a.o„a„as«.cp.s.e.nses.ea— 
„o„«o.afo.co„ec.flowco„«.ons.T.e.es...U„o.™,fthenow 

_.a.  vawes,  or  .ow„s«ea.  —  vawe  a„o.ea  ...=  co„„o, 

of  the  flow  conditions. 
,       ,fte.  .  00  seconds,  ten  samples  of  each  pressure  and  .en,pe«.u.e  are  read 

,       C--e.s™adewUM.e™„.o.co.o„e.ca.co»...S 

Led-V— s.eas.pc.a„..„co..e.ead.ns...es^- 
p.ocedu.,Uwasfou„d.hat*«e.e,sh.pressu.e.ead.„^-e„b.,„ed 

for  each  counter  reading. 
C3.       Ten  samples  of  each  system  pressure  and  temperature  are  read  at 
sampling  rate  of  1864  samples/second. 

r6        The  gas-supply  solenoids  are  closed. 

LrCisretumedtoitsinitialpositionusingm— motor.tor.ue. 

.he  means  and  standard  deviations  of  the  mixing-section  conditions  are 

.  A   ft.r  the  run  and  are  wntten  to  the  output  file  along 

calculated  before  and  after  the  run,  ano 

with  the  probe  data, 
terminating  the  test. 


C7. 
C8 
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nt  the  control  procedure  is  as  follows: 
Cla     The  Fourier  analyzer  is  set  up 

Caa.     Solencd  valves  are  opened  us.ng  .he  -cey  „„„.,„„,„„„3. 

..,..o„.e.s.eano«p.ec.se.epea*U.     -^^^^^^ 

,•  •  TVip  test  Will  not  run  n  "'"^ 

„„„.o.d..c„.e..owc„„a.o..T^  ^^^^^^^^^^^^ 

— "^^'^"""^  ™l.lve.nowed.e»eeo„.o, 

secondary  valves,  or  downstream  butterfly 


„f  *e  flow  condmons.  ^^  ^^___^^^^^  ^^^  ^^, 

C3a.     After..OOsec„nds..e„san,p.esofeachp  ^ 

a,  a  sanrpling  rare  of  1S64  samples/second,  and  a  checKf 

"  """'  ,  sorted    For  a  100  kHz  bandwidth,  this  step  took 

C4a.     The  Fourier  analyzer  is  started. 

less  than  2  seconds.  ^,  ^ 

CSa.     Ten  samples  of  each  system  pressure  and  temperature  are 

sampling  rate  of  1864  samples/second. 

C6a.     The  gas-supply  solenords  are  closed.  ^^^^_^_^^^  ^^^ 

onH  Standard  deviations  of  the  mixing 
e,a.     ^— ^^^^'^""^^^,,,„„,^,arewrttte„totbeoutp«tf,lea,ong 
calculated  before  and  after  the  run, 

with  the  probe  data.  ^ 

stored  in  a  file  for  post-processing. 
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3.6   Signal  Analysis 

This  section  generically  describes  the  complex,  discrete,  Fourier  transform  routines 
used  by  the  Tektronix  Fourier  Analyzer  (Model  2642).  The  computations  for  the  fast 
Fourier  transform  (FFT)  algorithm  are  performed  in  16-bit,  fixed-point  arithmetic,  and 
uses  a  block  floating  basis  to  scale  the  data.  The  data  obtained  from  the  discrete  Fourier 
transform  (DFT)  of  a  frame  (usually  3200  points  in  this  research)  is  then  converted  from 
fixed  to  floating  point.  Before  summation  of  multiple  frames  for  averaging,  the  DFT  is 
converted  to  32-bit  IEEE  floating-point  format  (8-bit  exponent,  23-bit  plus  sign  mantissa 
plus  the  implied  hidden  bit). 

The  FFT  routine  computes  the  DFT  pair: 

N-l 

Inverse:     z(i)  =  J^  Z(A:)PF*  ,i=0,...J\f-l 

k=o  (3.1) 

A^-l 

Forward:     Z{k)=^z{i)W'^  Jc=0,...JSf-l 

.X)  (3.2) 


where  W=exp(-j2:i/N), 

z(i)  is  a  complex  time  history,  and 

Z(k)  is  a  complex  frequency  function  of  the  same  number  (N)  of  data 

values. 
3.6.1    Autospectrum 

The  single-sided  autospectrum  is  computed  from  the  convolution  of  the  complex 
DFT,  X(k): 

(f2\ 


ui^y 


X\k)X{k) 

(3.3) 
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where  PT=NT=the  length  of  a  data  block,  in  seconds, 

k=0,l,2,...,(N/2)-l  ,  and 

T=  the  sampling  interval. 
The  factor  T^/Pt  arises  because  it  is  omitted  from  the  DFT  computation  for  simplicity  in 
scaling.  The  spectrum  is  single-sided  because  the  textbook  FFT  is  symmetric  about  zero 
and  therefore  has  "negative"  power  frequencies.  To  avoid  the  negative  frequencies,  only 
(N/2)- 1  points  of  X(k)  and  X*(k)  are  used.  This  doubles  the  power  for  the  positive 
frequencies  and  results  in  zero  power  for  negative  ones. 

In  all  of  the  data  collected  in  this  research,  multiple  frames  were  obtained  to 
reduce  the  statistical  uncertainty  of  the  measurements.  Additive  averaging  was  employed 
in  the  following  manner: 

(s   (kV 


'^L(^)-'^U-i)(^) 


M 


(3.4) 


where  Sxx„(k)  is  the  most  recent  autospectrum,  and 

M=total  number  of  frames. 
The  final  result  after  M  averages  is 


^w^  (3.5) 


The  normalized  standard  deviation  (or  standard  error)  of  the  averaged  rms 
measurement  is  given  by 

1 


e=±- 


2v^  (3.6) 
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Most  of  the  results  that  will  be  presented  are  averages  of  10  frames  which  have  a 
standard  error  of  16%.  However,  some  of  the  data  was  compared  to  40-frames  averages 
which  have  a  standard  error  of  8%,  and  it  was  found  that  10  frames  gave  repeatable 
results  in  terms  of  both  frequency  amplitude  and  rms.  These  results  will  be  discussed  in 
Chapter  5. 

3.6.2   Cross-spectrum 

The  cross-spectrum  is  computed  in  a  similar  manner  to  the  autospectrum  except 
that  the  cross-spectrum  is  complex.   The  cross-spectrum  of  a  size  N  frame  is 

''^^    P,  (3.7) 


and  the  average  cross-spectrum  of  M  frames  is  calculated  by 


(S^ik) 


M    ) 


(3.8) 


where  S^(k)  is  the  most  recent  cross-spectrum.    The  final  resuh  of  the  M  averages  is 
thus: 

.,  Af-l 

S'(k)=—yS    (k) 
^^  ^   M^   ^^  ^  (3.9) 

3.6.3   Correlation  Functions 

The  auto-  and  cross-correlation  are  covariance  functions  and  are  computed  using 
inverse  Fourier  transforms  (EFTs).  For  accurate  covariance  results,  the  acquired  data  is 
time-averaged,  zero-padded  to  avoid  wraparound  error,  and  has  no  window  applied  to  the 
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spectrum  oalculaUon^  The  speetrum  .s  modified  by  changmg  .he  s.gn  of  the  odd-indexed 
data  values.  k=l,3.5,...,(N/2)-l,  followed  by  an  IFT  to  obtam  the  covariance: 
.  Cross  -correlation :       s^(i)  =IFT\S^{k)] 
Auto  -correlation :       s^{i)  =IFT]S^ik)] 

where  the  subscnpt  "c"  distinguishes  these  ftinctions  from  the  fmal  corrected  Sanctions. 

For  both  functions,  the  lag  index,  i,  runs  from 

(3.12) 

The  auto-correlation  function  .  symmetnc,  and  hence  there  are  only  N  d.stmct  values  for 

this  function. 

Both  correlations  are  then  adjusted  to  account  for  the  linear  taper  that  is  induced 

by  the  fact  that  the  spectra  are  divided  by  N  instead  of  N  minus  the  lag  number.    Hie 

statistical  result  is  to  provide  an  unbiased  estimate  of  the  correlation  function: 


^^(i) 


N 


[N^] 


W(0   '  --f>--l'0'l'-'T-l 


(3.13) 


^64   Transfer  ^"d  Toherence  Functions 

The  averaged  transfer  function  is  computed  from 

^^  ^    S'Jk) 
and  the  averaged  coherence  function  is  computed  from 


(3.14) 


5"  (k)  S'Jk)       ,    -.      N  . 
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The  transfer  function  compares  the  gain  of  the  "response"   channel  to  the 

"reference"  channel,  that  here  would  be  calculated  from  two  pressure  transducers'  outputs. 

It  also  contains  phase  information.    Two  different  results  can  be  obtained  by  switching 

the  response  and  reference  channels,  and  the  results  can  be  compared  to  identify  which 

direction  certain  frequencies  are  travelling. 

The  coherence  function  is  similar  to  the  transfer  function  except  that  it  is  not 

affected  by  switching  the  reference  and  response  channels,  and  contains  no  phase 

information.  It  is  useful  for  identifying  the  relative  power  of  frequencies  present  at  two 

locations. 

3.6.5   Windows 

For  all  measurements  except  for  the  correlation  measurements,  a  window  is 
applied  to  the  data  in  the  frequency  domain.  The  benefit  of  using  a  window  is  to  reduce 
the  effects  of  discontinuities  in  the  analyzed  signal  introduced  by  the  FFT  process.  The 
function  of  the  window  is  to  reduce  the  values  of  the  ends  of  a  N  point  frame  to  zero,  and 
avoid  wraparound  error.  No  windowing  is  used  for  the  correlation  measurements  because 
the  window  is  applied  in  the  frequency  domain,  and  would  therefore  distort  the  time 
record,  producing  an  incorrect  analysis. 

The  Manning  window  is  applied  to  the  data  when  a  window  is  used,  and  is 
implemented  for  each  FFT  operation  in  a  calculation: 

Z'(k)=^Z(k)-^Z(k-l)  (^-^^^ 


where  Z'(k)  is  the  windowed  FFT. 
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It  is  obvious  that,  compared  to  using  no  window,  windowing  affects  the  rms  power 
and  amplitudes  of  the  spectral  lines.  This  can  be  overcome  by  applying  an  amplitude 
scaling  factor  of  2  (reciprocal  of  the  first  weighting  term),  or  an  rms  scaling  factor  of 
2.667  to  the  windowed  data.  In  this  research,  the  rms  scaling  factor  is  always  used  so 
that  the  total  power  values  are  correct.  As  amplitude  correction  is  not  employed,  the 
effect  on  spectral  amplitude  on  a  pure  sine  wave  for  example,  would  be  an  amplitude 
reduction  of  1.22  compared  to  the  true  value. 

A  further  effect  of  the  Hanning  window  is  to  reduce  the  effective  duration  of  the 
record  being  transformed  to  half  of  that  represented  by  the  samples  themselves.  This 
reduction  can  be  overcome  by  overlapping  the  frames  by  50%,  and  is  used  here. 

3.7  Boundary-Layer  Conditions  at  the  Mixing-Section  Inlet 
The  boundary  layer  thicknesses  at  the  mixing-section  inlet  were  calculated  to 
quantify  the  mean  inlet-conditions  and  to  validate  the  impact-pressure  measurement.  The 
thicknesses  were  calculated  using  the  compressible  Thwaites  method  (Thwaites,  1949)  for 
momentum  thickness  of  Rott  and  Crabtree  (1952),  that  is  valid  for  the  compressible 
boundary  under  arbitrary  pressure  gradient  with  no  imminent  separation.  The  necessary 
assumptions  are  :  (a)  adiabatic  wall,  (b)  the  Prandtl  number  of  the  gas  is  unity,  and  (c) 
the  viscosity  ji  is  a  linear  function  of  temperature. 

These  calculations  were  performed  for  the  trailing  edge  thicknesses  of  the  primary 
and  secondary  on  either  side  of  the  nozzle  block.  The  momentum  thicknesses  were 
assumed  to  be  zero  at  the  start  of  the  convergence  sections  leading  into  the  mixing 
section,  and  the  calculations  are  shown  in  Table  3.1.  The  details  of  the  method  are  given 
in  Appendix  B. 
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The  initial  thicknesses  on  the  side  walls  were  not  calculated,  but  are  approximately 
the  same  as  the  ones  calculated  presented  in  Table  3.1  due  to  the  similarity  in  the 
freestream,  axial  Mach  number  profiles. 

Using  the  calculated  value  for  6*  to  get  the  effective  primary  flow  height  2(b-6*), 
the  predicted  primary  Mach  number  can  be  estimated  using  this  value  in  the  one- 
dimensional,  isentropic,  gas-dynamics  relation  for  area  ratio.  Also,  it  can  be  calculated 
that  the  measured  value  should  be  increased  by  about  0.2%  to  account  for  the  static 
temperature  and  pressure  measurements  being  made  on  the  side  wall  of  the  primary 
plenum.  This  results  in  a  corrected,  measured  primary  Mach  number  of  1.79.  This 
agrees  with  the  primary  Mach  number  calculated  from  the  measured  nozzle-throat  and 
exit  heights,  and  accounting  for  the  displacement  thicknesses  in  Table  3.1,  to  within  the 
tolerances  of  the  height  measurements. 
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Table  3.1    Boundary-layer  thicknesses  of  primary  and  secondary  supplies. 


Location 

d{mm) 

8'  {mm) 

8  (mm) 

8 

e 
s 

e 

Primary  Throat 

0.0085 

0.027 

0.051 

0.54 

0.17 

3.2 

Primary  Trailing  Edge 

0.0085 

0.045 

0.073 

0.62 

0.12 

5.3 

Secondary  Trailing 
Edge — Inner  Wall 

0.037 

0.099 

0.20 

0.50 

0.19 

2.7 

Secondary  Outer  Wall 

0.039 

0.103 

0.21 

0.50 

0.19 

2.7 
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Figure  3.3        Flow  conditioning  sections,  settling  chambers  and  nozzles.   Flow  is 
from  right  to  left. 
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—  16.15  mm 


Figure  3.4       Nozzle-block  dimensions. 
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Figure  3.5        Mixing  section  and  impact-pressure  probe/pressure  transducer/stepper- 
motor  assembly.   Mixing-section  extension  is  not  shown. 
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Figure  3.6       Wall  static-pressure  plate  used  for  manometer  and  strain-gauge 
transducer  measurements. 
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Figure  3.7       Brass  transducer-sleeves  and  plugs  used  for  time-accurate  wall-pressure 
measurements. 
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Figure  3.9       Top  rail  of  mixing  section  showing  fitted  plugs. 
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Figure  3.10      Plugs  for  top  and  bottom  rails  of  the  mixing  section. 
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Figure  3.11      Exhaust  system. 
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Figure  3.12      Schematic  of  instrumentation  and  computer  interface. 
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Figure  3.13       Custom-built  power  supply  for  Kulite  transducer. 
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CHAPTER  4 
EXPERIMENTAL  RESULTS-mean  MEASUREMENTS 


4.1    Experimental  Objectives 

The  primary  objective  of  the  experimental  investigation  was  the  measurement  of 
time-accurate  wall-pressures  of  the  potential  and  shearing  flows  in  a  two-dimensional, 
single-nozzle,  supersonic  ejector  discharging  into  a  constant-area  mixing  tube.  Using 
these  measurements,  a  nonintrusive  technique  for  shear-layer  detection  and  measurement 
is  developed.  However,  in  order  to  validate  the  time-accurate  measurements,  and  to  gain 
insight  into  the  flowfield,  mean  measurements  were  taken. 

To  determine  the  mean  flowfield,  Schlieren  flow  visualization  was  employed,  as 
well  as  measuring  cross-sectional  impact-pressure  profiles  with  a  traversing  probe  at  a 
number  of  streamwise  locations  of  the  mixing  tube.  Mean,  static  wall-pressures  were 
measured  using  a  manometer  bank.   These  results  are  presented  in  this  chapter. 

4.2  Test  Conditions 
The  complexity  and  difficulty  of  measuring  and  interpreting  the  turbulent  wall- 
pressure  in  a  compressible  flow  was  reason  enough  to  choose  a  single  flow  condition. 
From  previous  experiments  reported  in  the  literature,  it  is  known  that  the  mean  flowfield 
is  affected  by  at  least  the  following  parameters: 

(a)  primary  and  secondary  molecular  weights 

(b)  primary  and  secondary  specific  heats  ratios, 

(c)  velocity  ratio,  r- 


(c)        density  ratio,  s  =— 

Pi 
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(d)  convective  Mach  number,  A/. 

(e)  mass  flow  ratio,    E^  =— 

A  nominal  primary  Mach  number  of  Mi=1.8  was  chosen,  which  is  near  the  limit 
for  the  available  supply  pressure  of  the  facility  for  atmospheric  blowdown.  The 
secondary  Mach  number  was  chosen  to  be  Mj^O.SO.  The  primary  and  secondary  static 
pressures  at  the  nozzle  exit  were  matched  at  103.5  kPa  to  minimize  leakage  from  the  test 
rig  to  the  atmosphere.  Another  advantage  of  choosing  the  static  discharge  pressure  to  be 
just  above  atmospheric  compared  to  sub-atmospheric,  is  that  the  range  of  turbulent  scales 
present  in  the  flow  is  increased  due  to  the  increased  Reynolds  number. 

Multiple  experiments,  which  are  described  in  the  following  sections,  were 
conducted  at  a  single  flow  condition,  and  parameters  of  the  flow  are  given  in  Tables  4.1 
and  4.2.   Calculation  procedures  for  these  values  are  given  in  Appendix  B. 

4.3    Optical  Visualizations 

Schlieren  photographs  of  the  mixing  section  with  the  knife  edge  at  90°,  60°,  30°, 
and  0°  are  shown  in  Figs.  4.1,  4.2,  4.3,  and  4.4,  respectively,  and  shadowgraphs  are 
presented  in  Fig.  4.5.  The  photographs  show  the  flow  from  the  throat  of  the  primary 
nozzle  to  just  over  27  cm  (x/2b»21.3)  as  indicated  by  the  scale  in  each  photograph.  The 
scale  zero  is  positioned  at  the  nozzle  exit. 

As  the  visualizations  show  the  flow  upstream  of  the  nozzle  exit,  the  mixing-section 
inlet  static  pressure  could  not  be  measured.  Therefore,  the  flow  conditions  were  set  using 
the  streams'  stagnation  pressure  and  the  mixing-section  exit  static  pressure.  These 
pressures  were  known  from  previous  runs  where  it  was  possible  to  measure  the  mixing- 
section  inlet  pressure.  The  stagnation  pressures  were  always  matched  to  the  previously- 
recorded  ones,  but  the  mixing-section  exit  pressure  was  always  about  3.5%  lower.  The 
lower  value  can  be  attributed  to  the  differing  wall  friction  of  the  optical  window  plates 
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(the  glass  mounting-compound,  particularly)  as  compared  to  the  aluminum  or  brass  plates. 
However,  it  is  believed  that  this  made  little  difference  to  the  flowfield. 

4.3.1  Waves  in  the  Flow 

Compression  and  expansion  waves  originating  from  the  throat  of  the  nozzle  are 
visible  in  all  the  photographs.  Ideally,  fmite-strength  waves  would  not  be  generated  from 
the  nozzle  throat,  but  there  is  no  way  to  avoid  finite  waves  since  the  flow  is  not  ideal. 
These  waves  impinge  and  reflect  from  the  upper  and  lower  shear-layer  as  they  propagate 
downstream.  The  shadowgraph  (Fig.  4.5)  shows  that  these  extend  to  approximately 
7.5  cm  (x/2b«6)  downstream,  whereupon  elliptical  structures  appear  along  the  centerline 
of  the  flow  (Fig.  4. 1).  Measurements  of  the  static  wall-pressure  that  are  presented  later 
in  this  chapter,  show  that  the  waves  cause  a  maximum  static-pressure  change  of 
AP/P"0.03  in  this  region,  and  it  can  therefore  be  concluded  that  they  are  weak. 

4.3.2  Growth  Rates 

The  visual  growth-rate  for  the  shear  layer  can  be  determined  firom  the  Schlieren 
photographs  with  the  knife  edge  at  90°.  The  precision  of  this  method  is  subject  to  the 
sensitivity  of  the  Schlieren  system  and  the  local  refractive-index  gradient  of  the  shear 
layer.  The  refractive  index  is  a  linear  function  of  the  local  flow-density  and  is  affected 
by  the  local  gas-composition  and  the  mixture's  thermodynamic  state,  and  hence  the 
Reynolds  number.  In  this  case,  it  is  felt  that  the  stream  Reynolds  numbers  are  high 
enough  to  permit  good  accuracy  in  determining  the  visual  growth-rate. 

The  growth  rate  was  obtained  by  fitting  straight  lines  to  the  mean  tangents  of  the 
upper  shear-layer  extending  from  where  growth  starts,  Xgf«10  mm  (Xg/2b«0.8),  to  where 
the  secondary  edge  of  the  shear  layer  does  not  grow  into  the  secondary  as  fast,  x«75  mm 
(x/2b«6).  Averaging  photographs  of  three  different  runs  results  in  a  visual  growth-rate 
of 
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(■^"^o)  (4.1) 


with  very  little  spread  (<2%)  between  the  individual  measurements.  The  accuracy  of  the 
measurements  is  estimated  to  be  ±10%. 

The  flattening  out  of  the  shear-layer/secondary  boundary  downstream  of  X"75  mm 
(x/2b«6)  implies  that  the  shear-layer  starts  to  undergo  a  structural  change  at  this  point. 
This  suggests  that  the  change  in  boundary  conditions,  manifested  by  the  large  elliptical 
cells  located  on  the  centerline,  are  the  cause  for  the  growth-rate  reduction.  This  will  be 
considered  further  in  the  next  section. 

4.3.3   Visible  Structures 

Several  structures  can  be  seen  in  the  visualizations.  The  ones  in  the  primary  core 
have  already  been  mentioned  above.  They  are  structures  located  between  approximately 
80  mm  (x/2b«6.3)  and  130  mm  (x/2b«10.3)  and  appear  in  Figs.  4.1  -  4.5;  they  are 
elliptical  in  shape  with  their  major  axis  along  the  centerline.  It  is  proposed  that  these 
structures  are  cells  emanating  from  the  viscous  interaction  of  the  wave  structure  in  the 
primary  core  with  the  shear-layer.  Their  distinct  elliptical  shape  suggests  that  the 
boundary  may  bound  large  eddies  in  the  shear-layer. 

In  Fig.  4.4,  it  can  be  seen  that  the  waves  appear  to  penetrate  the  shear  layer 
starting  at  x«55  mm  (x/2b»4.3).  However,  the  wave  cannot  be  penetrating  the  shear  layer 
as  the  width  of  the  penetration  is  thicker  than  the  wave,  and  the  penetration  extends  into 
the  subsonic  region.  This  must  be  the  result  of  the  waves  inducing  vorticity  in  the  shear 
layer. 

In  Figs.  4.3  and  4.4,  a  thickened  wave  structure  can  be  seen  centered  at  x«65  mm 
(x/2b«5.1)  which  extends  up  out  of  the  shear  layer.    This  suggests  that  secondary  fluid 
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is  being  entrained  at  a  locally  increased  rate  such  as  has  been  observed  in  underexpanded 
jets  (Gutmark  et  al.,  1990).  Also,  Tarn  (1972)  in  a  Mach  1.5  axisymmetric  jet,  and 
Gutmark  et  al.  (1990)  in  a  Mach  1.5  rectangular  jet  both  observed  that  one  of  two  distinct 
noise  source  locations  is  located  at  x/2b«=6.  Tarn  proposed  that  the  noise  source  is  the 
result  of  large-scale,  long-wave  instabilities  with  scale  length  equal  to  the  jet  diameter, 
which  may  lead  to  rapid  destruction  of  the  jet.  This  is  in  agreement  with  the  present 
observations,  as  within  15  mm  (x/2b»1.2)  the  waves  have  disappeared  from  the  core  flow. 
This  leads  to  a  possible  conclusion  that  the  energy  from  the  finite-strength  waves  is 
transferred  to  the  shear  layer  or  radiated  as  sound,  and  the  process  is  the  precursor  to  the 
formation  of  the  elliptical  structures  discussed  above. 

From  Figs.  4.3(a)  and  4.4(a),  it  is  possible  to  measure  the  mean  spacing  of  the 
wavy  structures  €  before  the  end  of  the  primary  core.  Similar  measurements  have  been 
made  in  single,  subsonic  shear  layers  by  several  investigators.  Brown  and  Roshko  (1974) 
found  that  the  ratio  €/6^  was  approximately  3  over  a  wide  range  of  velocity  and  density 
ratios.  Spencer  and  Jones  (1971)  found  the  most  frequent  spacing  (found  to  be  almost 
equal  to  €  by  Brown  and  Roshko)  was  ^Jb^=5A,  while  Koochesfahani  et  al.  (1979)  found 
£/b^=3.9.  Papamoschou  (1986)  found  with  an  estimated  accuracy  of  25%,  that  for 
compressible,  heterogeneous  shear-layers,  £/bpn  varied  from  1 .2  in  the  highly  compressible 
case  (Mj=1.5)  to  a  high  of  2.5  for  a  case  with  Mj«0.65.  He  suggested  that  for  the 
incompressible  shear  layer,  bjbp^=0.60  translating  into  €/8pi,=1.8  for  Brown  and  Roshko's 
data  and  €/6pi,=3.1  for  Spencer  and  Jones'  data. 

The  measured  value  from  the  visualizations  for  the  present  case  is  €/6pi,«2.3,  which 
was  determined  by  averaging  several  values  at  different  streamwise  locations  and  for 
different  test  runs.   The  accuracy  of  the  above  value  is  estimated  to  be  ±30%. 
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Downstream  of  x=140  mm  (x/2b«ll)  to  about  x=230  mm  (x/2b«18).  Fig.  4.4(b) 
shows  three  light  circular  shapes  of  size  on  the  order  of  half  the  mixing-section  height. 
This  suggests  that  an  increase  in  scale  occurs  in  this  region. 

4.3.4   Laminar-to-Turbulent  Transition 

The  flow  along  the  trailing  edges  of  the  nozzle  is  believed  to  have  been  laminar 
based  on  the  calculation  results  presented  in  Table  3.1.  These  results  can  be  written  in 
terms  of  the  trailing  edge  boundary-layer  momentum  thicknesses  and  are  presented  in 
Table  4.3.  The  Reynolds  number  based  on  trailing  edge  momentum  thickness  and 
freestream  conditions  is  ReQ=pU9/[x.  Also  presented  is  the  Reynolds  number  based  on 
shear-layer  thickness  at  the  start  of  visual  growth,  ReMp^^USAlii^,  and  the  number 
of  momentum  thicknesses  to  visual  growth. 

The  values  in  Table  4.3  for  the  primary  are  in  excellent  agreement  with 
Bradshaw's  (1966)  criterion  that  distance  to  full  development  in  a  free-shear  layer  takes 
a  total  distance  of  VxlO'vAJj  for  the  range  50Q<Reg  <1000.  The  value  of  Re^  is  also  in 
good  agreement  with  Goebel  and  Button's  (1990)  criterion  that  the  fiiUy  developed  value 
for  this  Reynolds  number  is  10'. 

4.4   Static  Wall-Pressure 

The  static  wall-pressure  was  measured  by  three  methods:  with  a  30-column, 
mercury-filled  manometer,  with  strain-gage  pressure  transducers,  and  with  a  Kulite 
miniature  pressure  transducer,  all  of  which  are  described  in  Chapter  3. 

Measurements  were  made  on  both  the  left  and  right  walls  with  the  manometer,  and 
on  the  left  wall  for  the  other  two  measurements.  The  manometer  measurements  were 
recorded  by  taking  photographs  of  the  manometer  bank.  Approximately  fifteen 
photographs  (3  or  4  for  each  test  run,  hence  4  or  5  runs/plate-location)  were  averaged  for 
each  measurement,  and  the  95%  confidence  intervals  were  calculated.      To  make 
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measurements  covering  the  full  length  of  the  mixing-section,  it  was  necessary  to  position 
the  measurement  plate  in  two  positions  along  the  wall  on  each  side  of  the  mixing-section. 
Hence,  four  sets  of  measurements  were  taken  to  arrive  at  a  complete  set  of  measurements. 
The  two  locations  on  one  side  of  the  wall  produce  an  overlapping  measurement  region 
between  x/2b=13.9  and  20.6.  When  the  measurement  plate  is  positioned  upstream,  it  ends 
at  x/2b=25.5;  when  positioned  downstream,  it  starts  at  x/2b=13.9.  Figure  4.6  shows  a 
typical  photograph  of  a  manometer  reading. 

The  strain-gauge  transducer  measurements  were  taken  from  sets  of  ten  readings 
made  during  the  separate  runs  for  impact-pressure  traverses  (the  traverses  are  described 
in  the  next  section),  and  were  located  on  the  wall  at  y/2b=0  of  the  traverse.  The  Kulite 
transducer  measurements  were  also  taken  during  different  runs  and  averaged  from  ten 
readings.  The  short,  brass,  plate  sections  used  for  the  Kulite  transducer  measurements  are 
different  from  the  long  aluminum  plates  used  for  the  other  two  types  of  measurements, 
but  allow  measurements  at  different  values  of  y/2b. 

Figure  4.7  shows  a  comparison  of  the  three  types  of  measurements  made  on  the 
left  wall  of  the  mixing  section.  There  is  good  agreement  between  the  manometer  and 
Kulite  measurements  at  y/2b=-1.60  and  y/2b=1.08,  but  more  discrepancy  at  y/2b=0, 
especially  up  to  x/2b=20.  This  distance  approximately  corresponds  to  the  length  of  the 
supersonic  region  along  the  centerline.  The  smaller  amplitude  waviness  of  the  Kulite 
centerline  measurements  in  this  region  suggests  that  the  aluminum  wall  used  for  the 
manometer  and  strain-gauge  measurements  introduced  waves,  probably  caused  by  a  flow 
interaction  with  the  relatively  large  1.6  mm  diameter  pressure  taps.  In  comparison,  the 
Kulite  transducer  has  multiple  holes  of  0.15  mm  diameter,  and  the  gap  around  the 
transducer  is  even  smaller.  The  Kulite  measurement  values  along  the  centerline  are  lower 
than  those  away  fi'om  the  centerline  up  to  the  end  of  the  supersonic  region,  and 
correspond  to  the  primary-core  structure  observed  in  the  visualizations.     The  wide 
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variations  in  the  measurements  at  x/2b=7  in  Fig.4.7(a)  correspond  to  the  initial  location 
of  the  large  elliptical  cells  discussed  in  Section  4.3.2. 

The  manometer  curve  in  Fig.  4.7(a)  shows  what  appears  to  be  a  distinct  wave 
pattern  from  x/2b=17  to  21  if  the  points  are  joined  with  straight  lines.  This  region  is 
located  where  the  pressure  plate  overlaps  when  positioned  in  the  upstream  and 
downstream  locations.  It  can  be  concluded  then,  that  these  waves  are  not  real  product  of 
the  ejector  flow,  but  are  a  manifestation  of  the  wave  production  and  pressure  loss 
associated  with  the  transition  of  the  flow  from  one  plate  to  the  other. 

The  difference  between  the  strain-gauge  and  manometer  readings  are  shown  in  Fig. 
4.7(b).  The  waviness  suggested  in  the  supersonic  region  is  probably  caused  in  the  same 
manner  as  explained  above  for  the  manometer  reading.  There  is  little  diff^erence  between 
the  readings  except  that  the  streamwise  pressure-gradient  appears  to  start  dropping  more 
rapidly  for  the  strain-gauge  measurements  at  x/2b"29.5.  This  suggests  that  the  impact 
probe  may  have  more  of  an  effect  on  the  pressure  field  downstream  of  this  location,  but 
little  effect  upstream. 

Figure  4.8  shows  a  comparison  of  the  pressures  on  the  left  and  right  walls  of  the 
mixing  section.  The  very  close  agreement  between  the  measurements  on  the  two  walls 
suggests  that  the  static  pressure  field  is  very  nearly  two-dimensional. 

4.5   Impact-Pressure  Measurements 
4.5.1   Normalized  Profiles 

Impact-pressure  profiles  were  recorded  at  twenty-five  streamwise  locations  in  two 
ranges  which  correspond  to  the  upstream  and  downstream  positions  of  the  measurement 
rails.  The  ranges  are  from  x/2b=0  to  20.55,  and  from  18.66  to  34.47  in  steps  of 
1.581x/2b.   The  data  was  recorded  in  steps  of  Ay/2b=0.010. 
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P„„e  4.9  show,  so^e  of  .he  i-npact  traverses.  F.gure  4.9(a)  shows  da.  for 
._s  whe.  .he  pn.a.  eo.  .s  a,.os.  u„.s.u.he.  exeep.  fo,  *e  wea.  waves 
aescnheaabove.He.easa,„,*ewaves.e„^hca„he<,ua„.meaby.e,a,ses.oh.e„e 

^  /P  .(^ax)  =0.975  .n  .he  core.  W,.h  no  s.«c  pressu.  ehange.  .h.s  corresponds  .o  a 
2^1.  L  p.or,.es  a.  ve.  sn.oo.  eve^he.e  e,se  wh.eh  de—.es  «  • . 
„,.  now  cuaU.  .s  h,sh.  When  *e  ..ave.es  a.e  .epea.ed  fo.  .he  profiles  .  P..  9(a), 
tt,e  resul.  are  .denUca,  .o  wi.h,n  U,e  accuracy  of  U.e  measuremen.  .echn„ue. 

TT^e  profiles  are  shgh.ly  shifted  .0  *e  lower  side  of  .he  n,™ng  seCon.  wh.ch 
also  no..ceaUe  ,n  .he  Schl.eren  pho.o^aphs  by  observes  .ha.  *e  upper  and  lower 
_,,,.e„..al-fiowsd.sappeara.d,fferen.loca.,o„s.Ananon,a,y~^^ 
(,;.^ve)  and  lower  (y/2Wve)  shear  laye.  ,s  seen  on  *e  inle.  profile,  x/2h=0,     .he 
p„„.yWshear-,ayer  houndanes.  The  cause  of  .h,s  ,s  a..nhu.ed  .o  d.fferences  ,n  .he 
ra.lin..ed,e  boundary-layer  flow  of  .he  no.le.  This  ,s  supponed  by  d,e  d.sappearance 
:;  .:  anonra.  before  .he  nex.  .averse  a.  ..b=..S  where,  fro.  fl.e  Schl.eren 
pho.ographs.  .ransi.ion..o-«.rbu,ence  .s  known  .o  have  already  occurred. 

Figure  4  9(b)  shows  .averses  downs.eam  of  *e  end  of  .he  priman,  po.enua,  core. 

These  look  dis.,nc.ly  "bu^py"  compared  .o  .hose  u^.rean,,  wh.ch  sugges.  .he  presence 
of  ..ge.  Slow  n,ov,ng,fiu,ds.ruc.ures.  They  n.us.be.argeandslow„ovu,gbecause.he 

.es.,nse-..n.e  for  .e  n,easure„e„.  ,s  no  bene,  .ban  .he  ,000  Hz  read.ng  fre,ue.y  o 
.He  pressure  .ransducer,  and  .he  larges.  .bumps'  appear  .o  have  a  charac.ens.,c  ,eng.h  of 
_a  5  .n,  and  a  charaCensUc  .,n,e  based  on  d,e  probe  sp«d  of  -0.1  s.  Thrs  charac.er,s.,c 

•  .•      f  tv,^  tnie  flow  and  therefore  these  measurements 
time  is  too  long  to  be  charactenstic  of  the  true  flow, 

nt.     However  they  are  representative  of  the  flow,  and 
are  not  true  mean  measurements.    However,  tney 

useful  information  can  be  obtained  from  them. 

TO  see  how  representative  these  "bumpy"  measurements  are,  multiple  traverses  at 
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H  sen  because  *ev  d.sp.ayed  the  leas.  repea.bn..y  of  a,,  the  trave.es 
..at.o„.  «.e  Chose    -^  ^^^  ^^^^^_^^  ^^  ^^,^,3  ,  ,,  , 

measure.    Ftsure  4, 0(a). hw^^_^^^^^,^,.,^ 

also  the  stream  locatton  v,here  the  shear  lay  ^  ^ 

.,       Hhv  the  tails  ofthe  traverses  having  values  of  P„/Pp.Cmax) 
as  evidenced  by  the  ta.ls  .        .  p  m  (max)-0-24.  It  can  be  seen  that 

V    ■  1  ,  .„  tlie  mixins  section  of  e^e,A"'^>^ 
,0  the  value  at  the  inlet  to  the  mK.ng  ^^^  ^^ 

.      r      ,  difference  between  the  two  traverses,  and  suggests  m 

*-  '=  '  ^■^■''^"'  ""^^7     bably  three-dintensiona.   Figure  4  ,0(b)  shows  three 
«,is   location  is  unsteady  and  probably  three  ^  ^.^^  ^ 

„s  at  ..b=.0.  that  ,s  the  ^^^  -rep..^  _  ^^^^^^    ^^ 

repeatability  for  this  case  ,s  far  better  than  for  xy2b=,5.  ^^^^  ^^ 

i^c  reoeatability  than  the  upper  layer,  suggesx    g 

'-  ^r:z:zz^  -*  ■"  -  — ^"'  "^  '^^ "" 

"~  Coin  wher  .e  shear  layers  start  to  interact  with  the  waUs. 
~rrr     l-,ng.hetraverseresultstootherco..onno»P— 

-;;:r::::r::cc:rrr-. 


described  in  Appendix  B. 
/^  «;  9    Simi'^"^  Profiles 


l^^^-^^SlilarifiarolilSS  ^^.tmcted  for  the  traverses 

„here  the  pnmary  potential  core  still  exists.  ^ 

...ity  profiles  IS  that  .e  assumpUon  that  the  stahc  P^  ^  ^  '^^  „, 

.  ,e  center  of  the  ..ing  .be  is  not  always  -^  ^        4  4.    It  is  also 

This  has  already  been  demonstrated  in  Section  4.4 
finite-strength  waves.    Uis  has  aire    y  ,„,  velocity  can  be  compared 

,nter  the  argument  that  similanty  profiles  for  velocity 
necessary  to  counter  the  argum  ,,^sformation.    These 

to  similarity  solutions  derived  using  an  Illing^orth-Stewartso 
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„e40ds  canno.  compare  .o  a  solution  of  .he  Navier-S.o.es  e,ua.io„s,  Cosed  «,th  a  s.te- 
„,.He-an  .u.bu,e„ce  ™ode..  and  computed  on  a  fas.  d,g,.a,  co™pu.er^  Fro.  ,he  con,pu.er 

„ind  fta.  fte  impac-pressure  similar,.,  solutions  are  presen.ed  here. 

Figure  4.12  shows  a  plo.  of  .he  s,m«a,i.y  parame.ers  .  agains.  n  *a.  are  defined 


as 


^pit   ^pU, 95% 

'^^'^ Ip  (4-3) 


where. 


and, 


>';c%=>'fp.'^%) 


P       =JL\p  ,(max)-?Jmin)l+P  /min) 

f        v/oK_n  tn  9  49     The  raw  data  profiles  have  been 
and  are  plotted  for  traverses  from  x/2b=0  to  9.4y. 

u       1      .fP  /P   rmax^=0  5    The  error  function  curves  shown  in  the 
centered  where  the  value  of  Ppi/Fpi.(max;-u.:j. 

figure  are  given  by 

7r(7])=0.5+0.55e7/(T2.527j+1.13) 

The  las.  ^averse  ,s  after  *e  prima,.  po.en.,al-core  has  ended  and.  al*ough  i.  follows  «,e 
error  curve,  ,.  ,s  no.  as  smooft  as  the  upstream  profiles. 
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„  .  .„.ere..„s  .o  no.  *.  *e  spread  .n  the  s.^Uanty  prom.  fo.  *e  upp.. 

«     ,„„    This  is  contrary  to  intuiuon  as  one 
generally  on  the  postoe  side  of  the  error  function.  Th,s.  con    ty 

Lm  expect  the  profiles  on  either  side  of  the  centerl.e  to  he  .irror  i.ase   of  each 
.her.  Th,s  can  he  explained  in  temis  of  the  vorticty  thic.ess  .owth-rate,  denned 

"    dU\  (4.7) 


,^he,nferredfro..he,™pact-pressnreprofilesthatthevort,citythic..ess.o..-.. 
.  .he  lower  shear-layer  at  ..«  and  1.3S  iS  less  than  that  for  the  .averses  ..^ 
aownstrea..     .  the  upper  layer  at  ..«..,  there  Is  a  distinct  d..nrhance  «. 

propagates  weU  into  *e  up^  secondary-no.   In  the  sanie  region  in  .e  lower  1  y. 
:eearea,sod„rerencesin*evortic.ty.hiC.esscon,paredtothee.ror  .notion  andthe 

aownstrean^proriles.  Hence,.tca„heco„cludedthatvor.ci.ySenerated,ro™.he    a.  n 
.,,e  of  the  no.le  results  in  nonsunilanty  of  the  velocity  profiles  at  leas,  up  o  x^^l>^. 
^is  distance  iS  .ust  over  four  ti.es  longer  than  x.  and  agrees  wi.  the  LOV 
measurements  of  Ooehel  and  Outton  (19.0)  who  found  equivalent  values  of  this  ratio 

3  for  M,=0.46,  and  6  for  M,=0.69. 

4  S  3   Growth-Rates 
The  edges  of  *e  shear  layers,  wh.ch  are  defined  at  the  posH.ons  y,,.  and  y,.. 

.own  in...  4.n.Mx/2b=0,.can  be  seen  thatthe  upper  Shear  layer  for^sonthe 

penary  s.de  of  the  no.le  tra.Un.edge  contrary  to  expectation.   ^.  can  be  expiated 
■    the  Schlieren  photographs  that  show  an  expansion  wave  extending 
bv  the  observation  in  the  ScWieren  pnoiugi  p 

f  ih.  nn/yle  liD    An  oblique  shock  forms  at  the 
from  inside  the  nozzle  to  downstream  of  the  nozzle  lip.   An        q 

.,.    Streamline  flow  close  to  the  wall  passing  the  upper 
lip  and  crosses  the  expansion  wave.  Streamline 


are 
primary 
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lip  is  turned  towards  the  oblique  shock  emanating  from  the  lip.  Following  this,  the 
streamline  is  turned  further  in  the  same  direction  when  it  crosses  an  oblique  shock  formed 
by  the  reflection  of  the  expansion  wave  from  the  free  boundary.  This  results  in  the 
observed  shear-layer  position.  Also,  from  Fig.  4.13,  it  can  be  inferred  that  -20%  more 
air  flows  through  the  upper  secondary-stream  inlet  than  the  lower  one. 

Figure  4.13  also  shows  that  the  shear  layer  does  not  grow  until  after  the  second 
traverse  at  x/2b=1.58,  which  is  twice  as  long  as  that  determined  from  the  visual 
observations.  It  can  be  seen  that  the  primary/shear-layer  boundaries  narrow  until 
x/2b«9.5,  which  indicates  the  end  of  the  primary  potential-core.  The  secondary/shear- 
layer  boundaries  indicate  that  the  shear-layer  reaches  the  top  and  bottom  walls  around 
x/2b=20. 

Figure  4.14  shows  the  thickness  of  the  shear  layer  as  a  function  of  streamwise 
position.  A  linear  least  squares  fit  of  the  primary  potential-core  data,  excluding  the  data 
at  x/2b=0,  results  in  shear-layer  growth-rates  of 


xJupper) 
8'^,{upper)=0.0Sl,       -ii^  =  -0.19 


(4.8) 


xJlower) 
8'   (lower)  =0.079,       — =  -0.68 


(4.9) 


where. 


^5^-y^^sV^o  <"■"" 


It  can  be  seen  that  there  is  excellent  agreement  between  the  upper  and  lower  growth  rates, 
and  that  the  imperfection  of  the  nozzle  flow  has  no  measurable  effect  upon  the  upper 
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shear-layer  growth  rate.   The  average  growth-rate  of  the  shear  layers  can  thus  be  quoted 


as 


S' ,=0.080  (411) 


These  rates  can  be  compared  to  the  experimental  and  theoretical  work  of  Papamoschou 
(1986)  and  Dimotakis  (1986),  respectively.  The  value  for  the  growth  rate  for  the 
incompressible  free  shear-layer  is  given  by  Eq.  (B.29),  and  using  the  value  of  e=0.14 
suggested  by  Papamoschou  gives  the  value  S'^^  =0.165.  A  curve  fit  to  Papamoschou's 
experimental  impact-pressure  data  (Benjamin,  1990)  gives  the  ratio  of  the  compressible 
shear-layer  growth-rate  compared  to  the  equivalent  incompressible  shear-layer  growth-rate 
with  the  same  velocity  and  density  ratios,  and  specific  heat  ratios  as 

2 

flM  )=0.2+0.8e"^^'^^ 

(4.12) 

Figure  4.15  shows  a  plot  of  this  normalized  growth-rate  curve,  the  present  data  point,  and 
a  compilation  of  other  researchers'  data.  For  the  current  case  flM  \  =0.454  compared  to 
the  present  experimental  value  of  0.48. 

Figure  4.15  shows  that  the  experimental  data  generally  lies  on  or  below  flM  \. 
From  further  analysis  of  the  experiments  of  Benjamin  et  al.  (1993),  it  is  observed  that  the 
data  lies  on  the  curve  when  LjsLj,  and  below  the  curve  when  Li>L2  or  when  the  primary 
potential-core  extends  beyond  the  end  of  the  mixing  section.  Hence,  it  can  be  concluded 
with  good  certainty  that  compressible  shear-layers  in  two-dimensional,  constant-area 
ejectors  have  the  same  growth  rates  as  similarly-conditioned,  free  shear-layers  for 
x-Xgf<Li,  and  Li<L2. 
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4.5.4   Entrainment 

From  these  observations  on  growth  rate,  it  becomes  apparent  that  the  entrainment 
behavior  of  the  semi-contained  and  ejector  shear-layers  can  be  affected  by  containment, 
and  that  the  effect  can  be  a  reduction  in  growth-rate.  Dimotakis  (1986)  proposed  that  the 
entrainment  ratio  is  a  function  of  the  velocity  and  density  ratios  for  the  incompressible 
shear-layer,  specifically. 


^v/''^)  = 


1.0.68(1^1  (4.13) 


Em,oM- ^ 


l+0.68fj^l  (4.14) 


For  increasing  velocity  ratio  and  constant  density  ratio,  the  local  entrainment  ratios 
increase.  However,  for  increasing  density  ratio,  the  local  mass-entrainment  ratio  increases 
while  the  local  volumetric-entrainment  decreases. 

For  the  ejector  shear-layer  where  Li>L2  and  downstream  of  L2+Xgp  the  boundary 
condition  for  the  low-speed  side  of  the  shear-layer  becomes  the  freestream  condition  at 
the  edge  of  the  wall  boundary-layer.  This  boundary  condition  is  also  a  function  of 
streamwise  position.  Hence,  the  effective  secondary-density  and  -velocity  both  increase, 
and  results  in  locally  increasing  values  of  the  velocity  and  density  ratios. 

For  the  current  experimental  case  where  L2>Li,  and  in  the  region  downstream  of 
Li+Xgp  the  boundary  condition  for  what  was  the  high-speed  side  of  the  shear-layer 
disappears.  At  the  centerline,  there  is  no  freestream  fluid  to  be  inducted  into  the  shear- 
layer,  so  that  the  high-  and  low-speed  boundary  conditions  now  become  equal  to  the 
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upper  and  lower  secondary  potential  flow  conditions,  respectively.  In  an  ideal  flow,  these 
conditions  would  be  identical  and  result  in  unity  values  of  both  the  mass  and  volumetric 
local-entrainment  ratios.  However,  Fig.  4.9(b)  shows  that  the  profile  at  x/2b=12.6  is  not 
symmetric,  and  becomes  less  symmetric  by  x/2b=15.8.  This  leads  to  the  conclusion  that 
the  local  entrainment  process  is  asymmetric,  and  it  is  this  asymmetry  combined  with  the 
requirement  that  the  overall  entrainment  ratios  must  be  satisfied,  that  provide  the 
conditions  to  generate  the  large  vorticity  necessary  to  produce  the  unrepeatable  traverses 
at  x/2b=15.8  and  shown  in  Fig.  4.10(a). 

The  compressible,  fi-ee  shear-layer  entrainment  ratios  can  be  predicted  by 
multiplying  Eqs.  4.13  and  4.14  by  the  normalized  growth  rate  (Dimotakis,  1989).  This 
leads  to  predicted  values  for  the  shear-layer  entrainment  in  region  1  of  £^=0.42  and 
E„=0.26.  The  overall  entrainment  ratios  given  in  Table  4.2  are  £^=0.64  and  E„=0.40,  and 
are  about  50%  greater  than  the  calculated  values.  Figure  4.13  clearly  shows  a  distinct 
increase  in  the  angle  of  the  secondary  potential  flow/shear-layer  boundaries  at  x/2b«9,  in 
agreement  with  the  above  scenario. 

4.6   Contour  Plots  and  Visualizations 
Contour  plots  and  color  visualizations  of  the  flowfields  were  generated  from  the 
impact-pressure  traverse  data.   Two  different  software  packages  were  used. 

1)  FloGRAF  is  a  PC-based  graphics  software  package  developed  in  the  mechanical 
engineering  department  at  the  University  of  Florida.  For  contour  plots,  the  code  finds 
a  contour  level  by  searching  inward  from  the  boundaries  using  a  linear  interpolation 
scheme.  Once  a  contour  value  is  located,  it  is  followed  until  no  adjacent  cell  contains 
the  desired  value. 
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2)  FAST  (NASA  Ames  Research  Center.  1992)  is  a  comprehens.ve  workstation-based 
flow  analysts  software  toolktt  that  was  run  on  a  stngle-processor,  Stl.con  Graphics  Iris- 
4D  Crnnson  computer.  The  color  visualizatton  capabihty  of  FAST  is  used  to  provde 
enhanced  insight  into  the  current  experimental  results. 

To  valtdate  that  the  vtsuafeattons  generated  by  the  two  software  packages  were  accurate, 
contour  plots  generated  with  FAST  and  FloGRAF  were  compared  and  found  to  be 
identical. 

4  6  1    Contour  Plots 

contour  plots  of  the  experimental  flowfeld  developed  from  the  impact-pressure 
traverses  are  shown  in  E.g.  4.16.  All  these  plots,  except  for  the  stattc  wall-pressure,  are 
calculated  from  1575  data  potnts  on  a  grid  of  75  potnts  in  the  y-direct,on  by  21  potnts 
in  the  x-dtrection.  The  phystcal  dtmenstons  are  0.635  mm  by  20  mm.  The  stafc  wall- 
pressure  contours  are  calculated  from  132  data  potnts  on  a  gnd  of  12  points  ,n  the  y- 
direction  by  11  points  in  the  x-d,rect,on.  and  the  phystcal  dtmenstons  are  3.4  mm  by 

43.6  mm. 

TWO  regtons  can  easily  be  observed  in  all  the  plots  except  for  the  static  wall- 

pressure; 

1)  the  primary  potential-core  extending  to  x/2b-10,  and 

2)  secondary  potential-flow  extending  to  x/2b-l  5  on  the  lower  side.  On  the  upper 
side,  unmtxed  flow  may  be  traveling  along  the  wall  all  the  way  to  x/2b.30.  even 
though  compression  is  taking  place. 
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The  static  wall-pressure  contour  plot  shows  that  the  static  pressure  is  nearly  constant 
everywhere  in  x-z  planes  except  in  the  region  2<x/2b<5  within  the  primary  potential-core, 
and  in  the  shear-layers  for  0<x/2b<l. 

Increased  detail  of  the  dynamic  pressure  is  shown  in  Fig.  4.17.  To  get  the 
increased  detail,  the  contours  for  these  plots  were  calculated  using  five  times  the  density 
of  data  in  the  y-direction  compared  to  the  plots  in  Fig.  4.16. 

Figure  4.17(a)  clearly  shows  the  path  of  the  momentum  deficit  manifesting  from 
the  boundary  layers  at  the  trailing  edges  of  the  nozzle.  TTie  lower  momentum-deficit 
travels  nearly  horizontally,  while  the  upper  one  moves  up  at  an  angle.  The  reason  for  the 
upper  momentum-deficit's  path  is  thought  to  be  in  comiection  with  the  imperfect  nozzle 
flow  discussed  in  Section  4.5.3.  Both  momentum-deficits  are  convected  to  the  end  of  the 
primary  potential-core,  and  then  diffuse  into  the  surrounding  flow. 

At  x/2b«4.5  and  x/2b«9.5,  distinct,  symmetric,  "bulges"  can  be  seen  on  the 
secondary  side  of  the  shear-layers.  Also,  at  x/2b«6,  symmetric  bulges  can  be  seen  on  the 
primary  side  of  the  shear-layers.  These  are  believed  to  be  indicative  of  high,  local 
entrainment  in  agreement  with  the  large  elliptical  structures  seen  in  the  Schlieren 

visualizations. 

On  the  centeriine  at  x/2b«10,  a  steep  negative-gradient  of  dynamic  pressure  is 
observed  in  the  x-direction,  the  resuU  of  sudden  entrainment  around  the  end  of  the 
primary  potential-core.  In  Fig.  4.17(b),  a  large  vortical  structure  is  visible  centered  at 
x/2b«12.5  on  the  centeriine.  These  observations  suggest  that  there  are  very  high  local 
shear  stresses,  that  are  necessary  to  satisfy  the  entrainment  ratio  criteria  discussed  in  the 
previous  section.  Outside  of  this  vortical  structure,  lines  of  constant  dynamic  pressure  are 
strained  so  that  "fmgers"  appear  in  the  field.  Between  the  fingers,  small  vortices  develop 
that  appear  as  the  shear-layer  reaches  the  wall.  The  most  upstream  vortices  are  located 
at  (x/2b,  y/2b)=(14.3,  1.2),  and  at  the  same  x-  location  adjacent  to  the  upper  and  lower 
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walls.  This  is  the  result  of  secondary  potential-fluid  being  inducted  into  the  shearing 
fluid,  and  remaining  largely  irrotational. 

Figure  4.17(c)  shows  that  at  y/2b«25,  another  line  of  eddies  extends  from  the 
upper  and  lower  walls  indicating  further  entrainment  of  secondary  fluid.  This  suggests 
that  secondary  fluid  has  "slipped"  along  the  walls  for  about  10  nondimensional  units. 
While  this  fluid  is  slipping,  the  fluid  in  the  shear  layer  undergoes  diastrophy  (irrotational 
fluid  is  "corrupted"  with  vorticity  through  the  action  of  viscosity  as  it  cascades  to  the 
Kohnogorov  scale)  which  causes  its  strain  field  to  relax.  When  the  strain  in  the  shear 
layer  is  suitably  relaxed,  the  shear  layer  can  entrain  the  fluid  slipping  down  the  wall,  as 
is  seen  at  x/2b«25.  As  the  fingers  or  eddies  do  not  appear  ftirther  downstream  of  this 
location,  it  can  be  concluded  that  the  fluid  mainly  undergoes  diastrophy  for  x/2b>25. 

4.6.2   Color  Visualizations 

Figure  4.18  shows  color  visualizations  of  the  impact-pressure  traverses.  These 
plots  are  calculated  from  9275  data  points  on  a  grid  of  371  points  in  the  y-direction  by 
25  points  in  the  x-direction.  The  physical  grid  size  is  the  same  as  in  Fig.  4.16,  and  the 
contour  levels  can  thus  be  read  from  that  figure. 

It  can  be  seen  that  the  length  and  shape  of  the  primary  core  is  clearly  identifiable 
and  similar  for  all  of  the  plots.  For  10<x/2b<16,  the  velocity  field  diffuses  most  rapidly, 
followed  by  Mach  number  and  density.  The  stagnation  and  dynamic  pressure  fields 
diffuse  slowest  in  this  region.  In  the  region  15<x/2b<22,  the  effect  of  the  unrepeatable 
traverses  show  up  as  the  zigzagging  of  the  contour  levels.  Hence,  the  contours  in  this 
region  are  representative  of  the  flow,  but  are  not  true  mean  contours.  Further  downstream 
of  this  (x/2bs:20),  the  flowfield  parameters  diffuse  fairly  steadily.  Also,  the  density  field 
differs  from  the  other  fields  in  this  region  as  the  density  increases  instead  of  decreases. 


Table  4. 1   Experimental  flow  conditions. 
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Parameter 

Primary 

Secondary 

Gas 

Air 

Air 

M, 

1.78 

0.30 

A.{mm^) 

321±1 

965±2 

Y/ 

1.4 

1.4 

^0.,  (^^«) 

576  (±0.5%) 

110  (±0.5%) 

P.  (kPa) 

103.4  (±0.5%) 

103.4  (±0.5%) 

ToAK) 

298  (±0.75%) 

298  (±0.75%) 

T.{K) 

182 

293 

m.  (kg/s) 

0.31 

0.12 

V,(mVs) 

0.155 

0.099 

U.(m/s) 

482 

103 

9,,ihlm') 

6.77 

1.29 

p.'ikg/m') 

1.98 

1.23 

QiikPa) 

229 

6.52 

c.  {mis) 

271 

343 

^c, 

0.62 

0.62 

}ji.  iixNs/m  ^) 

12.6 

18.1 

Re'.^ifnm-^) 

75650 

6990 

103 


Table  4.2   Calculated  nondimensional  and  flow  parameters. 


Parameter 

Value 

Em 

0.40 

E. 

0.64 

\r 

0.214 

\s 

0.623 

U^{mls) 

315 

U^  {mis) 

292 

AU{m/s) 

379 

Pa^h/m') 

1.60 

Ca^  i^nls) 

307 

K 

0.62 

M 

1.24 

^^a^^^lm'^) 

15.35 

Re'^imm-') 

41320 

ReJxlO') 

0.957 

Re'^^irnm-') 

53630 

K 

0.165 

V 

^'p^ 

0.075 

K. 

0.091 

Table  4.3   Laminar-to-Turbulent-Transition  Parameters 


^              Re.              Re^ 

9j                       "2 

Re, 

1200              650             7.5x10' 

280           260 

1.6x10' 
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Figure  4.6        Typical  manometer  reading.   Three  reference  pressures  (one  for  each 
mercury  reservoir)  are  indicated  by  tubes  1,  11,  and  21  from  the  left. 
From  left  to  right,  the  rest  of  the  tubes  indicate  increasing  pressure  in 
the  downstream  direction. 
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Figure  4.7       Mean  static-pressure  measurements  on  the  wall.  95%  confidence  intervals 
are  shown. 

(a)  Comparison  of  the  three  types  of  static-pressure  measurements  made 
on  the  left  wall  of  the  mixing  section;  (b)  Difference  between  the  strain- 
guage  and  manometer  static-pressure  readings  on  the  wall. 
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Figure  4.8       Comparison  of  manometer  static-pressure  measurements  on  the  left  and 
right  walls  of  the  mixing  section. 
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Figure  4.9       Impact-pressure  traverses. 

(a)  0^x/2b57.91;   (b)  12.6:sx/2b:fi20.2. 
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Figure  4.10      Unrepeatable  impact-pressure  traverses, 
(a)  x/2b=15.8;   (b)  x/2b=20.6. 
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Figure  4.11      Nonnalized  centerline  values  of  velocity,  density,  and  dynamic  pressure. 
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Figure  4.12     Normalized  impact-pressure  plotted  against  normalized  shear-layer, 
thickness.   Error  function  shown  by  solid  line. 
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Figure  4.13      Shear-layer  edges  determined  from  impact-pressure  measurements. 
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Figure  4.14      Measured  shear-layer  thicknesses,  growth  rates,  and  virtual  origins. 
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Figure  4.15      Comparison  of  normalized  growth-rate  for  ejector  shear-layers  (current 
experiment,  and  Benjamin  et  al.  (1993))  and  free-shear  layers. 
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Figure  4.16      Contour  plots  generated  from  impact-pressure  measurements, 
(a)  Ppi/Po^;  (b)  P0/P02;  (c)  p/p^;  (d)  U/U^;  (e)  M. 
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Figure  4.16--continued.  (f)  Q/Qij ;  (g)  P/Pi,2 
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CHAPTER  5 
EXPERIMENTAL  RESULTS— TIME-ACCURATE  MEASUREMENTS 


5.1  Experimental  Objectives 
The  chief  objective  of  the  experimental  investigation  was  the  measurement  of 
time-accurate  wall-pressure  of  the  potential  and  shearing  flows  in  a  two-dimensional, 
single-nozzle,  supersonic  ejector  discharging  into  a  constant-area  mixing  tube.  Using 
these  measurements,  a  nonintrusive  technique  for  shear-layer  detection  was  developed. 
An  experimental  test-matrix  of  one-  and  two-point  time-history,  wall-pressure 
measurements  were  performed,  allowing  calculation  of  the  following  time  and  frequency 
domain  results: 

root-mean-square  pressure  (second  moment),  Pn^ 


1 


^    7-^cc    1  •'o 


skewness  of  zeroth-derivative  of  pressure  (third  moment),  S 

1  jjm-put 


5  =  limlA_hllA_  (5.2) 

Pirns 


T-*«T  n3 


kurtosis  of  zeroth-derivative  of  pressure  (fourth  moment),  K 

rrms 
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pressure  correlation,  Rpp 
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(.,    ,./    ^  44'^'^!^ 


R. 


pp 


\2b'  2b'  ) 


2b   2b 


(       I        I 
x+x    y+y 


[2b  '   2b 


,/+T 


'  X     y 
,2b'  2b  '^"^ 


x+x    y+y 


2b      2b  ) 


(5.4) 


pressure  spectral-density,  S 


pp 


pp 


[2b'  2b 


,G) 


277'' -» 


2b'  2b 


,T 


dr 


(5.5) 


transfer  function  of  pressure,  Hy 


^ 


[2b   2b     ) 


M2Z)   2Z? 


i2Z)    2Z?      J 


ij  =21, 12 


(5.6) 


where 


-Spp{0,0,<^) 


is  the  autocorrelation  at  (x",  y"),  and 


(5.7) 


x"=x,  y"=y  for  ij=21 
x"=x+x',  y"=y+y'  for  ij=12. 
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coherence  of  pressure 


y 


The  d.screte  calculation  procedure  for  the  above  e.uat.ons  was  outUned  .n  Chapter  3. 

SX^econd^Thiri^nlZm^^ 

^  9  1    Mixing- g---^'""  T"'gt  Measurements 

Extra  care  was  taken  to  ensure  that  the  ,n,e.  flow  measurements  were  as  cons.sten. 

.,.,      c     .o.h  „f  the  twelve  pressure-transducer  positions,  at  least 
and  accurate  as  possible.   For  each  of  the  twelve  p 

.       „f     -in  frame      4096-point      windows      (163,84U 
two     separate     measurements     of     40-frame, 

points/measurement)  were  made.  All  of  these  cases  were  measured  usmg  a  100  kHz 
bandwidth  with  a  resulting  data  collection  time  of  0.64  seconds.  The  rms  press.es 
between  separate  runs  of  the  same  case  differed  by  less  than  3%  for  all  12  positions.  For 
oonfideneeintervalealculations,thedatawasspl,tmto,0wi„dowsofl6,3S4points^The 

..oiceoftenwindowswasarrivedatarbit^lybycompanngseveralcaseswithdifferent 

confidence  interval  was  found  to  be  less  than  .5%  in  all  cases. 

Figure  5.1  shows  p.  normalized  by  local  values  of  P.,  and  Q.  as  we,,  as  plotted 
in  decibels.   Tlie  95%  confidence  intervals  are  plotted  for  p^.-   K  can  be  seen  that 
*ere  is  ,itt,e  difference  in  the  shapes  of  the  three  cuives  except  for  p^Q  a. 
:^/2b.-0  3.    This  may  be  caused  by  the  measurement  of  the  impact  pressure  at  a 
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different  time  from  p_  or  could  also  be  an  interference  effect  of  the  impact  pressure 

probe. 

Values  of  p^Q  have  been  quoted  in  the  literature  by  several  researchers  and 
there  is  much  discrepancy  between  the  cited  values.  Mull  and  Algranti  (1956)  found  that 
PrJQ  has  a  value  of  0.003  at  M=0.2,  and  linearly  decreases  to  a  value  of  about  0.0012 
for  M>0.55.    Measurements  by  Schewe  (1983)  in  incompressible  flow  recorded  values 
of  PrJQ  using  different  diameter  transducers.    Using  transducers  with  nondimensional 
diameters  of  19^d^=duyv^333  (1  mm^d^lS  mm),  Schewe  found  that  p^Q  had  a  value 
of  0.010  for  d^=19  and  decreased  almost  linearly  to  half  this  value  for  d^=333. 
Karangelen  et  al.  (1993)  confirmed  these  values  for  33^d^=duyv^66.  Burton  (1973)  made 
incompressible  flow  measurements  and  found  values  of  p^Q  for  smooth  plate,  favorable 
gradient  flow  of  0.0100,  and  rough  plate,  favorable  gradient  flow  of  0.0174.  Burton  also 
measured  adverse  pressure-gradient  flows  with  a  value  of  p^QM).008  for  a  smooth  plate 
and  p^Q«0.014  for  a  rough  plate.  Eckelmann  (1989)  reports  that  the  dependence  of  the 
ratio  p^Q  on  d*  is  not  unique.    The  reason  for  this  is  that  there  is  a  difference  in 
Reynolds  number  that  necessarily  occurs  when  two  sensors  of  different  diameter  have 
equal  values  of  d\   Kistler  and  Chen  (1963)  measured  p^Q-^.005  under  a  supersonic 

boundaiy  layer  for  Mach  numbers  in  the  range  1.33-5.0;  however,  these  measurements 

were  made  at  very  low  Reynolds  number. 

Pressure  measurement  values  for  the  current  research  have  the  following  ranges: 

O.OOl^p^Q^O.2,  0.002^p,^„s0.013,  and  140.5^dB_j)™,^156.4.   The  values  of  p^Q 

in  the  secondary  are  much  higher  than  any  other  reported  values  due  to  the  acoustic 

contribution  from  the  supersonic  primary. 

The  skewness  of  the  inlet  measurements  are  plotted  in  Fig.  5.2.  It  can  be  seen  that 

for  y/2b>0.6,  the  skewness  rapidly  becomes  negative  in  contrast  to  that  for  y/2b<-0.6 

which  is  almost  constant  with  a  small  peak  towards  the  center  of  the  secondary  flow  with 
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a  value  of  0.25.  The  reason  for  this  rapid  change  is  believed  to  be  caused  by  the 
imperfect  flow  from  the  positive  trailing-edge  of  the  nozzle,  which  was  discussed  in 
Chapter  4.  The  skewness  in  the  primary  flow  has  a  value  close  to  zero.  It  can  be  seen 
that  there  is  a  distinct  demarcation  between  the  primary  and  secondary  flow. 

The  inlet  kurtosis  is  shown  in  Fig.  5.3.  On  the  positive  side  of  the  mixing  section, 
the  kurtosis  rises  rapidly  for  y/2b>0.3  from  a  value  of  3.4  to  5.1.  Again,  this  is  a  result 
of  the  imperfect  nozzle  flow.  It  can  be  summarized  that  the  general  trend  of  the  kurtosis 
is  that  it  is  nearly  constant  with  a  value  of -3.5,  and  decreases  to  ~3.0  at  the  outer  edge 
of  the  secondary  flow. 

Further  insight  can  be  gained  from  the  probability  density  functions  (pdf  s)  of  the 
pressure  measurements,  and  are  shown  in  Fig.  5.4.  It  can  be  seen  that  the  pdfs  all  lie 
close  to  the  Gaussian  distribution,  however  the  pdfs  for  y/2b^0.81  lie  furthest  from  the 
Gaussian  curve  and  the  other  data  points.  This  again  attests  to  the  different  flow 
characteristics  present  in  this  region. 

The  pdfs  can  be  more  closely  inspected  by  calculating  the  probability  W  that  the 
wall  pressure  amplitude  p  exceeds  a  threshold  value  p^.  These  calculations  are 
summarized  in  Table  5.1  and  compared  to  the  results  of  Schewe  (1983). 

Table  5.1  shows  that  W(p<-pJ  and  W(p>pJ  for  y/2b=0.81  differ  significantly 
from  the  values  of  W(p<-pth)  and  W(p>P(h)  at  the  other  values  of  y/2b.  However,  the 
values  of  W(|p|>p^  in  the  three  current  cases  are  almost  equal.  The  major  differences 
for  y/2b=0.81  occur  for  pVPnns<l,  P«/Pnns<2,  and  pu/p„„3>l.  Pti/Pnns>2,  the  effect  being  that 
there  are  more  negative  pressure  events  than  positive  ones  aty/2b=0.81,  and  more  positive 
pressure  events  than  negative  ones  at  y/2b=-1.87  and  y/2b=0. 

In  comparison  to  Schewe's  measurements,  it  can  be  seen  that  there  is  good 
agreement  up  to  the  threshold  Pa/Pnns=-1'  and  values  less  than  Schewe's  for  Pfl/Prm3>2  and 
Pu/Pmis*^~2.    Also,  there  is  an  increased  number  of  pressure  events  for  the  current  flow 
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u         ™    2<o/D    <-l     However,  all  cases  in  the  region 
compared  to  Schewe's  m  the  reg,o„ -2<p/p.  „„  „,,e  are  fewer  events 

l<|p,/p.<2  have  probaHHty  less  *a„  the  Gauss.an  value.   Also  there 

.jLe— .notaeteetthes-estsoaleso...^^^^^^^^^^ 

,ereas,ng  .ansducer  s.e.  H  w^  _^^^  ^^^  ^^^_^^^^^  ^  ^^  ^^^^^^^^ 

and  negative  pressure  events  for  ,pj  ^^^^  ^  ^^^,^  ^^  ^^^  „^,^^ 

are  truly  different  from  the  incompressible  case. 

on  e  ntegral  length  of  the  ,„ner  wall-layer,  and  hence  correspond  to  .ncrease    rnner- 
of  the  integral  ic  5  Simultaneous 

„„.er  layer  .nteracfon  relafve  to  .nco^press.ble  boundary-layer  fl 

confirm  this  scenario. 

<;  9  9    rrntrrlinf  Mpasinrements 

^^^-^^;;::::;;:;;::„u„r  of  data  and  co.pu.Uon  Ume  re..ed  to  store  . 
^„,,,,  .,  ,,.  „™e„t  dat.  the  rest  of  the  ^easure^en.  and  results  preserve 
Z  are  based  on  10  fra.es  of  .0,.  ..nts  .0,.0  „easurenren,..s,ho„)  w.. 
l,uar.erthenu.berta.enforthe.nletnreasure.e„.«snu».rwasfoundtog.ve 

■tv,  .n  averaee  95%  confidence  level  for  p^  ot  ±5.i/o. 
adequate  accuracy  w.th  an  average  9  ^^^  ^^^^^  ^^  ^  ^_^^^^__ 

A  plot  of  p^.  is  shown  m  F.g.  5.5.    It         be 
.ncrease  ,    p     .st  downstream  of  the  .nlet.  T..S  ,s  caused  by  acoustrc  radrafon  f.n, 
rn    rellt  ng  fro.  the  boundary-layers  on  the  no.,e  walls,  and  fron,  rad.afon 
the  wake  g  3hear-layers.  From  the  inlet  downstream 

emanating  from  the  mitial  development  of  the  1 J 

.0  ..b-n  there  are  a  senes  of  pea.  and  troughs,  whrch    p      .^I 
measurements,  were  rdentifed  as  correspondmg  to  waves  .n  the  flow. 
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I  .e  pel  a„.  .ou^s  .0  „c.  —  .  wouM  .  .pe.ea  ..  p.e  waves  .n  a 

„fx/2^20*ep     co„,.nues.o,ncreaseun.,U.pea^a.x/2b-26v,he.the 
Downstream  or  xi  10=10,  mc  Pn^s  v 

u       f  1  r.ee  Fk  4  16)  and  then  starts  to  decrease  more  rapidly 
flow  has  a  Mach  number  of  1  (see  Fig.  4. 10;,  a 

than  it  approached  the  peak.  i.th.tthe 

Pigure  5.6  shows  *e  centerHne  skewness.  The  im,ned,a.e  .mpress.on  .s  .hat 

t     ,;-!k<l8     Lame  confidence  intervals  are 
confidence  intervals  are  much  greater  for  x/2b<18.    Large 

.   .   „^5  7  10  and  17   The  las.  two  are  known  to  be  the  locat.ons  of  the  end 
observed  at  x/2b«5,  /,  lu,  anu  1/.    ^  ..,11 

e  and  the  beginn.ng  of  the  in.eraction  of  the  shear-layer  w,.h  the  wall, 
of  the  primary  core  and  the  beginning 

.    /0K_n  asain  supports  the  premise  that  large-scale 
respectively.    The  observation  a.  x/2b-17,  again,  supp 

„.i„g  IS  presen.  ,n  this  region.   The  reason  for  the  observations  at  x;2b=5  and  7  a 
Ipeld  tl  associated  with  noise  sources  of  the  pr^ary.  Xam  0.2)  Shows  far-field 

measurement  results  of  p„  for  an  axisymmetnc  Mach  1.3  .et  that  has  apparent  noise 

,.  x/D-6  and  x«.20.    He  also  notes  that  the  frequency  of  the  noise  sources 
sources  at  x/D«o  ana  x/i--  ^"  „  •   _a 

,    ,  ^n  VHz  as  x/D  increases.  The  value  of  the  skewness  is  -0 

decreases  significantly  from -50  kHz  as  x/i;  inc 

for  x/2b<20  and  then  decreases  with  increasing  rapidity  to  a  value  of  -.6  at  x/2  . 
P,gure  5.7  shows  the  centerline  kurtos.s.  Urge  values  of  confidence  intervals  are 
o^crved  at  the  inle,  .2b»5,  7,  10  and  17,  the  same  locations  at  which  such  observa.o. 
«cre  made  for  the  skewness.  From  the  obse^ation  that  both  the  skewness  and  kurtosi 
.hibit  Wide  confidence  intervals  at  these  locations  compared  to  other  measuremen 
iccafions  it  can  be  concluded  that  the  measurements  are  either  being  directly  influenced 

■      or  that  unsteady  acoustic  sources  cause  this  effect.  As  it  will  be  seen 
by  the  outer  region,  or  that  unsteaoya^. 

Iw  the  spectral  density  is  ve.  repeatab.e  at  all  frequencies,  and  this  most  probably 
„,es:utunsteadyacousticsources..fitisassumeatha.theeffectivetransducerdiameter 

,  larger  than  *e  thickness  of  the  inner  layer  (d->33),  then  the  wide  confidence  inte^als 


129 

for  the  skewness  and  kurtosis  must  originate  from  the  outer  layer  (0.2<y/6<l),  and 
therefore  from  the  main  flow.  Once  again,  this  suggests  the  presence  of  vigorous  large- 
scale  mixing  at  or  near  these  locations. 

5.2.3    Measurements  at  v/2b=-1.64 

Figure  5.8  shows  Pn^„  plotted  against  x/2b  for  y/2b=-1.64.  There  is  a  small 
local  peak  at  x/2b«2  which  corresponds  to  the  shear-layer  transition  to  turbulence.  The 
peak  value  is  -0.05  and,  like  the  centerline  peak-value  position,  is  located  at  x/2b«26 
which  is  the  streamwise  location  where  the  centerline  Mach  number  becomes  subsonic 
(see  Fig.  4.16).  The  skewness  and  kurtosis  are  shown  in  Figs.  5.9  and  5.10,  respectively. 

Comparisons  of  Pj^„,  S  and  K  with  the  centerline  values  are  shown  in  Figs. 
5.11,  5.12  and  5.13,  respectively.  The  turbulent  pressures  have  a  similar  overall  trend, 
but  differ  in  that  the  centerline  has  small  peaks  and  troughs  up  to  x/2b«18.  The 
skewnesses  also  show  a  similar  overall  trend,  except  that  the  centerline  values  have  peaks 
and  troughs  all  the  way  to  the  end  of  the  mixing  section.  For  x/2bsl7,  the  skewness  at 
y/2b=-1.64  is  always  negative  while  the  centerline  skewness  has  both  positive  and 
negative  values.   The  comparison  of  kurtosis  shows  a  similar  overall  trend. 

These  results  show  that  the  streamwise  locations  of  shear-layer  transition  to 
turbulence,  and  the  end  of  the  supersonic  core  on  the  centerline  can  be  confirmed  by 
measurements  parallel  to  the  centerline  that  cross  the  secondary  stream. 

5.2!4   Whole-Field  Measurements 

Measurements  were  taken  over  the  entire  wall  of  the  mixing  section.  The 
measurements  were  collected  at  132  locations  on  a  grid  of  12  evenly-spaced  points  in  the 
y-direction  (-1.87:sy/2b^l.08)  by  11  evenly-spaced  points  in  the  x-direction 
(0^x/2bs34.5).  Figure  5.14  shows  a  color  plot  of  p^,  Pno/Pw,  and  p^^  in  acoustic 
decibels  which  are  defined  as 
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dBa  =20\og^ 

P'ef  (5.9) 

where  Pref=20  ^Pa. 

The  plot  of  Pnj^  shows  the  primary  potential-core  in  dark  blue  that  extends  to 
x/21>«9.5  in  agreement  with  the  impact-pressure  measurements.  From  x/2b«6  just  outside 
the  primary  core,  high  values  of  ^^^  are  seen  to  begin.  This  region,  shown  in  green, 
grows  to  the  centerline  and  the  outer  wall  by  x/2b«15.  It  is  known  from  the  mean 
measurements  that  this  is  the  location  where  the  shear  layer  reaches  the  wall.  It  can  also 
be  seen  that  the  highest  p,^  gradients  are  located  near  to  these  locations,  where  it  has 
already  been  observed  that  the  skewness  and  kurtosis  along  the  centerline  have  large 
confidence  intervals,  as  well  as  being  noise-source  locations. 

The  plot  of  p„^„in  Fig.  5.14(b)  shows  similar  information  to  that  of  Pn„g  in  Fig. 
5.14(a).  Hence,  either  of  these  plots  can  be  used  to  identify  the  locations  of  the  ejector 
regions.   The  plot  of  p,^  in  dBa  emphasizes  the  region  up  to  x/2b«10. 

Figure  5.15  shows  color  plots  of  the  skewness  and  kurtosis.  These  show  clear 
peak  values  on  the  centerline  at  x/2b"7  of  S= 1.85  and  K=l  1.5.  These  values  indicate  that 
there  is  a  high  probability  of  large- valued  pressure  events,  and  that  the  frequency  of  large- 
valued  negative  events  is  less  than  that  of  positive  ones.  Figure  5.16  shows  the  pdf  for 
the  peak  value,  and  it  is  clear  that  there  are  many  positive  events  far  from  p/pnns=0.  The 
range  of  values  is  -3sp/p^s6.2.  Table  5.2  shows  the  breakdown  of  the  pdf.  From  the 
increased  number  of  events  for  |p|<l  compared  to  the  number  in  the  same  region  at  the 
inlet  (see  Table  5.1),  it  is  concluded  that  the  noise-source  effect  is  apparent  due  to  a 
relative  increase  in  inner-outer  layer  energy-transfer.  The  effect  in  the  inner  layer  is  a 
relative  increase  in  large  positive  pressure  events.  It  is  interesting,  yet  unexpected,  to  note 
the  numerical  agreement  between  the  results  for  W(|p|>3)  and  W(|p|>4)  of  the  current  data 
with  Schewe's,  especially  since  the  skewness  is  of  opposite  sign  in  this  case.   This  again 
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(see  section  5.2.1)  points  to  the  different  behavior  of  the  wall  pressure  in  incompressible 
and  compressible  flow. 

5.3   Autospectra 
Autospectra  were  calculated  from  the  same  data  as  presented  in  Section  5.2,  with 
the  inlet  autospectra  calculated  from  40  frames  of  4096  points,  and  the  rest  from  10 
frames   of  4096  points.      All  autospectra  plots   are  of  spectral   density   O  against 
dimensionless  wavenumber  or  frequency  (Hz),  and 


p. 


2 
rms 


(5.10) 


where  BW=bandwidth  (Hz)  and  NP=number  of  points/frame.  A  bandwidth  of  100  kHz 
was  used  for  all  one-point  measurements.  The  dimensionless  wavenumber  is  plotted  as 
k=(ji8/U,  where  the  shear-layer  characteristic  length  scale  6  is  defined  as  follows: 

(a)  if  the  measurement  is  in  the  secondary  flow  the  nearest  shear-layer 
thickness  is  used, 

(b)  if  the  measurement  is  in  the  shear  layer,  then  that  shear-layer  thickness  is 
used, 

(c)  if  the  measurement  is  in  the  primary  then  the  average  thickness  of  the  two 
shear-layers  is  used,  and 

(d)  for  flow  downstream  of  where  the  shear-layer  reaches  the  outer  wall,  the 
mixing-section  height  is  used. 
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5.3.1    Mixing-Section  Inlet  Measurements 

Autospectra  of  separate  runs  on  either  side  of  the  primary  and  within  the  primary 
are  plotted  in  Fig.  5.17.  It  can  be  seen  that  for  each  pair  of  autospectra,  there  is  excellent 
agreement.  Figure  5.18  shows  results  from  just  outside  of  the  shear-layers  and  in  the 
lower  shear-layer  at  the  inlet.  The  axes  have  been  scaled  so  that  the  slopes  are  at  the  true 
angle.  Measurements  for  Figs.  5.18(a)  and  (d)  are  located  in  the  secondary  flow,  Fig. 
5.18(b)  is  in  the  lower  shear-layer,  and  Fig.  5.18(c)  is  in  the  primary  flow.  In  all  of  the 
plots,  the  sharp  peaks  indicate  the  presence  of  acoustic  emissions  which  will  be 
considered  below. 

All  of  the  plots  show  a  small  k'  range  at  low  wavenumber.  There  is  a  second  k' 

range  in  the  measurement  in  the  shear  layer  (Fig.  5.18(b))  between  the  k'^'^  and  k'"'^ 

ranges.     The  second  k'  range  is  located  at  a  wavenumber  of  just  over  1  order  of 

magnitude  higher  than  the  first  k'  range.   If  it  is  assumed  that  the  two  ranges  originate 

from  the  same  eddy  travelling  between  the  primary  and  secondary  overlap  layers  in  the 

boundary  layer,  then  the  wavenumber  of  these  ranges  should  be  approximately  equal 

when  normalized  by  the  boundary-layers'  freestream  velocities  and  thicknesses.   Hence, 

the  ranges  should  differ  by  the  ratio  of  their  boundary-layer  thicknesses  divided  by  r. 

Using  values  for  6  from  Table  3.1  and  the  value  of  r  from  Table  4.2,  it  is  found  that 

—  • — —13  which  is  close  to  the  observed  difference  in  (o6/U  between  the  two  ranges. 
r  8^ 

Therefore,  it  can  be  concluded  that  the  measured  spectrum  contains  information  from 
eddies  travelling  through  primary  and  secondary  regions  of  flow.  It  also  suggests  that  the 
conclusion  by  Thomas  and  Bull  (1983)  that  the  contributions  from  various  source  regions 
will  fall  off  rapidly  with  their  distance  from  the  point  under  consideration  may  not  always 
be  applicable.  Indeed,  depending  on  the  investigator's  requirements,  the  questions  that 
needs  to  be  addressed  before  designing  a  turbulent-pressure  experiment  is,  "what  feature 
do  you  want  to  measure,  and  from  which  region  of  the  flow?" 
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At  all  the  inlet  locations,  it  is  apparent  that  the  significant  acoustic  emmisions  tend 
to  be  located  between  turbulent  ranges,  a  distinct  example  being  between  the  k'  and  k'^'^ 
regions  in  the  shear  layer  (Fig.  5.18(b)).  Also,  the  measurement  of  the  k'^'^ 
(turbulence-mean-shear)rangeathigherwavenumberthanthe^--(turbulence-turbulence) 

range  is  an  unexpected  result  in  view  of  the  calculations  of  George  et  al.  (1984). 
However,  their  calculations  were  based  on  unbounded  shear-flow,  and  it  is  this  difference 
that  may  be  the  cause  of  this  observation.    If  the  freestream  shear-flow  is  assumed  to 
behave  as  George  et  al.  predict,  that  the  turbulence-turbulence  interactions  dominate 
everywhere  except  near  the  eddy-producing  region,  then  it  can  be  hypothesized  that  the 
measurement  technique  records  the  turbulence-turbulence  and  turbulence-mean-shear 
interactions  from  the  freestream  and  boundary-layer,  respectively.  The  turbulence-mean- 
shear  interaction  generated  in  the  boundary  layer  dominates  the  high  wavenumber  regime 
because  of  its  closer  proximity  to  the  pressure  transducer  than  the  turbulent-turbulent 
interaction.    This  is  supported  by  Kraichnan's  (1956b)  prediction  for  turbulent  wall- 
pressure  under  a  boundary  layer  in  chamiel  flow,  that  the  turbulence-mean-shear 
contribution  to  p_  is  5.6  times  that  of  the  turbulence-turbulence  contribution. 

5  3.2   Centerline  Measurements 

Thirty-six  centeriine  autospectra  are  shown  in  Fig.  5.19  that  range  in  location  from 
x/2b=-0.32  to  x/2b=35.00.  Slopes  for  the  k',  k^'\  and  k'"'  ranges  are  shown  in  the  toi> 
left  chart  of  each  page.  Agam,  the  axes  have  been  scaled  so  that  the  slopes  are  at  the  true 

angle. 

There  are  five  major  characteristics  of  the  spectra  that  are  observed.  The  first  is 
the  strong  presence  of  a  large  acoustic  peak  at  o,6/U~l  from  x/2b=^-0.32  to  x/2b-4.  Its 
spectral  location  is  initially  between  the  k^'^  and  k'"^  ranges.  By  x/2b-4,  the  *"'-  range 
has  disappeared  and  the  diminishing  acoustic  disturbance  lies  on  a  k'^'  range.  It  takes  on 


134 


i.  ™axtau™  spectra,  energy  de„s,ty  value  a.  x/^b-O.  which  sugges.  .hat  ,t  arises  fro. 
U,e  nozde  Up  area.   Th.s  w,ll  be  discussed  further  in  Section  5.5. 

The  second  feamre  is  a  dismrbance  that  occurs  at  «*/U-^.l  at  the  start  of  the  k 
range  and  can  be  first  seen  clearly  at  x/2b=0.60.  Between  x/2b=,  .79  and  x/2b=2.40  there 
i,  an  order  of  magnimde  increase  in  ,.s  spectra,  energy  dens.ty.  wh.ch  remains  constant 
unttl  the  sutroundtng  spectral  energy  dens.ty  reaches  the  same  level  at  x/2b-5.    Tlte 
phys,ca,  locatton  of  th,s  peak  suggests  that  ,t  ,s  assocated  w.th  the  shear  layer  transttton- 

to-turbulence. 

The   third   feature   is   the   appearance   of  a   sharp   slope   at   X/2WD.60   and 
0  0Ka,6AJ<0.1.   The  slope  steepens  and  shrinks  in  wavenumber  range  until  x/2b-3.59, 
after  which  .t  ,s  broken  up  by  a  broad  spectra,  peak  at  x/2b=4.,9.  For  x/2b^.l9.  there 
are  large  d.sturbances  ,n  the  spectral  energy  denstty  of  wavenumbers  at  values  of  th,s 
peak  and  lower.  At  X/2W.96,  >t  can  be  seen  that  two  peaks  appear  and  pers.s.  through 
x/2b-7  57    This  is  the  location  of  the  peak  positive  skewness.  peak  positive  kurtosts,  U,e 
star,  of  the  regton  of  mcreased  entra.nmen.  on  the  pnntary  s.de  of  the  shear-layers,  and 
the  location  of  a  major  noise  source.  The  two  peaks  reappear  at  x/2b=9.96.  .he  end  of  the 
primal  core  and  remain  distmctive  un.,1  k/2^,2.,3.  the  center  of  the  large  vortex  (e.g.. 
see  F,g.  4.17(b)).   It  ts  surmised  that  the  two  peaks  artse  due  to  the  in.eractton  of  the 

Upper  and  lower  shear-layers. 

The  fourth  observafon  is  that  the  lowest  wavenumbers  show  generally  increasmg 
or  decreasing  trends  at  different  centeriine  locattons.  Tlte  most  consistent  trend  ts  seen 
at  x.15.13  where  the  low  wavenumbers  generally  keep  increastng  This  is  *e  reg,on 
downstream  of  where  the  secondary  side  of  the  shear  layers  reach  fte  wall. 

The  fifth  observation  is  ftat  of  the  trend  in  the  spectra,  decay  of  the  pressure.  As 
has  been  a,ready  mentioned,  by  x/2b^  .he  initia„y  present,  sma,,  *'-  range  has 
dtsappeared  leaving  just  a  ,'"  range.  At  x/2b=6.96.  the  location  of  the  center  of  a  notse 
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source,  a  small  k'"^^  range  reappears  and  remains  small  until  x/2b«20,  after  which  it 
becomes  more  obvious  and  is  very  clear  at  x/2b=35.00.  At  x/2b«16,  the  k'^^^  range  starts 
to  shrink  and  a  k''  range  starts  to  develop  in  its  place.  At  x/2b=35.00,  the  k''  range  almost 
dominates  the  whole  wavenumber  range.  These  major  changes  in  the  spectral  behavior 
are  probably  the  result  of  the  switch  from  adverse  to  favorable  p^  gradient  at  x«27  (cf 
Fig.  5.5).   The  mean-pressure  gradient  is  adverse  for  x/2b>15  (cf.  Fig.  4.8). 

5.3.3   Measurements  at  y/2b=- 1.64 

Figure  5.20  shows  the  autospectra  calculated  from  these  measurements.  The 
autospectra  shown  in  the  secondary  flow  in  the  range  -0.32^x/2bsl0.02  show  similar 
general  features.  In  the  region  0.1<{jo6/U<10,  an  approximately  k'^^^  range  is  observed. 
It  is  highly  disturbed  by  acoustic  peaks  radiating  from  the  primary  jet.  The  acoustic 
radiation  will  be  discussed  in  more  detail  in  Section  5.5.  At  x/2b=1.19  and  co6/U=10,  a 
large  acoustic  peak  is  observed  that  sits  atop  a  spectral  peak.  At  (o5/U>10,  the  spectrum 
decreases  at  just  higher  than  a  k'"^^  slope. 

Between  x/2b=10.02  -  16.05,  the  spectra  change  gradually,  but  significantly.  The 
acoustic  disturbances  disappear  and  three  spectral  peaks  are  observed  at  co6/U~0.1, 
(o6/U~5  and  (o6/U~50.  Also,  the  spectral  slopes  change  to  k''  at  low  wavenumber,  k'^^^ 
in  the  middle  range,  and  remains  at  just  over  k'"^^  in  the  high  range.  Multiple  spectral 
peaks  have  been  observed  in  the  spectrum  of  oceanic  turbulence  (Monin  and  Yaglom, 
1987,  p.  526)  and  are  characterized  by  an  appreciable  influx  of  external  energy,  such  that 
in  the  gaps  between  separate  inertial  ranges  the  energy  influx  and  dissipation  play  a  minor 
role.  The  oceanic  spectral  peaks  correspond  to  global  scales,  tidal  scales,  and  wind 
(surface  waves)  scales.  It  is  hypothesized  that  the  current  peaks  can  be  attributed  to  large 
scales  of  the  shear  flow.  Unlike  the  oceanic  energy  influxes,  it  is  proposed  that  the  influx 
of  energy  from  the  shear  layer  at  co6AJ~5  is  large  enough  to  invalidate  the  Kolmogorov 
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hypothesis  for  the  inertial  subrange,  resulting  in  a  decrease  in  the  spectral  slope  decay 
from  k'^^^  to  k'^'^.  The  k'^'^  range  persists  until  x/2b=19.63  and  then  gradually  increases  to 
a  k'  range  by  x/2b=26.56. 

In  Chapter  4,  the  region  15<x/2b<20  where  the  k'^'^is  observed  was  found  to  be 
unstable,  and  it  was  hypothesized  that  an  overall  entrainment  ratio  adjustment  takes  place 
in  this  region  through  increased  vorticity.  Since  vortex  stretching  is  the  only  known 
spectral  energy-transfer  mechanism  (Tennekes  and  Lumley,  1987),  vortex  stretching  in 
this  region  must  be  increased.  This  implies  that  the  spectral  energy  flux  T  and  the 
dissipation  rate  s  increase  (T=e  in  the  inertial  subrange).  Further,  this  implies  that 
diastrophy  is  increased  and  leads  to  the  conclusion  that  large  eddy  strength  also  increases 
in  this  region. 

5.4   Two-Point  Spectral  Measurements 

The  motivation  for  the  two-point  measurements  was  to  identify  separate  areas  of 
the  flow  where  in-phase  or  anti-phase  frequencies  exist.  This  allowed  determination  of 
whether  the  frequency  originated  from  acoustic  radiation  or  from  a  turbulence  interaction, 
and  from  what  physical  feature  of  the  flow. 

All  of  the  two-point  spectral  measurements  were  calculated  from  40  frames  of 
4096  points  (100  kHz  bandwidth  unless  noted)  measured  at  each  transducer  location.  The 
data  is  presented  by  phase  difference,  cross-spectral  density,  transfer  function,  and 
coherence  which  were  obtained  using  Equations  (5.5),  (5.6),  and  (5.8). 

Figure  5.21  shows  the  results  from  a  measurement  on  the  inlet  plane  in  the 
secondary  flow  on  either  side  of  the  primary.  Peaks  1  and  8  are  in  phase  while  peaks  2  - 
7  are  out  of  phase.  The  out  of  phase  peaks  are  typical  of  the  helical  (Kelvin-Helmholtz) 
mode  instability  measured  on  opposite  sides  of  a  supersonic  jet  (Grutmark  et  al.,  1989, 
1990).   Figure  5.22  shows  a  measurement  on  the  inlet  plane  in  the  primary.   The  largest 
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peak  is  at  500  Hz  compared  to  the  highest  peak  in  the  secondary  at  312.5  Hz  (of.  Fig. 
5.21).  All  of  the  peaks  are  in  phase  except  for  the  peak  at  3.625  kHz,  the  highest 
powered  jet-noise  frequency  observed  in  Fig.  5.21. 

The  results  in  Fig.  5.23  were  calculated  from  a  20  kHz  bandwidth  measurement 
along  the  centerline  just  before  the  end  of  the  primary  potential  core.  It  can  be  observed 
that  5  peaks  ranging  from  187.5  Hz  to  587.5  Hz  are  close  to  being  in  phase.  The  peak 
at  487.5  Hz  is  the  most  in-phase  peak  with  (j)=-l°  which  corresponds  to  the  500  Hz  peak 
with  <J)=0°  seen  at  the  primary  inlet  (cf.  Fig.  5.22).  This  result  signifies  that  the 
characteristic  frequency  of  the  primary  core  is  -500  Hz. 

Figures  5.24  and  5.25  show  transverse  measurements  across  and  below  the 
centerline,  respectively,  at  the  streamwise  location  of  the  end  of  the  primary  potential 
core.  At  the  centerline.  Fig.  5.24  shows  that  the  dominant  peak  is  750  Hz,  (|)=12°  with 
a  first  harmonic  at  1.4375  kHz,  <j)=0°.  This  is  the  characteristic  frequency  of  the  flow 
mixing  around  the  end  of  the  primary  core,  and  the  reason  for  the  generation  of  the  large 
eddy  centered  at  x/2b~12.5  (cf  Fig.  4.17).  Figure  5.25  shows  many  of  the  peaks 
observed  in  the  primary  and  secondary.  However,  the  coherence  spectrum  shows  that  the 
687.5  Hz,  <j)=6°  peak  is  the  most  coherent.  The  first  and  second  harmonics  appear  as 
peaks  4  and  5.  Peaks  6,  7,  8  and  9  arise  from  the  acoustic  jet  noise.  The  characteristic 
frequency,  3.625  kHz,  is  seen  to  be  out  of  phase  across  the  centerline  and  in  phase  below 
it,  in  keeping  with  the  helical  mode  instability  mentioned  above. 

Figures  5.26  -  5.29  show  transverse  spectra  at  the  streamwise  location  just 
downstream  of  where  the  shear  layer  reaches  the  wall.  Figures  5.26  -  5.28  show  that  the 
peaks  are  almost  in  phase  for  measurements  up  to  half  the  mixing-section  height.  Also, 
the  dominant  low  frequencies  move  125  Hz  for  the  half  mixing-section  measurement  (cf 
Fig.  5.28).  Figure  5.29  shows  that  the  measurement  over  either  side  of  the  centerline 
indicates  that  the  low  frequencies  move  even  lower  and  the  coherence  increases  as  they 
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go  lower.  Again,  the  peaks  over  3  kHz  are  out  of  phase  and  can  thus  be  attributed  to  jet 
noise,  and  hence  do  not  play  a  part  in  large  scale  mixing.  The  low  frequencies  observed 
at  this  streamwise  location  indicate  that  the  integral  length  scale  of  the  turbulence  is 
increasing  toward  the  characteristic  dimensions  of  the  mixing-tube. 

5.5  Narrowband  Measurements 
From  the  autospectra  and  cross-spectra,  it  can  be  observed  that  there  are  three 
distinct  frequencies  that  appear  at  many  locations  in  the  flow.   Figure  5.30  shows  color 
contour  plots  of  autospectral  line  values  taken  from  the  measurements  presented  in 
Section  5.2.4.   The  four  plots  are  as  follows: 

(a)  The  312.5  Hz  line  is  representative  of  the  dominant  low  frequencies  characteristic  of 
the  integral  longitudinal  scale  of  the  flow.  Within  the  primary  potential  core  up  to 
x/2b~3,  there  is  very  little  large-scale  influence.  Outside  of  this  region,  the  pressure 
level  is  almost  uniformly  high  except  for  two  locations  at  x/2b~10  on  either  side  of 
the  centerline  where  there  are  significant  local  minimums.  The  symmetry  of  these 
minimums  suggest  a  standing-wave  pattern. 

(b)  The  3.625  kHz  line  is  the  dominant  line  of  the  screech  noise  corresponding  to 
turbulence  convected  through  shock  waves.  Using  Eq.  (2.45),  the  predicted  screech 
frequency  is  fe=2.8  kHz  which  is  close  to  the  observed  value.  Alternatively,  using  the 
observed  frequency  value,  the  Mach  number  predicted  from  Eq.  (2.45)  is  1.6.  In  light 
of  the  fact  that  Eq.  (2.45)  was  derived  for  discharge  into  a  quiescent  fluid,  the  current 
observation  is  reasonable.  The  observed  frequency  is  plotted  from  locations  where  it 
can  be  clearly  identified.  The  core  region  is  clearly  seen  as  having  a  very  low  value 
until  x/2b~7  where  a  transverse  standing  wave  pattern  is  located.  Assuming 
symmetry,  there  are  4  peaks  and  5  troughs  across  the  mixing-section,  which 
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corresponds  to  a  wavelength  of  2b,  the  height  of  the  primary  nozzle-exit  and  the 
characteristic  length  for  the  screech  frequency  calculation. 

(c)  The  plot  in  (c)  shows  the  49.625  kHz  line  where  an  acoustic  or  broadband  peak  was 
present.  This  frequency  is  that  of  the  acoustic  noise  radiated  from  the  region  starting 
at  the  trailing-edge  boundary-layers  of  the  nozzle,  to  the  shear-layer  transition-to- 
turbulence.  Spherical  waves  of  this  type  of  radiation  have  been  observed  by  Lowson 
and  Ollerhead  (1968).  The  highest  value  appears  close  to  the  wall  starting  at  x/2b~3, 
which  is  on  a  line  at  the  Mach  angle  of  the  primary  flow.  The  peak  is  then  smeared 
along  the  wall  as  the  radiation  can  be  convected  at  a  range  of  Mach  numbers,  from 
Ml  down  depending  on  the  location  in  the  freestream  or  boundary-layer.  From  this 
result,  it  is  now  possible  to  explain  the  peak  at  half  this  frequency  that  is  observed  in 
many  of  the  autospectra  (cf  Fig.  5.17).  Curie  (1955)  showed  that  a  resultant  dipole 
field  at  the  solid  boundaries  is  the  limit  of  Lighthill's  quadrupole  distribution,  and  the 
dipole  field  will  oscillate  at  half  the  frequency  of  the  quadrupole  field  as  shown  in 
Chapter  2. 

(d)  This  plot  shows  the  region  from  plot  (c)  where  only  an  acoustic  peak  was  present 
(x/2b;s7). 

It  should  be  noted  that  consideration  of  transverse  duct  modes  were  also  investigated 
for  the  mixing  section.  The  lowest  frequency  mode  is  predicted  to  be  the  first  mode  in 
the  y-direction  with  a  frequency  of  -3.4  kHz.  Although  this  is  close  to  the  observed 
3.625  kHz  frequency,  it  was  rejected  as  being  the  reason  for  the  observation  because  the 
3.625  kHz  frequency  was  observed  to  be  in  anti-phase  on  either  side  of  the  jet  (cf  Fig. 
5.21),  in  accordance  with  the  screech  tones  predicted  by  Powell  (1953a). 
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Table  5.2         Calculation  of  the  probability  W  that  the  wall  pressure  p  exceeds  the  threshold  Pt^ 
at  (x/2b,  y/2b)=(6.89,  0) 


Location  '     -* 


i 
0 
1 
2 
3 
4 


-f''\w{-P->-^\ 


(6.89, 0)  a 


yPrms         Prf 

% 


48.4 
10.0 
0.39 

0 

0 


49 
13.5 
2.7 
0.6 
0.2 


(6.89,  0) 


I  Prms        Prt 


~nns        *rrm$  I 


51.6 
12.6 
3.5 
1.0 
0.36 


51 
13.5 
2.3 
0.4 
0.1 


(6.89,  0)  a 


I  Prms         PrmM 


100 
22.6 
3.9 
1.0 
0.36 


100 

27 

5.0 

1.0 

0.3 


Gaussian 


100 

32 

4.6 

0.27 

0.01 


a=Schewe  (1983). 
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Figure  5.2        Skewness  at  the  mixing- 


•section  inlet,  x/2b=0. 
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Figure  5.3        Kurtosis  at  the  mixing-section  inlet,  x/2b=0. 
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Figure  5.4       Probability  density  function  of  wall-pressure  at  the  mixing-section  inlet, 
x/2b=0. 
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Figure  5.5        Turbulent  wall-pressure  along  the  wall  centerline,  y/2b=0. 
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Figure  5.6       Skewness  along  the  wall  centerline,  y/2b=0. 
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Figure  5.7       Kurtosis  along  the  wall  centerline,  y/2b=0. 
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Figure  5.8       Turbulent  wall-pressure  along  the  wall  at  y/2b=-1.64. 
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Figure  5.9        Skewness  along  the  wall  at  y/2b=-1.64. 
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Figure  5.10      Kurtosis  along  the  wall  at  y/2b=-1.64. 
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Figure  5.11      Comparison  of  turbulent  wall-pressure  along  the  wall  at  y/2b=0  and 
y/2b=-1.64. 
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Figure  5.12      Comparison  of  skewness  along  the  wall  at  y/2b=0  and  y/2b=-1.64. 
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Figure  5.13      Comparison  of  kurtosis  along  the  wail  at  y/2b=0  and  y/2b=-1.64. 
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Figure  5.16     Probability  density  fbnction  of  wall-pressure  at  the  location  of 
maximum  skewness  and  kurtosis,  (x/2b,  y/2b)=(0,  6.89). 
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Figure  5.18      Wall-pressure  spectral  decay  slopes  at  the  mixing-section  inlet,  x/2b=0. 
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Figure  5.19      Development  of  the  wall-pressure  spectra  along  the  centerline,  y/2b=0. 
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Figure  5.19--continued 
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Figure  5 . 1 9— continued. 
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Figure  5.20      Development  of  the  wall-pressure  spectra  along  the  line  y/2b=-1.64. 
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Figure  5.20~continued. 
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Figure  5.21      Two-point,  transverse  pressure-spectra  at  x/2b=0.   Measurement 

locations  are  in  the  upper  and  lower  secondary  streams  at  y/2b=-1.36 
andy/2b=1.09. 
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Figure  5.22      Two-point,  transverse  pressure-spectra  at  x/2b=0.   Measurement 
locations  are  in  the  primary  flow  at  y/2b=0  and  y/2b=-0.27. 
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Figure  5.23      Two-point,  longitudinal  pressure-spectra  at  y/2b=0.   Measurement 

locations  are  just  before  the  end  of  the  primary  core  at  x/2b=8.83  and 
x/2b=9.15.   Measurement  bandwidth  is  20  kHz. 


164 


200 

150 

100 

50 

0 

-50 

-100  - 

-150  - 


-200 


10-' F 


10' 


*   ■ '  t  t  I  I  I 


10 
t(Hz) 


10 


d 


10 


10 


10' 


1:  ISTiHit—M 
5:  3ia  kHt  ^ire' 


lo' 

f,(Hz) 


10 


Figure  5.24      Two-point,  transverse  pressure-spectra  just  downstream  of  the  end  of  the 
primary  core  at  x/2b=10.34.   Measurement  locations  are  either  side  of 
the  centerhne  at  y/2b=-0.27  and  y/2b=0.27. 
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Figure  5.25      Two-point,  transverse  pressure-spectra  just  downstream  of  the  end  of  the 
primary  core  at  x/2b=10.34.   Measurement  locations  are  on  the  lower 
side  of  the  mixing  section  at  y/2b=-1.08  and  x/2b=-0.27. 
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Figure  5.26      Two-point,  transverse  pressure-spectra  at  the  streamwise  location  where 
the  shear  layers  interact  with  the  wall,  x/2b=17.23.   Measurement 
locations  are  on  the  lower  side  of  the  mixing  section  at  y/2b=-0.54  and 
x/2b=-0.27. 
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Figure  5.27      Two-point,  transverse  pressure-spectra  at  the  stream  wise  location  where 
the  shear  layers  interact  with  the  wall,  x/2b=17.23.   Measurement 
locations  are  on  the  lower  side  of  the  mixing  section  at  y/2b=-1.87  and 
x/2b=-0.81. 
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Figure  5.28      Two-point,  transverse  pressure-spectra  at  the  streamwise  location  where 
the  shear  layers  interact  with  the  wall,  x/2b=17.23.   Measurement 
locations  are  on  the  lower  side  of  the  mixing  section  at  y/2b=-1.87  and 
x/2b=-0.27. 
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Figure  5.29     Two-point,  transverse  pressure-spectra  at  the  streamwise  location  where 
the  shear  layers  interact  with  the  wall,  x/2b=17.23.   Measurement 
locations  are  either  side  of  the  centerline  at  y/2b=-1.87  and  x/2b=1.08. 
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CHAPTER  6 
CONCLUSIONS  AND  RECOMMENDATIONS 


A  detailed  experimental  investigation  of  a  two-dimensional,  Mach  1.8  air-primary, 
Mach  0.3  air-secondary,  ejector  has  been  successfully  carried  out.  The  investigation 
showed  that  it  was  feasible  to  measure  and  interpret  the  wall  pressure  in  fine  detail  and 
utilize  the  fast,  nonintrusive  technique  for  identifying  and  locating  the  various  regions  of 
the  ejector  flow.  Additionally,  it  was  shown  that  mean  impact-pressure  and  Schlieren 
optical  measurements  could  be  utilized  to  confirm  shear  layer  growth-rate  and  entrainment 
details  of  the  flow.  It  was  also  found  as  a  result  of  the  measurements,  that  theoretical 
pressure  spectrum  functions  that  have  not  been  previously  confirmed,  do  indeed  exist. 
Turbulent  pressure  statistics,  one-  and  two-point  pressure  measurements,  the  mean 
flowfield,  shear-layer  laminar-to-turbulent  transition,  and  shadowgraph  visualization  were 
also  employed  to  establish  more  definition  of  the  flowfield. 

6.1  Conclusions  from  Significant  Results 
(a)  A  whole-field  plot  of  p^^  clearly  shows  the  primary  potential  core.  Additionally,  the 
edges  of  the  shear-layer  can  be  observed  in  region  1 .  The  structure  of  the  field  is 
qualitatively  comparable  to  the  mean  impact-pressure  flowfield,  and  the  lengths  of 
regions  1  and  2  are  observed  to  be  the  same  as  measured  from  the  impact-pressure 
data.  Along  with  the  conclusions  drawn  in  (f),  (g)  and  (h),  it  can  be  claimed  that 
time-accurate  wall-pressure  measurements  provide  a  fast,  nonintrusive  technique  for 
characterizing  the  various  regions  of  ejector  flow,  and  identifying  significant  flow 
events. 
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,.,  XHe  .™pac..p.ssu.  an.  v.ua,  .e.-,a,e.  .ow.h  n,.es  .„  ..o.  ,  of  .Ke  — d 
,ec.o.  no«  a.  *e  sa.e  as  p..c.e<.  .0™  previous  .eas.e.e„.  .n  compress* 
,eeshea,->aye.as,o„sas*epn.a.pote„.,a,-co.e„ds>.fore,orwhen,theshea, 

,ay.  .aches  .he  .aU.  A>so,  an  e^o.  function  prov.des  a  good  approx„na.,on      the 
impact-pressure  profile  in  the  shear-layers. 
(.)  -The  dtstance  to  transitton,  x.  was  estimated  .on,  .he  SchUeren  v,sua„za..o„s  .0 
occu.  a.  «.,^7.5X1C.  and  «....6.1.,  wh.ch  .s  tn  a,ee.en.  w.th  p.ev.ous 
nieasuvemente  of  compressible  free  shear-layers. 
(,)  The  developmen.  length  .  self-s,n„lan.y  for  .he  .rnpac.  pressure  was  found  to  be 

^.«.5,  also  .n  a^eentent  w.ft  nteasnrenren.  of  compressible  free  shear-layers, 
(e)  A.4e,n.e.plane>eprin,aryflowhasskewnessof-0.25a„db.r.os,sof-.5.  The 

values  ind,ca.e«.eoccurrenceofn,orepos,.,vepressure-peaU*annega.,vepressure- 
pealcs  wh,ch  ,s  opposUe  .0  that  of  .he  reporred  l.terature  for  incompress.ble  flow. 
(,  .  regron  1  a.  5...b.S.  the  shear-layer  edges,  deten.,ned  fr^m  unpac.  pressu« 
measurements,  show  «  .he  shear-layer,  low  speed  boundanes  f,a„en  ou,  a.d 
,,^-speed  boundaries  .ncrease  m  slope.    However,  .he  shear-layer  grow^-rates 
.mam  approximately  constant.    Th,s  .s  co.eCured  .0  resuU  from  .ncreased  loca 
en.rammen.  at  .e  low-speed  boundary,  and  by  v.scous  effects  a.  the    .gh-spee 
boundary.  Th.s  phenomenon  .s  .so  obse^ed  tn  the  color  contour  plot  of  — . 
wall-pressure  as  a  d.tinc.  loca.  tncrease  m  p„  emana.ing  from  *e  cen.er  of  the  shear 

observed  locally  around  the  centerUne  of  the  mixing  seCion. 
(,)  Bowns^eam  of  region  ,  bufbefore  .he  end  of  reg,o„  2,  large  edd.es  are  observed  on 
.he  m.ins-sect.on  cen.rUne  m  .e  .mpaC-pressure  flowfreld  *a.  are  drfferen.  from 
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b        d  m  free  shear-layers^  Also  obsewed  in  .his  region,  is  a  sudden  increase 

rianlL^e  .ow-s.ed  .undar.  of.e  shear-iayer,  as  weU  as  s,^.«ca^.  ^ 
.d  troughs  .n  the  .rbuienr  pressure  aions  the  center.,ne...,sco,ec.redthUhese 

lJonsarearesu,tofthenecessar.nowre,u.re™enttosat,s.theo.r      n^s 
.d  voiu^etnc  entra.nnrent  ratros,  wh.ch  .„  thrs  case  are  not  e,ua,  to  the  >oca  s^ 

,,er  en.ra,nn,ent  ratros  .n  regron  , .  Using  the  deta.  of  the  ■«— -" 
Its  for  reference,  .t.sohv,ousthattheco.orcontour,otofp„aUod,sp,aysth.s 

same  flow  phenomenon. 
,H)Tra„sve.e,  two-po.t  spectra,  nreasurements  sho.  a  P—  ^  ;;-^-; 
fre,uenc,  pea.  a.  the  end  of  regron  >  and  where  the  shear-,a,er  reaches  the  «a,.. 
le  pi:  correspond  to  edd.es  and  the  iarge-scaie  ...ng  process  ohse.ed  a^ 
„,J,us,onCf,,aswenas.o.hereg.ona,ongthe.aUat.<.2h.Z.show„h,.he 

tapact-pressure  measurements  and  p„  color  contour  plot. 
(,  The  sUewness  and  .urtos.s  of  the  pressure  exh«  large  confrdence-interval  pea^  at 

Lnrsttwo  locations  houndthe  sudden  entra.nmentarea,nreg.on.next,s.heend 

of  the  pnnta.  potent.-cor.  and  the  las.  ,s  the  star,  of  shear-layer  .n.eract,on  .. 

f  i^<;<7  S  and-2^K:$ll,  andall 
the  wall.  The  large  confidence  intervals  range  from  -1  .S.2.5  and 

f  «  ivith  fairlv  narrow  coniidence 
of  these  locations  exhibit  local-maxtmums  of  p_  w*  fa.rly 

intervals. 
0)Poradd.onal,nsi^..ntothephys,csof.heahove-men.,o„edflowphenonre„a^- 

found  that  pressure  autospectra  could  he  rel.ahly  calculated.  Kesults  show  *a.  the 
Lospec  Jchange  s.gn.ficantly  m  the  vanous  regrons  of  the  now,  and  clearly  show 
3pectral  slopes  of.'.  *-.  W."'^  in  drfferen.  regions,  representing  the  boundary- 
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.ay.  ove.ap  res-on,  *e  .u.bu,ence-»rbu,e„ce  ,„.e.c.,o„,  and  *e  .^uienc^ean- 

shear  interaction,  respectively. 
«„.sbeUevea.haU.sMhenrste.pen™e„t.wWch.he.'"^apec.a,s,opehasbee„ 

obse^ed.  perhaps  as  a  result  of  the  h.gh  Reynolds  number. 
0)  The  autospecra.  measurement  show  tha.  *e  ejeCor  configurafon  provdes  a  flow 
which  conrams  all  .he  currently  known  Aeoretica.  pressure-spectral  ranges. 

There  are  at  leas,  two  areas  of  study  .hat  faU  ou.  of  *,s  research  .ha.  would  benefit 
,om  further  .nveshgation.  These  areas  are  the  strucre  of  .he  mean  flowfteld  and  .he 

Structure  of  the  turbulent  pressure  field. 

u  «  finwfield   it  is  suggested  that  impact- 

To  more  completely  charactenze  the  mean  flowfield,  it  is      gg 

p_e  traverses  be  made  ,n  the  .d,rect.„  m  several  y.planes.  T^ee  d_l 
v.suaUzat,o„  would  then  reveal  .he  mean  structure.  To  vahdate  .he  use  of  .he  rmp  C- 
pressure  probe,  se.ecttve  use  of  a  nonmtrus.ve  dlagnostrc  techn.ue  such  as  LBV  or  a 
pla„arlmasmgteehm,ue.srecommendedPollow,nsthesesteps,.t.srecom— . 

omputatrona,  flu.d  dynam.c  (CPB)  calculatron  and  ca,.bra..on  for  .he  .eCor   ou.  . 
emblduponw,.hah,shcer..n.ofsuccess.  U,sof„o.etopo,n.ou.thats  ould*e 

.e  .mpacpressure  measurement  techn.ue  be  well  val.dated  for  .Ms  m»ns    ow  t 
,  ,  a  very  cr,.,cal  measuremen.  for  CFD  vaUdatlon  and  calibratton.  especally 
heterogeneous  and/or  non.soenerset.cm,xing-.aye,s.Th,s,s  because  the  impact-pressure 

.  a  «on  of  S...C  pressure,  temperature,  molecular  we.ght,  and  constant-pressure 

Cet^or  conf.surat.on  prov.des  a  How  wh.ch  contains  .1  .e  currently  .own 
theoretrcal  pressure-spectral  ranges.    T.,s  ,s  novel  and  ..  ,s  sugges.ed  that  de.,M 
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these  ranges  be  instigated.  It  is  also  recommended  that  conditional  averaging  of  the 
turbulent  wall-pressure,  and/or  velocity  in  the  near-wall  region,  be  investigated  to  study 
the  details  of  turbulent  bursting  in  compressible  shear-flow. 

Other  possible  considerations  for  experimental  verification  and  expansion  are 
variations  of  independent  experimental  parameters  including:  primary  Mach  number, 
secondary  Mach  number,  convective  Mach  number,  entrainment  ratios,  mixing-section 
inlet  area-ratio,  dissimilar  gases,  and  nonisoenergetic  flows.  Variation  of  the  physical 
configuration  to  an  axisymmetric  and/or  variable-area  mixing  section  is  also 
recommended. 


APPENDIX  A 
BOUNDARY-LAYER  MODELS 


The  boundary  layer  momentum-thicknesses  were  calculated  using  the  compressible 
Thwaites  method  for  momentum  thickness  of  Rott  and  Crabtree  (1952),  that  is  valid  for 
the  compressible  boundary  under  arbitrary  pressure  gradient  with  no  imminent  separation. 
The  necessary  assumptions  are  :  (a)  adiabatic  wall,  (b)  the  Prandtl  number  of  the  gas  is 
unity,  and  (c)  the  viscosity  ^  is  a  linear  function  of  temperature. 

From  assumptions  (a)  and  (b),  the  temperature  distribution  inside  the  layer  is  given 
by  the  Crocco-Busemann  relationship: 

I^-U^Mi\-}^\  (A.1) 


T  2 


U 


that  describes  a  smoothly  varying  temperature  distribution  having  values  of  static 
temperature  at  the  freestream  boundary  and  stagnation  temperature  at  the  wall.  The  static 
temperature  is  given  by 

Il=ul±M^  (A-2) 


According  to  the  Ulingworth-Stewartson  transformation,  new  coordinates  x^,  y;  are 
introduced 
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3T-1 


-litj'"'" 


(A.3) 


-KPlf^ 


To]       {  ^' 


dz 


(A.4) 


where,  for  a  given  x,  the  pressure  is  constant  across  the  layer.   The  velocity  distribution 
outside  the  boundary  layer  must  also  be  transformed,  by  putting 

(A.5) 


The  same  relationship  holds  inside  the  boundaiy  layer  between  the  values  of  u,  and  u: 


u.l^.,  z,.)  = 


^    u(x,z) 


(A.6) 


From  Eqs.  A.2,  A.3,  and  A.5,  it  follows  that 


dU; 


dx.    \T) 


^T.pdU 
dx 


(A.7) 


The  compressible  momentum  thickness  is  given  by 


H'-T, 


U  \   p  U 


U)  pU 


dz 


(A.8) 


and  in  the  transformed  case,  this  becomes 
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U- 

1- 

V 

fj 

U; 

-!-dz. 

u.    ' 


(A.9) 


From  Eqs.  A.6,  A.7  and  A.4,  these  two  momentum  thicknesses  are  related  as  follows: 

■'       2(T-1)  Q 

(A.  10) 


9. 


T 

\0i 


The  transformed  properties  can  now  be  evaluated  using  the  incompressible  method  due 
to  Thwaites  (1949): 


&.!=^^oUr'furdx^ 


(All) 


where  a=0.45,  and  b=6.    If  dxj,  Uj,  and  0;,  are  replaced  by  Eqs.  A.3,  A.5,  and  A.  10, 
respectively,  the  result  is 


0  2=av„ 

X  0 


f^oPHl.-/^f^ot^(^L... 


[T 


u-f 


[T, 


V'dx 


(A  12) 


or  for  the  case  when  y=\A, 


0/=O.45Vo 


[T) 


in        2 


^liTi 


U^dx 


(A.  13) 


This  result  permits  the  direct  computation  of  the  compressible  case. 
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The  displacement  and  boundary  layer  thicknesses  are  calculated  from  the  following 
equations,  respectively,  as  given  by  White  (1974,  p.607)  for  an  insulated  plate  and  unity 
Prandtl  number: 


S*  1 


(A.  14) 


2«  «    2as/^     [/nS-V^J 


1-1 


(A.15) 


where 


a 


^m)? 


(A.16) 


APPENDIX  B 

FLOW  PARAMETERS  FOR  TABLES  4.1  AND  4.2  AND 

SUPERSONIC  PITOT  PROBE  DATA  REDUCTION 


The  Mach  number  is  defined  from  the  ideal,  one-dimensional  isentropic,  compressible 
flow  analysis: 


Af= 


(p  \1± 


P) 


1 


(B.l) 


N  Y-1 

where  Pq  is  the  stagnation  pressure,  P  is  the  static  pressure  and  y  is  the  gas  specific 

heat  ratio. 

The  ideal  mass  flowrate  is  given  by 


m=A 


Y  ^0 

M 

^Ao 

2 

2(Y-1) 

(B.2) 


where  Tq  is  the  stagnation  temperature,  A  is  the  inlet  area  of  interest  and  R  is  the  gas 
constant.   The  mass  entrainment  ratio  is  defined  as 


W- 


w. 


(B.3) 


The  gas  constant  is  given  by 


R^"' 


MW 


(B.4) 


181 


182 


where  MW  is  the  molecular  weight  and  R^  is  the  universal  gas  constant,  where  the 
universal  gas  constant  is  defined  as 

/ 


R  =8314.39 


The  static  temperature  is  given  by 


°K-kg-mole 


(B.5) 


T=- 


I        v-1        -,\ 

I       2 


(B.6) 


The  sound  speed  and  average  sound-speed  for  an  ideal  gas  are  ,  respectively: 


c=s/yRf 


(B.7) 


^1+^2 


The  stagnation  and  static  densities  are  found  from  the  ideal  gas  equation: 


(B.8) 


Po  = 


RT^ 


(B.9) 


The  average  density  is  defined  as: 


RT 


P1^P2 


(B.IO) 


(B.ll) 


The  stream  dynamic  pressure  is  defined  as 
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(B.12) 


(B.13) 


The  volumetric  flowrate  is  calculated  from 

P 
and  the  volumetric  entraimnent  ratio  is  defined  as 

E  =—  . 
"    V^  (B.14) 

The  velocity  is  found  from  the  Mach  number  and  sound  speed: 

U=Mc  .  (B15) 

The  average  velocity  is  defined: 


u    =— — - 

a^  2 


(B.16) 


The  convective  velocity  is  defined  as: 


l+is  (B.17) 

The  primary  and  secondary  convective  Mach  numbers  are  found  fi-om  their  respective 
definitions: 

(B.18) 


K,' 

^1     ^c 

^1 

^.r 

u^-u. 

c. 

(B.19) 
A  single  convective  Mach  number  for  the  flow  is  sometimes  used  and  is  defined  as 

(B.20) 


ave 
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Mach  number  can  also  be  defined: 

M=-^ (B.21) 

c 

ave 

^„e«c  mode.  ««.  the  Unnard-Jones  .2-6  poten.a.  ft.nct,c„: 

A  rK,  ic  the  value  of  the  collision 
where  a  ,s  *e  hard-sphere  diameter  m  Angstroms  and  Qv  ,s  the 

integral  given  by 


A        C        E 

ilv~-^-^^-7i:  •  (B.23) 


.    1  ^f.^A^  R=0  14874  C=0.52487,  D=0.77320, 
The  values  of  the  constants  are  A=1.16145,  B  0.14874, 

E=2.16178  and  F=2.43787.   T.  is  given  by 


T    ^        o3<r<ioo  _^^. 


Th^  values  of  hard  sphere  diameter  and 
and  e^/k  is  the  characteristic  gas  temperature.   The  values 

characteristic  temperature  for  air  are  3.711  Angstroms  and  78.6  K. 

The  average  molecular  dynamic  viscosity  is  defined  as 


/^aw=-Y-  •  (B.25) 


The  dimensional  Reynolds  number  is  defined  as 


/A 


(B.26) 


the  average  dimensional  Reynolds  number  is  defined  as 
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Re 


I      _  "  ave     a 


f^a 


(B.27) 


and  the  jet  Reynolds  number  is  defined  as 

pi^i  m 


RV- 


/^i  (B.28) 

The  theoretical  incompressible  two-dimensional  mixing-layer  growth  rate 
(Dimotakis,  1986)  is 


5'o=^ 


1-r  ^ 


,l*r^, 


U^- 


\-is 


(B.29) 


1 +2.9(1 +r)/(l-r) 
where  e  is  a  constant. 

For  calculations  of  stagnation  pressure  from  probe  pressure  when  the  local  flow  is 
supersonic,  the  Rayleigh  supersonic  pitot  formula  is  used: 


s,  probe 


2        j       [yA         Y-IJ 


(B.30) 


where 


pit 


s,  probe 


,    2    J 


(B.31) 


Ppi,  is  the  stagnation  pressure  measured  by  the  probe  (i.e.,  behind  a  normal  shock)  and 
Ps,  iHobe  is  the  wall  static  pressure  at  the  streamwise  location  of  the  probe  tip. 


APPENDIX  C 
PRESSURE  TRANSDUCER  PARAMETERS 


Table  C.l        Micro  Switch  pressure  transducers,  model  number  142PC30D 


SERIAL 

No. 

RANGE 
(AkPa) 

CALIBRATION 

CURVE 

(AkPa) 

LE^ARITY/ 
REPEATABILITY 

(%FSO) 

USE 

4100020 

0-207 

32.61  *V-40.29 

±0.75/±0.15 

Pi 

4960005 

0-207 

32.68*V-43.82 

±0.75/±0.15 

P2 

4960012 

0-207 

32.64*V-39.51 

±0.75/±0.15 

P5 

4960020 

0-207 

32.68*V-40.97 

±0.75/±0.15 

P3 

4960025 

0-207 

32.82  *V-40.70 

±0.75/±0.15 

P4 

9430003 

0-207 

32.39*V-39.80 

±0.75/±0.15 

P6 

9430013 

0-207 

32.83 +V-40.63 

±0.75/±0.15 

P9 

9430020 

0-207 

32.50*V-39.86 

±0.75/±0.15 

P,i 

9820004 

0-207 

32.67*V-39.31 

±0.75/±0.15 

Ps 

9820006 

0-207 

32.39*V-40.04 

±0.75/±0.15 

P,2 

9820011 

0-207 

32.64*V-38.28 

±0.75/±0.15 

P7 

186 


187 


Table  C.2     Other  pressure  transducers  utilized. 


PRESSURE 
TRANSD. 

RANGE 

CALroRATION 
CURVE 

ACCURACY 

(%BFSL) 

USE 

Omega 
PX181 

101.3-2169 
kPa 

P^b+206.8*V 

±0.3 

Poi 

Validyne 
DPI  5 

±.8274 
AkPa 

P,„,+10.34*V 

±0.5 

P02 

Validyne 
DP15 

±12.41 
AkPa 

P.„b+1.241*V 

±0.5 

DP2 

Validyne 
DP15 

±689.5 
AkPa 

P^b+41.37*V 

±0.5 

p 

^  O.probe 

Validyne 
DP15 

±344.7 
AkPa 

P,n,H+10.34*V 

±0.5 

p 

^  s,probe 

Validyne 
DPI  5 

±68.95 
AkPa 

P,„,+6.895*V 

±0.5 

Ps2,bk 
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